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1 Introduction 



In this work we are concerned with quantum electrodynamics (QED) at finite temperatures. 
The main result is that this theory is infrared (IR) divergent already in the cross sections 
to fourth order, a result in contradiction to the literature where calculations are carried out 



for all temperatures (||l[, g, |3| and |J - At the end of Sec. |3.4.4| we give a detailed 
review of the literature and point out where the shortcomings lie in those calculations. 

If we restrict ourselves to small temperatures k B T <C m c ~ we reproduce the IR-finitc 
results found in the literature (e.g. [[ll)) ), but in a much more straightforward way. 

We also present a rigorous discussion of the thermal corrections to the magnetic moment 
of the electron in this low temperature approximation and reproduce the value of [p^| . After 
having pointed out these results we now present some introductory remarks on the methods 
that will be applied. 

For the field theoretic aspect of QED at finite temperatures we utilize the causal approach 
to quantum field theory (QFT) because this provides us with a framework which is free of 
ultraviolet (UV) divergencies and allows a stringent investigation of the IR-behaviour (pi), 
p^|). We will present an introduction to this theory in Sec. where we will especially 
emphasize how the UV-divergencies and subsequent regularization and renormalization are 
avoided by carefully treating the involved mathematical objects as what they are, i.e. dis- 
tributions and not functions: the origin of the UV-divergencies in the standard framework 
lies in ill-defined products of distributions. 

In addition, the use of the causal method avoids the so-called pinched singularities which 
arise through products of ^-distributions. They stem from the Feynman-rules using thermal 
propagators. The causal method avoids such ill-defined products by construction. 

To include finite temperature we use the ideas of thermo field dynamics (TFD). The 



basics of this theory (p3fc, \24j, |2§) will be presented in Sec. \L2. 

In Sec. |^ we calculate all second order and some third order graphs we will need to 
investigate the IR-behaviour. The second order self energy contribution will be used to 
discuss the thermal corrections to the electron mass. There we reproduce the results of 
the literature (fiof, @, etc.). The vertex graph will be the basis for our calculation of the 
corrections to the magnetic moment /i c - of the electron. These corrections are unfortunately 
too small to be tested by current experimental techniques. The only experiments we are 
aware of that perhaps could detect finite temperature field theoretic effects in the not too 
distant future are measurements of the Casimir force (see J2(J, p7j and references therein) 
and of energy level shifts in highly excitated atoms (Rydberg states), see [^8) (and more 
recent work on this topic). 

Finally we present all terms that have to be considered for the IR-behaviour of the theory 



in the adiabatic limit (which will be explained in Sec. 2.1.1) and investigate if their divergent 



parts cancel or not. This is done in Sec. 3.4.1 and 3.4.2 with the before-mentioned result 



A short discussion of this result and the relation to the literature can be found in Sec. 



3.4.4. To conclude, we present some ideas that possibly could help to formulate an IR-finitc 
thermal theory in Sec. ||. 

Throughout this paper we use units with h = 1 = c if not mentioned otherwise and de- 
fine the Fourier transform as follows: 

f(p)-=J^ Jd A xf{xyv\ (1.1) 

To keep this in mind is important for the correct factors (2-7r). A remark on the notation: 
we often simplify the notation by the definition m := m c - and thus write for the electron 
mass either m c - or just m. 
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2 Quantum Field Theory at Finite Temperatures 



In this section we will present the formalism we are going to use to do QFT at finite 
temperatures. For the field theoretic aspect we utilize the so-called causal approach to 
QFT. To incorporate the temperature we make use of TFD. Both these formalisms will be 
introduced and explained in the following. 

2.1 The Causal Approach to QFT 

The traditional formalism of QFT is plagued with the well-known UV-divergencies one 
copes with by various regularization- and renormalization-procedures. The reason for these 
divergencies lies basically in the fact that one performs mathematically ill-defined operations 
by using Feynman rules or performing the usual time-ordering procedure that leads to them, 
respectively. There distributions (Feynman propagators or 0-functions and products of field 
operators, resp.) are multiplied as if they were ordinary functions, a procedure, which is in 
general not defined. This entails the appearance of the above-mentioned divergencies, which 
one has to get rid of by regularization and renormalization of physical parameters to extract 
useful information from the theory. 

Instead of causing divergencies by wrong manipulations and having to cope with them 
afterwards, a better strategy is to avoid them from the very beginning, carefully paying 
attention to the mathematics that underlies the theory. 

This has been done by Epstein and Glaser in the early seventies plj and then further 
developed by Scharf et al. 

This so-called causal (since causality plays a crucial role in the formalism we will soon 
present) or finite (because there do not arise any UV-divergencies) QFT is a perturbative S- 
Matrix theory. Using some basic physical assumptions as causality and Poincare-invariance 
the S'-Matrix is constructed inductively order by order. In order to do that the basic fields 
of the theory (described by well-defined free fields on Fock-space), their (anti-)commutation 
relations and how they couple in first order (which corresponds to the interaction Lagrangian 
in the standard formalism) have to be known. These free fields and the S'-Matrix are 
operator-valued distributions. To extract physical information the relevant objects (e.g. S- 
matrix elements) have to be smeared out with test functions <?, suitably falling off at infinity 
(c.g.[] g 6 6>(M)). These test functions cut off the long range part of the interactions. To 
end up with physically meaningful results the so called adiabatic limit g — > 1, which restores 
the total interaction has to be considered. This limit may not exist (since 1 ^ <S(M)). 

If g € <S(M), S(g) is well defined and free of UV- or IR-divergencies. The latter can 
show up in the limit g — > 1, e.g. in QED (because of the massless photon, which implies 
a long-range interaction like the Coulomb-Potential) for S-matrix elements. But even there 
(in QED at zero temperature) the adiabatic limit exists if the right physically measurable 
quantities (inclusive cross sections) are considered, making QED at T = a physically 
meaningful theory (in the cross sections at least up to fourth order - up to this level the 
existence and uniqueness of the adiabatic limit has been proven, see p2| ). 

This last remark alludes to another aspect of the causal approach. One has to be modest; 
many results known in the standard theory are not yet rigorously verified. For example, it is 
difficult to prove non-perturbative statements. Likewise a rigorous proof of the cancellation 
of IR-singularities to all orders a la Grammar- Yennie [^9| is not yet known - to mention just 
two things. On the other side, by calculating on sound mathematical grounds one always 
knows exactly what is going on - and there is no other possibility to achieve this knowledge. 



1 M denotes the Minkowski space, i.e. the manifold R 4 furnished with the Lorentz scalar product based 
on the metric tensor = diag(l, —1, —1, —1). 



5 



2.1.1 The Causal Construction of the S-Matrix 



Here we want to present the causal formalism. For more details than we give here we refer 
the reader to [§§, |o| and |l). 

In quantum mechanical scattering theory the S*-Matrix can be written under certain 
assumptions as a time-ordered exponential of the interaction Lagrangian (cf. e.g. Sec. 0.3. 
in ^2|). This and the work of Stiickelberg and Bogoljubov, Shirkov^j motivated Epstein 
and Glaser to start in field theory with the following /oraal| power series in the coupling 
constant: 

00 1 f 

S(g) = 1 + ^ — / d A xi...d A x n T n (xi, ...,x n )g(xi)...g(x n ). (2.1) 
n=i n ' ■* 

We cannot make any statement about the convergence of this series - that is a main content 
of the expression formal. But every T n will turn out to be a well defined operator-valued 
distribution. After smearing it out with the test function g € <S(M), it is a well defined 
operator on Fock-space 0. From the construction we will describe in the following, it is clear 
that T n is of the form 

Tn (^1 3 ■ ■ • 7 x n ) — ^ ^ t n (x\ x n , . . . , x n — 1 x n ^j . (Dk {x\ , . . . , x n ) . , (2.2) 

k 

where : Ok ■ is a normally ordered product of free field operators (a Wick monomial), and 
the numerical^] distribution tjj depends only on relative coordinates as a consequence of 
translation invariance. 

The so-called n-point distributions or rt-point functions T n will be constructed inductively 
order by order from T\ which has to be known from the beginning and is given as a Wick- 
ordered product of free fields on Fock-space^], i.e. expressed by well-defined objects only. 

In standard field theory T n would be constructed as follows: 

T n (x\, ...,x n ) = T{T 1 (x 1 )...T 1 (x n )} 

: = e ( x °n(i) - x n(2)) ' •■• ' ®( x n(n-i) ~ x u( n )) T i( x n(i)) • •■■ ■ Ti(xn (n )), (2.3) 
n 

where T denotes the time-ordering operator, T\ the usual interaction Lagrangian and the 
sum runs over all n! permutations II of the n variables x±, ...,x n . This leads to the above- 
mentioned UV-divergencies; the reason is that the time-ordering by multiplication with 
O-functions and subsequently the Feynman rules are usually not well defined procedures. 
The causal construction, which we will present now, is essentially a method to do this 
time-ordering in a well-defined way. 

The expansion of the S'-Matrix in (2.1) can be formally inverted^]: 



S(g) 1 = 1 + ^ I d i xi...d i x n T n (x ll ...,x n )g(x 1 )...g(x n ). (2.4) 

n=l U ' J 



2 For a short historical introduction and references we refer to the introductory remarks in [ p^ . 

3 Throughout this paper we use the notation 'formal' to denote infinite series without having control of 
convergence, which nevertheless can be treated in the well-defined setting of 

4 Furthermore, instead of treating the formulae to come in the context of formal power series, here it is 
always possible to take the sums up to a convenie nt, la rg e eno ugh but finite value m £ N instead of oo 



would then be used as a definition, motivated by ( p.q ); (2.1C) etc. can then be derived by calculations with 
polynomials. 

5 Here 'numerical' means that d contains no field operators; but it does not mean that it is necessarily a 
scalar. 

6 E.g. for QED at T = 0: Ti(x) := ie :^(x)^^(x): A M (x). 

7 Some remark concerning the typography: the literature on the causal method usually denotes T n by 
T„. Since we already want to use the ' ~ '-symbol to denote a different object in the thermal field theory 
context in agreement with the literature of this latter area, we decided to change the notation in the causal 
context. 
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By definition (c.f. (2.1) and (2.4)), T n and T n are symmetric in x%, ...x n . For the disordered 
set of n points in M we write 



X = {xj GM I j = 1, ...n}. 
If we write S := 1 + T, then formally 



(2.5) 



S- 1 = (1 + T)- 1 = 1 + ]T(-T) r . 



r=l 



Now we collect all terms of the same order in the last sum and get 

n 

T n (X) = ^(-ly^T^CXi)...^.^), 

3 = 1 Pj 

where the second sum runs over all partitions of X into j disjoint subsets 



(2.6) 



(2.7) 



I = X 1 U...UX i , X h ?9, \X h \=n h , h£{l,...,j} 



(2. 



We want to emphasize that the products of distributions in ([2/71) are well-defined: they are 
direct products of distributions, because the factors have disjoint sets of arguments. 

Some other relations between the T n and T n , which for example will be used to obtain 
Q2.23| ) - ( pT25| ), can now be derived with the help of 1 = Sjg)S(g)- 1 (pTl| ), ( fO] ) and collection 
of all terms of the same order. For that we refer the reader to [|22[ . The consequences of 
unitarity, translation-invariance and Lorentz-covariance of the S'-Matrix for the properties 

of T n and T n can be found there as well. 

We now turn to the most important property, that is causality. Consider two test func- 
tions <7i, g2 G S(M), for which there is a space-like surface that separates their support. Then 
there exists a Lorentz-system, in which all points of supp(<?i) are earlier^ then all points of 
supp(<72)- In this case we write supp(gi) < supp^)- We now require for the S'-Matrix 



S(gi + .g 2 ) = S(g 2 )S(gi) Vgi,g 2 G 5(M) if supp(gi) < supp(g 2 )- 



(2.9) 



Combining this with (2.1) and (2.4), we get after permuting the integration variables in a 
suitable way (see (22)) 



-^n \&1 ? •*•) X n ) T m (xi , ■ ■ ■ , Xm^Tfi—m (x 



m+1 ) 



c„) for {xi,...,x m }>{x m+1 ,...,x n } and 



T n (jCl j . . . , Xn) — T m \ X\ j . . . j Xfji) T n—jYi {Xjxi-^x 3 . . . , X n ) fol {xi , . . . , X m j <C 3 ■ ■ ■ j j ■ 

(2.10) 

This property shows that the T„ are time-ordered products; especially if all the arguments 
Xj have different temporal components ar^, we can permute them to order them in time: 

X-t ^> Xn ^ 



> x® n . Repeated application of ( 2.10 ) then yields 
T n (xi, —,x n ) = Ti{xi) ■ ... ■ Ti(iE n ). 



(2.11) 



Now we turn to the general inductive step to construct T n on condition that we know 
all T m , 1 < m < n — 1, and that they have the properties mentioned above. Using (2.7) we 

%0 < 3/0 in this 



s Definition: Given x, y S 
system. 



x is earlier then y in a certain Lorentz-system 
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can calculate T m , 1 < m < n — 1. Then we define the following distributions 



A n (#1 7 ■ ■ ■ 7 ) ■ ^ ^ -^"ni (-^)^n-rti ) ; 

-^n (^1 ' ' ' ' 3 ) ■ — ^ ^ T n — ni (Y 7 X n ~)T 7ll (^) , 

where the sum runs over all partitions 

P 2 :{x 1 ,...,x n . 1 }=XUY, X^0, 
into two disjoint subsets with \X\ = m, \Y\ = n — l — m. We also define 

(^1 1 ■ ■ ■ 7 ^n) • -^n (*^ 1 ? ■ ■ ■ 7 *^7i) A n 0^1 J ■ ■ ■ 7 ^n) • 



(2.12) 
(2.13) 

(2.14) 
(2.15) 



Then we define two other distributions, just extending the sums in ( E.12D and ( |2.13| ) to all 
partitions P^, including the empty set X = 0: 



^-n 1 ) ■ ■ ■ ; ) • — ^ ^ -^ni (^j^n-ni (^ *^n) 



(2.16) 



1 ? ■ ■ ■ 7 *^ti ) • — ^ ^ ^n-ni ^n)^ni C^O 

The distributions A n , R n are not known to us since they contain the unknown T n . 
know the difference 

D n = R n ~ A n = R n — A n . 

Now we will show, and this is the crucial step of the causal construction, how R n 
can be determined separately, thus making possible the construction of T n by 

T n = A n — A n = R n — R n . 

In order to achieve that, we investigate the support properties of D ni A n and R n . 



done in 22 . The result is the following: Let 

V + (x) := {y e M \{y - xf > 0, y° > x } 
denote the closed forward cone of x and 

V~(x) := {y e M \{y - xf > 0, y° < x } 

its closed backward cone. The n-dimensional generalizations are 

Tt(x) :={(x 1 ,...,x n )eM n \x j eV ± (x)Vj = l,..., 

Then it can be proven inductively that Vn > 3 

supp(Z)„(xi, ...,x n -i,x n )) C T n (x n ) U T„ (x n ) 
supp(i?„(xi, ... ))cr„ T ( 
supp(^l„(xi, ...,X n -l,X„)) C T„ (x„). 



>}■ 



(2.17) 



But we 



(2.18) 
(or An) 

(2.19) 
This is 



(2.20) 

(2.21) 

(2.22) 

(2.23) 
(2.24) 
(2.25) 
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Therefore, we say that l D n has causal support '[]. A n and R n are called the 'advanced' and 
'retarded' distribution, respectively, because of their supports being located in the backward 
and forward light cone, respectively. It follows now from the equations ( 2. 23] ), ( p. 24 ) and 



( 2.25| ) that A„ or R n and therefore T„ (with ( 2.19] )) can be constructed if one knows D n 



and has split it into two parts according to these support properties. 

Because of the special role of x n in this procedure, T„ as calculated in ( 2.19| ) in general 
has to be symmctrizised with respect to the arguments. 

This splitting is the crucial point in the whole causal construction, and we will see that 
it has to be done very carefully and that just multiplying by 8-functions, which would be 
the most direct method, will in general fail. It is exactly this splitting-procedure that causes 
the UV-divergencies if done wrongly. We will present the splitting, its problems and how it 
is done correctly right after some other remarks. 

Up to now we described the general inductive st ep fr om u-1 to n. To start the induction 
we have to know T\ , to construct D-i according to ( 2.15| ) and to prove explicitly that £>2 has 



causal support. Then D2 has to be split into R2 and A2 according to the causal support 



properties. After that T 2 is constructed using ( 2.19 ). 

Since we will have to perform the adiabatic limit g — > 1 explicitly in our calculations, we 
present here how it is carried out in the most suitable way: we choose a fixed test function 
go G <S(M) with go(0) = 1 and consider the limit 

g (x) := g (ex), with e -> 0. (2.26) 

This so-called scaling limit describes a method to perform the adiabatic limit g — ► 1 explic- 
itly. Then we have in momentum spaced 

m = (2^ / d ' x ^ x y kX = (2^)2 / d 4 xg (ex)e ikx = ±g (?) (2.27) 
and in the limit 

limg(fe) = M^f*j = (27r) 2 <5(fc). (2.28) 

PROOF: This is a distributional equation. Therefore we consider g £ <S(M) as a regular distri- 
bution j £ 5(M)'. Then we have \/tp G 5(M) 

lim {g,ip) = lim (g,<p) = lim / d i kg(k)(p(k) = lim / dr kgo(ek)<p(k) = / d i k(p{k) = 

S-i-l 9->l 9^1J f^Oj J 

= (l,<f>) = (1 )¥ >) = (2tt) 2 (8, 9 ). qed. 
Here we used the theorem of Lebesgue in the fourth step. 

This leads to another useful formula: 

1 . fk 



e 4 Ve 



d^kgik) = / d 4 fc-j5o - = / d 4 k(2TT) 2 S(k) = (2tt) 2 . (2.29) 



In explicite calculations, which we will do mostly in momentum space, this is incorporated 
as follows. We are interested in expressions (e.g. a matrix element) of the form 

d 4 k 1 ...d 4 k n g(k 1 )...g{k n )f(k 1 , k n ), (2.30) 

where T(fci, k n ) usually contains some normally ordered field operators and numerical 
distributions (cf. |2.2| ). Using ( |2.27 ) this equals 



|d 4 fc 1 ...d 4 fc„5 (^)...3 (^)f(fc 1 ,...,fc„). (2.31) 



9 Since supp(D n ) lies in the n-dimensional generalization of the light cone with respect to x n , i.e. each 
Xi, 1 < i < n lies in the region that can causally influence or be influenced by x n . 

10 We caution the reader that we adopt the usual definition of the Fourier transformation here, differently 
from the corresponding formula 3.11.2 in tea]. 
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Now we perform the transformations fcj — > eki for i = 1, ...,n , and get 

d i k 1 ...d 4 k n g (k 1 )...g Q (k n )f(ek 1 , ...,ek n ). (2.32) 



In this way the adiabatic limit is usually performed: the variables ki, in which the whole 
expression is smeared out are replaced by eki, and the test function g is replaced by go. 
Then some integrations usually have to be performed, and the limit e — > has to be taken 
in the end. However, sometimes the adiabatic limit can be performed 'trivially' in some of 
the variables fcj. This means that all the involved expressions are of such a structure that 
e can be set zero at the very beginning 11 !, and the remaining fc^-integration leaves us with 



factors (27r) 2 (use ( 2.29 )). This amounts to setting the considered ki equal to zero in ( 2.3C ) 
and replacing the corresponding integration J d A kig{ki) by (27r) 2 . This will be used in many 
of the calculations in this work. Another often used trick is to set the involved ki equal to 
zero where this is allowed, although the adiabatic limit is not trivial in this variable. An 
example is (p — k\) with p 2 — m 2 : this can be replaced by (p) . 

2.1.2 Splitting of Causal Distributions 

As mentioned in the precceding section the crucial step of the causal construction of the 
ii-matrix is the splitting^] of D n into an advanced and a retarded part 

By construction, D n has the following form (for an example see ( |C.3|) where we calculate 
D 2 for QED at T>0) 

D n (xt, ...,x n ) = ^d k n (xx - x n , ...,£„_i - x n ) : O k {x\, ...,x n ) :, (2.33) 

k 

where : Ok ■ is a normally ordered product of external free field operators (a Wick monomial), 
and the numerical distribution d^ depends only on relative coordinates as a consequence of 
translation invariance. It is also the numerical distribution d\, which contains the informa- 
tion on the support properties of D n , that is, being concerned with the splitting problem, 
we can forget about the operator structure. 

Given a tempered distribution d e S'(M® n = M m ), m = An, with supp(d) C F^(0) U 
r„ (0), the splitting problem can be formulated as follows: 

Is it possible to find a pair of tempered distributions {r, a] C S' (M® n ) with supp(r) C 
r^(0), supp(a) C r~(0) and d = r - a 1 

The answer is 'yes' 0. This will be shown in the following. 

A remark on the uniqueness of the decomposition d = r — a: Assume there are two solutions 
of the splitting problem {ri,ai} and {r2,a2}, d = ri — ai = r 2 — a 2 ; then the difference 
n — r 2 = &i — a 2 has support {0} C M 8 " and therefore it is of the following form: 



£ c a V*S(x), x 6 M°», V" := ®X 'T^ ' {2M) 

| Q | OX l - OX 4n 

where a is a multi- index and x = {x\, ...,x± n ) S M® n . 

So the splitting problem has an ambiguity in the solution at '0' only. We will see soon 
that it is exactly the behaviour of the distribution d at this point '0', which is crucial for 

11 Mathematically this means that we can apply the theorem of Lebesgue. 

12 A detailed discussion of this splitting can be found in |S| and [pfj| . 

13 Even the more general problem of decomposing a not necessarily causal distribution into two parts with 
certain properties is solvable under certain assumptions; see e.g. 
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the splitting. Depending on how 'singular' d is in this region, the splitting can simply be 
done by multiplication with 0-functions - or not; in this latter case d has to be split more 
carefully and a multiplication with O would generate divergencies. 

To have a notion of how singular a causal distribution is near '0' we introduce the 'singular 
order' of a causal distribution, which is a rigorous definition of the power counting degree 
of UV-divergence one knows from standars method jj. 

Definition: The distribution d £ S'(R m ) has the quasi- asymptotic s0 d £ S'(R m ) 
at x = with respect to p £ C°([0, oo), [0, oo)) if the limit 

lim p(S)S m d(Sx) = d (x) ^ (2.35) 

exists in S'(R m ). 
The equivalent definition in momentum space reads as follows: 

Definition: The distribution d £ S'(W n ) has the quasi- asymptotics d £ S'(W m ) at 
p = oo with respect to p £ C°([0, oo), [0, oo)) if the limit 

lim p(S)(dU\cp(p)\ = (d ,<p) + (0,0) (2.36) 

exists V<p £ S(W n ). 

The condition 0' is necessary to have a unique definition of the term 'singular order' as 
we will see soo n0. The quasi-asymptotics probes the behaviour of d and d in the vicinity 
of x — and p = oo, respectively (cf.p2[). 

We derive by a scaling transformation that (cf. |22| ) 

3lo £ R : lim = a" Va £ (0, oo). (2.37) 

Thus we call p(S) the power-counting function and give the following 

Definition: The distribution d £ S'(W n ) is said to have singular order ui if it has a 
quasi-asymptotics do at x = or its Fourier transform d has a quasi-asymptotics do 
at p = oo, respectively, with power-counting function p(S) satisfying ( [2. 37 )^. 



For examples and some remarks we refer the reader to [g2| and [Q • 

This notion of singular order uj is crucial for the splitting of the distribution d = r — a, 
because it can be shown that for uj < the splitting can be done essentially by multiplication 
with 0-functions and turns out to be unique^, whereas for to > one has to proceed more 
carefully and ends up with some ambiguities of the form of terms with support in {0}. Now 
we turn to the splitting problem more explicitly: 

Case 1 : u> < 0. In this case we have 

p{5) -> oo for 6^0, (2.38) 



14 An exemplary calculation, which shows how the singular order determined by sloppy power counting on 
the one hand and by the rigorous procedure presented here on the other can differ, can be found in [134]. 
15 We use the notation introduced by the authors who introduced this concept: [pa]. 
16 See footnote 

17 If we do not require do ^ in (2.35), and the corresponding condition in (2.36), respectively, any u)' > ui 
could be chosen as singular order of d. An exception is u) = — oo: then do = is possible and allowed. To 
have a definition where this case is in 



corp oratc d as w ell the singular order has to be defined as an infimum 



of all the values for which the limits (2.35) and ( p.36[ ), respectively, - leaving out the condition 0' - exist. 
18 If uj(r) < u)(d) and oj(a) < ui{d) is required . 
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which implies 

d(*),^£)^<*^->0 for 5^0. (2.39) 

Now we choose an arbitrary but fixed vector v = (vi,...,v n ) 6 r n (0). Then we define a 
{An — l)-dimcnsional hyperplane by 

n 

H := {x e M® n \vx = Y^ v o x i = °}- ( 2 - 40 ) 

j'=i 

H splits the causal support of d: all products Vj Xj 3XC either > for x € T„ (0) or < for 
cc G r„ (0). In addition we choose a monotonous C°° -function xo with 

'0 for t < 

\u(/) <( s e [0,1) for <t < 1 . (2.41) 
1 for f > 1 

This function is sort of a smooth version of a 0-function. Now it can be shown (see p2|) 
that the limit 

(vx \ 
— )d(x)=:Q(vx)d(x)=:r(x) (2.42) 
o J 

exists, defining the multiplication of d by a 0-function unambiguously^ and independent 
of v. Thus a solution of the splitting problem is given by 

r with supp(r) C T^(O) and a :— r - d with supp(a) C T~ (0). (2.43) 

Many calculations in field theory are best done in momentum space, therefore we give a 
formula in the form of a dispersion integral for f (for details and derivation: p2| , po| ): 

Case 2 : > 0. Now we have 

P(S) 



for 5^0. (2.45) 



Proceeding as in Case 1 would result in ill-defined and UV-divergent expressions. Never- 
theless it is possible to derive (see 22 ) a dispersion integral]^] for f in a more complicated 
way: 

f (p)=±_L (°° dt-, ^0 r, Vper* (2.46) 



19 Because of the claim in footnote |18|; adding distributions with point-support in {0} would violate this 
claim. 

20 Some intuitive remarks concerning this integral: The factor 1 _t± i0 is the Fourier transform of a 0- 
function and as a convolution in momentum space reflects the product with © in x-space. A factor t -^ i0 in 
momentum space is sort of integration in ic-space (remember: differentiation in x-space becomes multipli- 
cation by p in momentum space; therefore a division can be seen as an integration). Singular order ui > 
means that in momentum space d ~ p"; if we divide this (u> + 1) times by p, we get — , which has w = — 1 

and can therefore be split t rivia lly - thus -4r becomes convoluted with the Fourier transform of 0. This may 
give some intuitive idea of (2At). To derive this formula corresponding considerations for the test functions, 
with which d is smeared out, have to be incorporated rigorously. 
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Here the subscript '0' refers to a special property of this solution of the splitting problem, 
namely 



2^0(^=0 = 0, \a\<oj. (2.47) 

This means that fo is in a certain way normalized at the origin in momentum space. There- 
fore this special solution 7~o is called the 'central splitting solution' with normalization point 
0. There is also a formula for a more general solution with an arbitrary normalization point 
q g M® ra , but it is not very handy. 



The fact that for u > the ambiguity (2.34) of the splitting solution cannot be removed 



and remains a priori undetermined is very important. In momentum space (2.34) takes the 
form of a polynomial: 

to 

c *P a - (2-48) 

|o|=0 

These terms are called normalization terms and the coefficients c a normalization constants; 
they have to be determined by further physical conditions as Lorentz covariance, PCT- 
invariance, gauge-invariance, etc. 

For more remarks on the dispersion integrals, how analytic continuation is used to gen- 
eralize ( 2.46| ) to arbitrary p G M® n , which consequences the existence of perhaps infinitely 



many normalization constants has for a theory etc., we refer to |22| and |B0| 
Concluding remark: In this section we have presented the basic formalism of the causal 
approach to QFT and pointed out the main features. The main point is that there are no 
UV-divergencies in the causal approach. This is related to the splitting of causal distribu- 
tions. For singular order ui < it can be done trivially with a unique solution, and the 
results that are obtained for T n are the same as if they would be calculated using ordinary 
Feynman rules (if T n does not contain any subgraph with u> > 0; an example of this is QED 
in the case of tree graphs). For ui > the trivial splitting yields UV-divergencies - just as 
they appear using Feynman rules, whereas with a more careful procedure as described above, 
we end up with UV-convergent and well defined qu antit ies. Bu t in t his case the solution is 
not unique, there remain 'local terms' of the form ( 2.34 ) resp. ( 2.48j ) to be determined. 



Furthermore, the causal approach provides us with a well understood formalism to tackle 
the Hi-problem, which amounts to inquiring into the adiabatic limit g — > 1. There is no 
need to introduce an unphysical photon mass or other regularization schemes. This will be 
crucial in our calculations for QED at T>0. In the literature this limit (or the equivalent 
procedure in the formalism used) is usually done at the very beginning of the calculations, 
whereas the correct proceeding is to do it at the end. This causes relevant differences in the 
results, as we will see later. 

The causal formalism has up to now been used to investigate several theories (as QED in 
[2§ and Q, Yang-Mills in [§6| - @, the standard model in [|| - @, quantum gravity in 
[45] - EjJ, lower dimensional theories in [Q, |5^], etc.) and fundamental questions (as 
the ghost structure in |32| , the convergence of the 5-Matrix in [|53| , the fundamental role of 
gauge invariancc in QFT inp^] and J55|, etc.). For an introduction to the theory we refer 
the reader to S,fl30|l and |3lT 



2.2 Field Theory at Finite Temperatures 

In this section we make some remarks on field theory at finite temperatures and present the 
formalism we will use in our calculations. 

The basic aim of thermal QFT is to describe particles that are embedded in some system 
that serves as a thermal background, e.g. an electron in a black body radiation background, 
within a quantum field theoretic setting. In this work we will only consider thermal QED, 
i.e. photons, electrons and positrons in a black body radiation background, if temperature 
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is high enough including effects of e _ e + -pairs generated in the bath. We will not consider 
the implication of our results for other things of great interest as the quark-gluon plasma, 
etc. 

There are several approaches to thermal QFT. A very complete review of these methods 
and their relation among each other is given in |2j| , where an extensive list of references up 
to 1987 can be found as well. 

Since we want to incorporate temperature in the causal approach to QFT we presented 
in the foregoing section, we decided to use the ideas of TFD, the latter being an operator 
formalism without path-integrals and by its structure allowing a quite straightforward com- 
bination with the former^]. We demand from TFD to provide us with field operators, their 
(anti-)commutators and T\ to be able to use the causal formalism. In the following we will 
develop the ideas of TFD as far as we need them to achieve this goal. Then we will take 
this thermal T±, in our case for QED at T> 0, and start to calculate several things using 
the causal method. We compare the results with standard field theoretic calculations in 
QED at T> and establish several important differences, which will be traced back to some 
shortcomings in the standard calculations. 

We want to emphasize that these differences between our results and the literature do 
not depend on the fact that we use the causal method. The causal method differs from 
the standard calculation in the UV-region. But the problems in thermal QFT arise from 
the IR-behaviour, which can be treated using ordinary Feynman rules with an IR-regulator, 
since the involved graphs all have singular order lu < 0, and graphs with (lo > 0)-subgraphs 
have a very transparent structure (at least the ones we are concerned with) - cf. the remark 
at the end of the introduction to Sec. ||. One just needs a formulation of the theory wherein 
the IR-problem, say the adiabatic limit, can be treated in a rigorous and controlled manner, 
which is the case in the causal method. 



2.2.1 Thermo Field Dynamics 

We will present some short introduction to the ideas of TFD. This will be done in a rather 
loose manner, without any proofs or considerations concerning mathematical rigor, because 
we just need some basic features of TFD in an axiomatic way to build up a causal thermal 
QFT. All the things we will present concerning the algebraic background, etc. just serve 
to put the ideas we finally take over into some broader context. This part may be skipped 
without big consequences for the understanding of the calculations to come, since we will 



state the framework of our causal theory for T> in a rigorous way in section 2.2.2, quite 
independently of the motivating remarks that follow now. For more detailed introductions 
to TFD we refer the reader to ||, || and ||. 

A basic idea of TFD is as follows: We write the statistical mechanical expectation value 
of some quantity A in the mixed state p as an expectation value on a Hilbert space as in 
ordinary quantum mechanics: 

p(A) = (Q p \A\n p ). (2.49) 

More concrete: let {|n)} C TC be an orthonormal basis of eigenstates of a Hamiltonian on a 
Hilbert space Ti.: 

H\n)=E n \n), (m\n) = 5 mn . (2.50) 



21 Another reason to use TFD is the fact that it seems to be extendable to describe non-shell quasi-particles 
and non-equilibrium situations in a quite straightforward way (Jstj], Jj^]). In pj) . Sec. 2.5.6, it is shown 
that TFD and other formalisms are equivalent in some aspects as for example the perturbative calculation 
of thermal Green functions. 
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Define a state £H®H byQ 

|fi/J>:= J-$>-^|n>®|n) with ^:=tr(e-^)=^e-^, /3:=-^. (2.51) 

P n n B 

Then we have VA e £(W, H) 



(Qp\A»I\n p ) = =: (A). (2.52) 



(2.52) can be considered in an algebraic setting: Compare quantum mechanics in the 



standard vs. the algebraicPj formulation: 

standard: algebraic: 

state: ^> € H linear functional ui 

observable: A € C(H,H) < ► on a C* -algebra of observables A 

expectation value: v(A). 

Given a C* -Algebra A and a state ui on A, we can construct the following objects by the 
so-called GNS-construction (see |p8|): 

A representation 7r w of A on a Hilbert space 7i u , 
wherein 30 € 7i w : f2 is cyclic, and 

w(,4) = (0M„4)ft) VA e A 

If we now take an equilibrium state in the standard formulation, e.g. the grand-canonical 
ensemble, then we can define a linear functional on the algebra of observables A as follows: 

e -PH 

cop(A) := tr(pA) V A e A with p := t ^ e _ pH) ■ (2.53) 

Now we carry over the time evolution of the standard formalism to the algebraic one and 
get in this way a condition on wp restricting it to be an equilibrium state. This condition is 
the so-called KMS-condition (see Sj) and reads as follows: 

ujp(A t B) = up(BAt+i P ) with A t := e UH Ae- UH . (2.54) 

A state that obeys this condition is called KMS-state. 

If we now carry out the GNS-construction for a KMS-state, we get a formalism whose 
constituent objects have all the special properties necessary to apply the results of Tomita 
and Takesaki on modular algebras (TT-theory) (see |5^|, |^|, By that we finally 

have a formalism to describe states with temperature ^ 0, which formally can be built 
up like an ordinary Fock space quantum field theoretic one. serves as the vacuum 

state (the 'thermal vacuum'), and we can define and construct emission and absorption 
operators, field operators, propagators, etc. An important point is that the TT-theory 
for the GNS-construction on KMS-states demands a doubling of the degrees of freedom to 
correctly incorporate the time-evolution in the formalism. Roughly speaking one has the 
ordinary objects of the standard theory and in addition a copy of them, often referred to as 
the 'tilde-objects' (because they become marked by '~') and the time evolution is given by 
means of H — H. 

Historically the whole story developed quite differently. The doubling of fields in TFD 
was established mostly on more or less stringent physical considerations (see e.g [p5| ). Beside 

22 The 'doubling of degrees of freedom' that shows up here by H (g> H will be explained later. 
23 See e.g. [Rq|, fA. 
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that the doubling was necessary to get rid of certain singularities (the so-called 'pinched 
singularities', see e.g. p3| |, p. 205) that showed up in the theory in two- and more loop 
graphs as a consequence of products of (5-distributions (as we will see later, temperature is 
incorporated in the Feynman propagators by terms ~ S(p 2 — to 2 ) - thus ill-defined products 
of <5's are generated using ordinary Feynman rules). This doubling could also be the solution 
of the problem the author of [|50| sees with the existence of a perturbative thermal theory in 
general. Another argument for doubling the degrees of freedom is that a similar structure 
shows up in some other formalisms for thermal QFT. 

The connection between TFD and the algebraic results mentioned above has been estab- 
lished by Ojima in p4j| . 

2.2.2 TFD in the Causal Approach 

In this subsection we will build up the formalism to incorporate temperature in the causal 
QFT. This will give us a Fock space formalism as it is known from ordinary (T = 0)-QFT, 
we just will have two copies of every field, etc., as motivated in the foregoing section. 
The interpretation of the physical content of these 'tilde-fields' is somewhat unclear and 
controversial; in our opinion the calculations have to guide us in this respect. Thus we will 
for example see that the formalism urges us to interpret (^(k, (3) + Qp\ as a final state with 
one photon of momentum k absorbed from the thermal vacuum - so the argument that the 
'tilde-objects' are as unphysical as ghosts, as it is sometimes claimed in literature (e.g. in 
f6li), cannot provide us with the correct point of view. We will also see that 'tilde-fields' can 
show up in external legs of graphs. All these things will become clear in the course of our 
concrete calculations in section [| Some considerations concerning the physical significance 
of the 'tilde-objects' can be found in |25| . 
Now we will define our formalism for spinor QED at T>0: 

BOSONS: For the photons we take two copies of absorption and emission operators: 

o"(j?,/3), a AI (fc,/3) + and a»(je,p), a"(fc,/3)+ (2.55) 

which obey the commutation relations: 

[a" (k, (3), a" (¥, (3)+} = 5{k-k')5^ = [h> 1 (k,(3),a l/ (k / ,(3) + }, = otherwise. (2.56) 

The common vacuum shall be denoted by \Clp): (k , (3)\Q p) — and a M (fc, /3)\£lp) = 0. 
Now we define the field operators: 




A»(x) := —t / -^=1 \A + f(k)a^k,P)e~ lkx + VfQcja^k, P)e ikx 



±^/m^(k ) /3)+e- ikx ± \A + /(fcK(fc,/?) + e ifc * 




A»(x) := 7T ^ r / y/1 + f(k)a»(k, [3y kx + v7(fcK(£ Py- tk! 

h\k\ 



±^J(kja»(k,{3) + e lkx ± ^/l + f(k)^(k,(3)+e- lkx j. (2.57) 

Here f(k) :— e)3|fc ^ | _ 1 is the Bose distribution function, u € M denotes a time like normed 
(i.e. u M w M = 1) vector that describes the four-velocity of the background heat bathp] 

24 Incorporating this vector u in the theory, we establish a covariant formalism. Cf. the remark at the end 
of this section. 



25 According to the ideas of TFD as presented in Sec. 2.2.1 , the Fock space vacuum may be interpreted 



as a thermal equilibrium heath bath of temperature T = -r^- 
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(thus we have u — (1,0) in the rest frame of the bath). Furthermore, we define (3 := 
fc^yr, where k B is the Boltzmann constant (which can be set to 1 in suitable units) and T 
denotes the temperature (as a Lorentz scalarQ. From '±' the '+ '-sign is to be taken for 
the components^] /i = 1, 2, 3, '— ' for /i = 0. We emphasize that in ( [2.57 ) we implicitly have 
k 2 = and fco = |fc|. Incorporating this covariantly ( 2.57 ) reads as follows: 



A»(x) =— - / d 4 kS(k 2 )Q(k Q ){ yjl + f(k)a»(k,/3)e- ikx + v7M«"(M)e ,k 

, (2.58) 



± ^J{k)^(k,P) + e- ikx ± y/l + fWa^k,®- 



Commutation relations are given in App. [X| 

Fermions: As for the photons we take two copies of absorption and emission operators 
for the electrons: 

b,(f,0), b s (p,P) + and b.fafl, b s (p,l3) + , (2.59) 
and the same for the positrons: 

d.%0), d s (p,(3) + and d s (p,P), d s (p, P)+ . (2.60) 

they shall obey the following anti-commutation relations: 

{b s (p,0),b a ,(j?,P) + } = 6(p-i?)6 ss , ={b s (p,P),~b s >(p',P) + }, =0 otherwise, (2.61) 
{d a (p,0),d,(p,0)+} = 8{p-p')8 ss , = {d 3 (p,f3),d s (jS,/3)+}, = otherwise. (2.62) 

Again we define field operators: 
(2tt)2 I 



+ VI - f+(p)d s (p,p) + v s AP)e lpx + tVUip)d s (p,p)v s Ap)e tpx j, 

:=— ^ V f d 3 p\ y/T^TW s (p,/3)u s , a (f)^ +i^jAp)b s {p,P) + UsAP)z wx + 
(2tt)2 <-f J I 



v/l-/+(pK(p, 0) + v s , a {f)e- ipx - tVU(pjd s (p, ^v^e-*" 



(2.63) 



where f±(p) '■= el3(|p „| L ^ A , )+1 is the Fermi-Dirac distribution function with chemical potential 
±/Lt, and u, v are the ordinary four-spinors, solution of the Dirac equation, as used in the 
(T = 0)-theory. The anti-commutation relations for the field operators can be found in App. 

In the definition of the field operators and in the following we adopt the viewpoint of 
describing the electrons, positrons, photons and their tilde-conjugates by 6 different sets 

26 For the different possibilities how the temperature could transform under Lorentz transformations see 
|p2| and references therein and recent developments in this area, respectively. For the present calculations we 
adopt the most popular choice: the temperature is a Lorentz scalar. We did not investigate the consequences 
these other possible choices for the transformation properties of the temperature would have for thermal 
QFT. 

27 Compare with the quantization of the radiation field in VA, Sec. 2.11. 
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of emission and absorption operators in the same Fock space (with vacuum as it is 

ordinarily done in standard QED, with its three different 'particles'. 



Here we can already see from (2.57) that the incorporation of the temperature will cause 

0. 

In TFD the 



additional IR-problcms since f(k) diverges as for |fc| 



Now we want to give some motivation for the formulae (2.57) and (2.63) 
following relations between the ordinary emission and absorption operators of the (T = 0) 
theory and the ones for T>0 can be established: 



o"(M) 



ch sh\ / a^(fc) 
sh ch) U^(fc) + 



a^(fc) 
5"(*)+, 



ch -sh\ / a"(Jfe,/3) 
-sh ch ) \a»(k,/3)+ 



(2.64) 



bs(p) 



b s {p) 
ib s {p)' 



— s 



b s (p,P) 
ibs{P,PY 



(2.65) 



Here (ch) 2 := 1 + f(k), (sh) 2 := f(k) and c 2 := 1 - /_(p), s 2 := /_(p). 'ch', 'sh' can 
be written as a cosinus hyperbolicus and sinus hyperbolicus, resp., of a certain angle, 'c', 
's' as a cosinus and sinus, respectively. Such transformations as in ( [2.64| ), ( 2.65| ) are called 
Bogoliubov transformations. For the implications of this connection between T = and 
T > 0, generalizations of them or further results related to them we refer the reader to the 
literature (see e.g. j56|, p5| ). We will not consider these things here because they are of no 
use to us. We just take these relations to give some more motivation for our formulation of 
the thermal QFT: Given ( [2.64D and ( |2.65D , they can be inserted into the formulae for the 
field operators at T = 0, which read as follows (see any QFT book; but we adopt here the 
normalization of [E2|, Sec. 2.2 and 2.11): 




(2.66) 



T=0( 



(27T)i 



E 



d 3 p{ b s {p)usAP)e- ipx + d s (p)+v s , a (p)e 



IpX 



(2.67) 



and we get ( 2.63 ) and an expression quite similar^j to ( 2.57 ). Pay attention to the fact that 
unlike in the (T = 0)-theory here the 'positive' and 'negative frequency part', resp.' (' e ±l,ca: ') 
is not anymore associated with the absorption and emission part, respectively. For the 'tilde 
absorption' and '- emission' operators it is just the other way round. This will evoke quite 
different structures compared to T = in some calculations. 

To use the causal method we now need to know the various contractions of the thermal 
field operators (because of the Wick theorem we will use) and T\ for QED at T > to 
start the whole procedure. Beside that, we give some other useful formulae as contractions 
in momentum space, causal distributions, Feynman propagators and some relations among 
them. All this, except for T±, which we want do discuss now, can be found in appendix 



At the end of section 



2.2.1 



we mentioned and shortly motivated the doubling of the 
degrees of freedom in thermal QFT, especially the new Hamiltonian H := H — H, which 

28 Since we want the '±' to be associated with the emission operator part of the field for T > as well 
(cf. the argumentation in |E2j Sec. 2.11) we have to adjust these signs in a suitable manner. 
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has to be used. Based on this consideration we define T x for QED at T>0 as follows^: 



Ti(x) : = ie :^(x)^^{x): A^x) - ie :*(x)7 m *(.t): A^x). 



(2.68) 



From this equation we see that there are two types of vertices in the theory, one with 
'ordinary' fields, denoted by 'type 1', the other with 'tilde- fields', denoted by 'type 2'. 

A remark to avoid confusion: In our formalism of QED at T > we only have thermal 
objec tf^\. We only have the thermal absorption and emission operators from (2.55), (2.5£) 
and ( 2.60D , and the field operators built up of them as in (2.57) and ( 2.63| ). Any mentioning 
of (T = 0)-quantities was only done to give some motivations for our formalism. Thus 
we emphasize especially that any normal ordering (e.g. in Ti) is done with respect to the 
absorption and emission operators from ( 2.5£ ),(2.59) and ( 2.6C| ), and that expectation values 
have to be considered between states generated from the thermal vacuum 

Remark: In this section we have constructed a covariant formalism for QFT at T>0, 
here especially for QED. We see that covariance need not be lost in a thermal QFT as often 
claimed because of the preferred rest frame distinguished by the heat bath. Using a covariant 
formulation one usually has a better control of the calculations^]. In appendix [b| we give 
some short (not rigorous) calculations in this formalism, which may serve to strengthen the 
confidence in it. 

We want to emphasize that different authors in TFD sometimes use different conventions, 
which differ by factors i, — 1 etc. in the formulae, e.g. for the Bogoljubov transformations. 



3 Cross Sections to Fourth Order for QED at T > 

After having set up the framework of our theory in the foregoing section, we want to calculate 
all second order terms of QED at T > 0, which will allow a discussion of the thermal mass 
correction (using the self energy of the electron). We also need them to calculate some third 
order contributions as the vertex, which will be needed to discuss the thermal corrections 
to the magnetic moment of the electron. Furthermore, we will gather all contributions from 
first to third order to the cross section (for scattering of an electron on a C-number potential) 
to fourth order to investigate the IR-behaviour of the theory. 

The main result is that thermal QED, as it is usually formulated, is IR-divergent already 
in cross sections to fourth order if it is not restricted to small temperatures k B T <C rn e - (see 
Sec. p~2[ ): There is no complete cancellation of IR-divergencies as it is claimed in literature 



(0, || , etc.). Besides that we get a much simpler cancellation of IR-divergencies than the 
one in the literature (|1Q), etc.) for small T; we do not have to introduce thermal spinors 



and a Dirac-equation for T>0 in this approximation (see Sec. 3.4.4) 



F urt hermore, we calculate the thermal corr ection s to the electron magnetic moment (see 



Sec. p.3| ) and the thermal mass correction (Sec. p.l.Sj ). We calculate these two quantities fx e - 
and m e - for k B T -C m e - only, since the theory is not yet good for arbitrary temperatures 
because of the beforementioned IR-divergencies, and it does not make sense to extract much 
information before we have a theory that works for all temperatures at our disposal. 

A rather complete version of the calculations can be found here and in the appendices 
since our results depend heavily on these detailed somewhat tedious derivations. So we 



- — T— 

29 This loosely could be written as 'TT~® — T\ ', but we want to emphasize that it would rather be 
Ti +Ti , using our definition for Ti of QED at T = 0, which we take over from , and the ordinary definition 
of a 'tilde-operation' as given in the literature, because of the factor i in our Ti (the 'tilde-operation' takes 
operators O to operators O and numbers c S C to c*, the complex conjugate; it relates ordinary and 'tilde-' 
objects. We will not define it in a rigorous way or discuss it because we will not need it). 

30 Beside the spinors u, v. These are spinors as used in the (T = 0)-theory. We emphasize this because 
of some calculations in the literature where so-ca lled 'thermal spinors' and a thermal Dirac equation are 
claimed to be necessary to do QED at T>0 (e.g. []l(J). This is not the case as we will see. 

31 Some people claim for example that th^ differences in the value of fi c - in literature originate from 
problems with non-covariant calculations ([[ 
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want to provide a reader, who would possibly like or is forced to reproduce the results with 
as much material as possible. The calculations to come are done by the causal method as 



described in section 2.1.1, starting with T\ from (2.68). Doing the calculations we will see 
explicitly that ordinary Feynman rules could be used when lu < — 1 (cf. the remark at the 
end of Sec. 2.1). This is for example the case in the second order calculations. In the third 
order calculations we see in addition that g raphs with subgraphs with lu > have to be 
treated more carefully - as described in Sec. [2.1.2 . Thus it is not the case as often claimed 
in literature that the thermal contributions cause no additional UV-problems. Surely, all 
the subgraphs with lu > are temperature independent, but in combination with the rest of 
the graph they give T-dependent contributions. So the essential structure of possibly new 
UV-problems is already known from the (T = 0)-theory, but how they contribute to matrix 
elements, etc. depends on the thermal structure as well. 



3.1 Second Order QED at T > 
3.1.1 Calculation of D 2 

According to ( |2.15|) , (|2.12|) , (|2.13|) and (|J) we have 
D 2 (x,y) = B^(x,y)-A' 2 (x,y) = T 1 (y)T 1 (x)-T 1 (x)T 1 (y) 



-Ti{y)T\{x) + Ti(x)Ti(y) 
-A' 



+ :^(x) 7 ^(x): A„{x) :<Z(yh v *(y): A v (y) 

- ^(xh^(x): A„{x) :%h M %): My) 

- ^(xh^(x): A„{x) A v (y) 

+ &(x)'f&(x): A„(x) (y): A„(y) 



(3.1) 



dito with x «-> y 



(3.2) 



+ lablcd ^a(x)^ h (x):^ c (y)^ d (y): A^(x)A u (y) 

-ia b ic d ■Mx)Mx)--*c(y)* d (yy. M*)Mv) 

-"fabled &a(x)h(x):& e (y)94y): A^x)A v {y) 
+ lablcd *.(4W»*c(#i(»): M*)Mv) 



dito with x 



+R' 



(3.3) 



Now we use _ C(**), where C(**) := {^~\ * (+) } denotes the contraction, 

and the general Wick theorem to getp| 

:§ Q tf 6 ::tf c * d : = * a * 6 # c * d - f a * b C(f c * d ) - C(^ a ^ b )^ c ^ d + C^ a ^ b )C(^ c ^ d ) = 
=:* a *6* c^d- -C(# a # c ) :^ b ^ d : +C(* a # e2 ) :* b * c : +C*(* b * c ) :§ a * d : + 
+ C(y b * d ) :# a * c : +C(<i> a t> d )C(t> b <i> c ) + C(# a * c )C(* 6 * d ). 

(3.4) 



2 For shortness we denote both ty a ,c and ty a ,c by ^ a ,c, ^b,d and *S>b,d by ^b,d in ( p.4[ ) etc. 
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Since C(**) = 0, C(**) = etc., we finally have 



-c(* a tf d )c(* 6 tf c ). 

(3.5) 



So we will need the following four equations: 

:* a (x)* b (x): :* c (tf)*«l(l/):=:*«(a!)*6(!c)*c(y)*d(y): + (y)} :^(x)* c (y): + 

+ {^(z), :* Q (x)* d (y): (s), ^(y)}^^), *<+>(y)} 

:* o (a0¥ 6 (z): :* c (y)f d (y):=:* a (x)* 6 (a;)§ c (t/)f d (j/): +{¥<-> (s), f d +) (y)} :*&(x)* c (y): + 
:* a (x)* b (x): :* c (i/)*d(w):=:*.(a:)*6(sc)*o(»)*d(»): *i +) (y)} ^i,(»#c(y): + 



*i +) (y)} *a(x)*«j(»): ^ +) (y)}{*6 *i +) (y)} 



?(-). 



:* a (x)* 6 (x): ^ c (y)^ d (y):=^ a (x)^ b (x)^ c (y)^ d (y): +{*(->(a;), * d +) (y)} :* 6 (x)* c (y): + 

+ {*l\x), ¥+\y)} :V a (x)h(y): +{# (">(*), * d + \y)}{^ (x), ¥+\y)}. 

(3.6) 

For the photons we need 

M*)MV) ='M*)MV)' +[A^(x),A^(y)], (3.7) 



and the corresponding relations for AA, AA, AA. Inserting fl3.6|) and ( |3.7| ) in (3.3) we finally 
get the expression (C.3), see App. C.l, for D 2 . 

Now we have to show that D 2l eq. ( |C.3| ), has causal support. This is done in App. |C.2[ 
Because the temperature independent parts of D2 have causal support separately (see J22[) 
the temperature dependent parts have to have causal support as well. Thus we can split the 
T-dependent and -independent parts of D 2 separately. The T-independent parts are treated 
in ^2|, so we can restrict ourselves to the investigation of the T-dependent ones. This will 
turn out to be a great advantage because for these parts the splitting can always be done 
trivially; they all have singular order uj < as we will see soon^J. 

3.1.2 Calculation of T 2 

In this section we calculate T 2 . We will do that separately for every different process that is 
'contained' in D 2 (distinguished by different operator structures multiplying the numerical 
distributions) . This is possible because each term that describes a certain process has causal 
support on its own (cf. the discussion of the causal support of D 2 on page[8C|). For the longer 
calculations we will refer to appendix C.3. A remark on the notation: we often handle the 
different parts of D 2 separately. The l Aku Afci' on the left indicates from which terms the 
expressions on the right origin (cf. eq. ( |C.3| )). We also want to emphasize that sometimes 
inessential factors or coefficients are omitted in intermediate steps (e.g. during the splitting 
procedure). Furthermore, we often take track of the (T = 0)-parts just to have complete 
formulae. These parts are treated exactly in |22 and are not involved here since the T- 
depend ent an d -independent parts can be treated separately as mentioned above at the end 
of Sec. |3.1.l|. 



33 This is true in the n-th order as well: As long as the temperature dependent part of D n contains no 
(T = 0)-subgraphs with ui > it can be split trivially. 



21 



The 'Trivial' Graph 

For An, -Bn, Cn, Dn we have D 2 = in all four cases; 



To 



i?2 — i?9 



An 



/3j-itrivial 



e 2 : f (i/)t"*(i/)¥ ( X yf*(x) :: ^(y^x) 
e 2 : ^(yW^iy^ixh^ix) :: ^(yji^x) 
e 2 : *(y)<f&( 3/ )*(x)<f*(x) :: i^y)^) 
e 2 : *(y) 7 1/ *(y)*(x) 7 ^*(a;) :: i,(y)i M (x) 



~i?9 



(3.8) 



Electron-Electron Scattering (M0ller Scattering) 



B 



12 



7-)M0llcr 
u 2 



A-12 + A 



Cfn + B 



Df 2 + £ 



12 



- e 2 : *(x) 7 A '*( a :)*(y) 7Al vI/(y) : j ( /»£)(+) (a, _ y ) _ ^> (y ~ x)) 

+ e 2 : §(x) 7 ^(^)*(y)7^(y) : » ( - 2/) - " *)) 

v . ' 

=o 

+ e 2 : *(x) 7 ^(x)*(y)7^(y) : i ( - If) " "-D& } (|/ " *)) 



=o 



e 2 : #(a0 7 ^(z)%)7/^(2/) : * ( - 2/) - - *)) 



'°D(y-x) 



(3.9) 



Since this Z?2 just contains the causal Distribution D, which appears in the (T = 0)-theory 
and which has singular order ui — —2, cf. [ p2| , the splitting can be done trivially. Thus we 
have 



T. 



M0llcr _ 



= R-B! 



A a A- A 13 



T=0 



>D(x - y) 



= (>Dg(x-y)- ?D[-\x-y)= ^(x-y), 
B? 2 + Cf 2 : - (- <>D£\v x)) = („ - x), 
C£ + Bf 2 : - (- - *)) = ^f 2 (y - x), 



12 

/3 
12 

^12 ' ^12 



D w( x - y) 



- {- P D ^\y - x)) = % 2 F 2 (x-y). 



With all the factors this gives finally: 



T, 



M0llcr 



A 12 


+ <: 


- ? e 2 


*(x)7^*(a;)*(y)7 M *(y) : 


D\i{x- 


y) 


B 12 


+ 


+ ie 2 


*(a;)7' t *(a;)#(y)7 M t(y) : 


D F 12 (y- 


X) 


°12 


+ Sf 2 : 


+ ie 2 


^(x)^^(x)^(y)j^{y) : 


dUv- 


x) 


U 12 


+ <: 


-*e 2 


^(x)-f"^(x)^(y)j^(y) : 


dU*~ 


y) 



(3.10) 



(3.11) 



22 



Compton Scattering 



2# = 



A 2l 



4" 



A 



31 



'■- : *{y)Y % x {y - x)^(x) 



B, 



21 



31 



e z - : ^{x) 1 ^S 11 {x-y) 1 v ^{y) 



M x )Mv) ■ 
M x )Mv) ■ 



+ e\ : i>{y)Y P S 21 {y ~ x) >f*(x) :: A„(x)A v (y) 



D, 



21 



H-.^ix)^ PS 21 (x-y)j»*(y) 



=o 



C 21 + B 



LJ 21 



H:*{x)Y % 2 (x-y)j^(y) 



h-.^iy)^ % 2 {y-x) 1 ^{x) 



=0 



A^x)A u {y) : 
A^{x)A v {y) : 



(3.12) 



^{y)Y % 2 (V - xW%(x) :: Mx)A v {y) : 

*(x)r %2(x-yh u hy) ■■■■ M*)My) ■ 



Since /3 5 , 11 / 12 ( z ) = T ~°S(z) (see App. A.2 ) and the latter has singular order oj = —2 (see 
p2]|) the splitting can be done trivially as well. For T 2 C this leads to (here we choose xq > yo 
without^] restriction of the general validity of the expressions): 



A 21 



4, 



(^ u ( y -z))-t- P S ( +\y-x) = 

- PSll{y ~ *) ~ S ( +\y -x) = ^(y - x), 

^(x-y), 



A 31 


+ < 






y) = 


DO 




:0 - ^(y-x) 


= ^l(y- 


x), 


°31 




:0 - P S { 2 -\x-y) 


— — £21 


-v), 


pa 
^31 




:0 - ^(x-y) 


= — 5f 2 (a: 


-y), 


L/ 21 


+ ^1 


:0 - ^+ } (y-x) 


= %(y- 


x), 


^21 






- p S [ +\y- 


x) = 



(3.13) 



D 



31 



/3 



- p S? 2 (y -x)- p s { 2 p(y -x)= esi 2 (y - x), 

( P S 22 {x - y)) ret - f>S&\x -y)=- P S F 22 {x y). 



34 Since we are dealing with distributions it does not make sense to consider their value for a certain 
argument. Thus we can avoid x = y. By Lorentz transformations of these covariant expressions we can 
avoid xq = yo as well. 
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With all the factors we finally have 



To = 



A 2l 


+ A-n : 


— «e 




A a 
A 31 


+ <: 


+ ie 2 


*(x) 7 ^(-l) ^(x - y) 7 ^(y) :: A M (x)^(y) 


^21 




+ ie 2 


*(y)7 V ^(y-x^ix) :: ^(^(y) : 


°31 


+ 


-ie 2 


^(^y (-1) "sUx - y )Y*(y) ■■■■ M*)My) 


r> a 
B 31 


+ <: 


+ ie 2 


*{x)r{-l) P S ¥ 12 {x - yh v *{y) :: A^{x)A v {y) 




+ 


-ie 2 


^{y)l v S! 2 (y-x)rV(x) :: A^x)A v (y) : 


T) a 
u 1\ 


+ 


-ie 2 


^{y)i v sUy-x)rhx) ■■■■ M*)Mv) ■■ 


u 3\ 


+ <: 


+ ie 2 


^J 2 (z - Mx)A v {y) 



Vacuum Polarization 



D 



VP 



41 



-B4I + C41 



e 2 tr 



^^{y-x^^ix-y) 



e 2 tr 



j^s[t\y-xh^s[^(x- y ) 

^^\ y -x) 7 ^S[t\x- y ) 



e 2 tr 



T^+^y-xh^^Or-y) 
^^\ y -x) 7 ^sH\x- y ) 



e 2 tr 



T^^+^y-xh^^^-y) 
r^iy-x^^siPix-y) 



j^SitHv-xh^S^ix-y) 



A„{x)A u (y) 



AJx)A»(y) 



A fJt {x)A v { y ) 



A»(x)A v ( y ) 



(3.14) 



(3.15) 



In appendix C.3 we calculate w VP with the following result: oj vp = —00 for the temperature 
dependent parts. Thus we can split them trivially and arrive at: T 2 = R — R' ~ (for the 
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terms Afr+A^): 

[(T = 0)-part] rct + T >°s[-\y - x)^ T=0 S$(x - y) + 

+ T =°Smy - x)-f T>0 S { n\x - y) + ^+ } (y - x)^ ^\x - y) = 
= [(T = 0)-part] rct + T =°S[+\y - x)-f T =°S^\x - y) + 

+ T >°s[-\y x)^ T =°Slf(x -y)+ T >°s[t\y - x)>f T =°S^(x - y) + 
+ x)-f T> °s[-\x -y)+ T -° S [V(y ~ x)-f T >°S^(x -y) + 

+ T >°S[t\y-x)r T>0 S ( n \x-y) = 
= [(T = 0)-part] - T>0 S F (y - x)^ ^Sf^x - y) - 

T =»sU v - x)-f T >°s F (* - y) - T>0 sUv x)>f T > s F (* - y). 

(3.16) 

For the other terms we get 

B2 1 + C& : D 2 = — > T 2 = -R' 2 ~ 

~ PS£\V - x)r S[ 2 \x -y)=- f3 Sf 2 (y - x)>f f3 Sf 2 (x y), 
Cft + : D 2 = — » T 2 = -# 2 ~ 

~ ^(y - x) 7 ^ - y) = - % F 2 (y - x) 7 ^ ^ F (x - y), 

+ Z?^ : calculation as for the A-part gives 

T 2 ~ [(T = 0)-part] - T>0 S 2 F 2 (y - x)^ T =°S F 22 (x - y) - 

- T =% F (y - x)-f T >°S F (* - y) - r >°S F (y - zh" T >°S F (x - y). 

With all factors we thus have 



(3.17) 



rrnVP 
J 2 



A% + Aa-i ■ 



+ e 2 tr 



7"([(T = 0)-part] - T>0 S F (y - x) 7 ^ ^"^(a - y) 



- r =% F (y - a )y T >%(, - y) - T >»SUy - x)-f T >°S F (x y)) 



■ M x )My) ■ 



C4I + A41 : 



I>?1 



41 • 



+ e 2 tr 



-e 2 tr 



+ e 2 tr 



y( - ^ F (y - :r)y % F (x - y)) : ^(^(y) : 
7"( - Sf 2 (y - x)r P S F l2 (x - y)) : A,{x)A v {y) : 
7"([(T = 0)-part] - T >°S F (y - x)^ T =% 2 (x - y) 



T=0oF 



] s! 2 (y x)>f T> % 2 {x y) T>0 s! 2 (y - x)-f T >% 2 {x y)) 



■A^AAy). 



(3.18) 
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Vacuum Graph 



dY g = 



+ 



B', 



C2 



Q 



42 • 



B 



42 • 



u i2 -f i/ 42 



ie 2 tr 



ie 2 tr 



+ ie z tv 



ie 2 tr 



7" s{-\y x) lv OsitHx - y) D^(x - y) ~ 

- Y P S[p (y - x) 7u Pstf (x - y) (y - x) 
7" ^{y - x) lv ts£\x - y) D{+\x - y) - 

- 1" P S ( +\y - x) lv ?s[ 2 \x - y) ^(j - x) 
7" P s[-\y - x)<y u ^{x - y) ^+\x - y) - 

- l v P S^\y - x) lv f>stf(x - y) l*D£\y x) 
7 " W-\y - x) lv f>S%\x - y) f>D&\x - y) - 

- l" ^i? (y - x) lu f>S&\x ~ V) 4V (V - x) 



(3.19) 



We use the relations from appendix A. 2 to rewrite the propagators in the D-term in a suitable 



way and thus prove D" 2 + D% 2 + A% 2 + A% 2 = 0. Analogously we prove B% 2 + C% 2 = and 
C42 + B 42 = °- This g ives 



\0 







T. 



VG 



-R'o 



A% 2 + + D" 2 + D^ 2 : 


+ ie 2 tr 




X)lv 


S { n 


\x 


-y) 


D{t\y- 


x) 




+ ie 2 tr 


7^2%- 


X)lv 


°22 


\x 


-y) 


D&(y- 


x) 


i?42 + : 


+ ie 2 tr 


\^sit\y- 


x)j u 


0q(~ 

°12 


\x 


-y) 




x) 


riot, 1 ~r>0 . 
°42 ' -°42 ■ 


+ ie 2 tr 




X)lv 




\x 


-y) 


D&(y- 


x) 



(3.20) 



In appendix C.3 we show that the adiabatic limit in the vacuum expectation value of the 
vacuum graph is of the form \c, with c = const. 7^ 0. Thus the adiabatic limit does not 
exist. But we have yet some undefined normalization constants at hand from the (lj = 4)- 
splitting in the (T = 0)-case (see |22|], 4.1.27): we can split different parts of a distribution 
separately if they have causal support each - to do that we each time use the singular order of 
the part considered. This provides us with the retarded and advanced parts of these different 
parts. But the normalization freedom is given by the maximal singular order of all these 
parts, since the philosophy is to split all parts with the same operator structure together - 
that is with the maximal singular order of these parts^]. If we now choose Co = = C2 and 
C4 = — c the limit exists. Since c = for T — (as can be seen from the expressions in App. 



C.3) the result for T — is not changed. 



35 Because of the remarks on the uniqueness of the retarded and advanced parts preceeding eq. ( |2.3^ ) this 
docs not contradict the statements above. 
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Self Energy 



u 1 



^32 ' A 22 



-not _i_ /~ip 
-°22 T ^32 



C32 + -^22 



^32 + ^22 



^22 + ^32 



-^22 + -^32 



-O32 + -^22 



-e 2 : ^{x)-y v 



-e 2 : *(i/)7* 



+e 2 : ^{x)-f u 



-e 2 : *(x)7 1 ' 



+e 2 : (j,) 7 " 



+e 2 : (j/h" 



-e 2 : ^(x)Y 



^\y-x)^\x-y) + 

+ ^\y-x)^\y-x) 



fi sit\x-y)^D[i\x-y) + 



+ f, sk\x-y)f>D[i\y-x) 
P S^\y-x)^+\x-y) + 

+ f> Sit\v-x)f > D[+\y-x) 
Ps£ ) {x-y)PD i £{x-y) + 

+ f ) S { -\x-y)^(y-x) 
^s[t ) (x-y)^D{p(x-y) + 

+ ^s[-\x-y)^+\y-x) 
f ) S { 12 \y~x)^(x-y) + 

+ ^S[p(y-x) D{+ ) (y-x) 
f3 S [ 22 \y-x) f, D$(x-y) + 

+ fl S&\y-x)< i D&\y-x) 
f j S { 2 t\x-y)^{x~y) + 

= ^-\x-y)PD<£\y-x) 



7„*(x) : 



i^{y) ■ 



i^{y) ■ 



iu^{y) ■ 



7**0*0 : 



(3.21) 



The calculation of the singular order of the self energy graph (see appendix |C.3| ) shows that 
lu se — — 00. Thus we can again split trivially and arrive after proceeding as in section^ at 
the following expression for T 2 = R — R': 



[{T = 0)-part] + T > 5f 1 (y - x) ^D^x - y) + 



+ T =°SUy - x) T>0 D F n (x -y) + ^S* u (y - x) ^D^x - y) 



T>0cF /■„, \ T>O n F 



(3.22) 



This are the terms from A 22 + A^ 2 . The other terms of A and D are calculated analogously 
and the terms of B and C are easy to calculate since D 2 = 0. Finally we arrive at 
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T, 



SE 



^22 + ^32 



+ e 2 : fo) 7 " 



{{T = 0)-part) + 1 = u ^ 1 (y - x) 1 >^{y - x) + 



+ 



T>OqF 



>S* n (y - x) ' =^(y -x)++ >«S* n (y - x) > >^(y - x) 



A32 + &v. 



22 



((T = 0)-part) 



T=OnF 



>S* u (x-y)^D* n (x-y) 



T>OqF 



, S? 1 {x-y)*=»D* u (x-y) 



T>OcF 



>S* n (x-y)'>»D* 1 (x-y) 



ivV(y) 



B« 2 + G 



C32 + -^22 



S32 + c. 



+ ^32 



*(j/)7" 



- Psl(x - y) ?D F 12 (y - x) 7 „*fo) 

- esUx-y) DUy-^iuHv) 



^ 2 (y-x)^Df 2 (y-x) 



((T = 0)-part) + T =°S 2 F 2 (y - x) J >u ^ 2 (y - x) + 



T>O n F 



+ J >»S% 2 (y - x) J =»D* 2 (y - x) + * >»S% 2 (y - x) 1 >»D* 2 (y x) 



D% 2 + D 



22 



*(x)7 i ' 



((T = 0)-part) 



T=0qF 



] S* 22 {x-y)^D* 22 {x-y) 



T>0qF 



^ 2 (x-y) ^"D^ix-y) 



T>0oF 



>SZ 2 (x-y) ^D^x-y) 



y„V(x) 



(3.23) 



Remark: Up to now we have calculated all contributions to T 2 . It can be verified 
quite straightforward that all these contributions are symmetric under the transformation 
x — > y,y — > x (use the relations from App. |A.2| ). We see that their temperature dependent 
parts are the same as if they had been calculated using ordinary Feynman rules (they all 
have lo < —I), whereas the (T = 0)-parts have to be treated more carefully as it was done in 
P2| . To get some insight in the consequences of these formulae, to extract physical quantities 
etc. it remains to perform the involved integrals. We will only do that in the case of the self 
energy under the assumption that fcgT <C m c - . This will allow the discussion of a possible 
thermal mass correction and is done in the next subsection. 

3.1.3 Self Energy and Thermal Corrections to m e - 

Here we want to discuss a possible thermal mass correction by the self energy. This is usually 
looked at as follows: we consider the infinite summation of self energy graphs 

S F (x — y) + J d i xid 4: x 2 S F (x — xi)Y,(xi — x 2 )S F (x 2 — y) + 

+ J d A x\d i X2d A x^d A XiS F {x — x\)Yi{x\ — x 2 )S F (x 2 — x^j^ix^, — X4)S F {x^ — y) + ... =: 
=: S tot (x-y). 

(3.24) 



2S 



This can be interpreted as sort of a total propagator. Taking the Fourier transform of this 
gives 



1 ;S tot (p) 



(3.25) 



(2tt)2 

= j^S F (p)+(27r) 2 S F (p)t(p)S F (p) + (2nfS F (p)t(p)S F (p)t(p)S F (p) + ... . 

Since this is a geometric series we can sum it and get 

<> F 1 1 

Q tot — — = (<l OR) 

1-(2tt) 4 S5 f (S-F)- 1 - (27r 4 )S (27r) 2 (^ - m + iO - (2tt) 2 E) ' K ' 

These calculations are not rigorous, e.g. we set g = 1 at the beginning. If we do it with 
j ^ 1, such a summation is impossible (just consider the Fourier transform of S F (xi — 
X2)g(xi)g(x 2 ) + S F (xi - 2/i)£(yi - y 2 )S F (y 2 - x 2 )g(x 1 )g(x 2 )g(yi)g(y2) + ••■)■ 

Another difficulty origins in the fact that we have two types of vertices. In the consid- 
erations here we investigated the sum of terms with non-tilde vertices only. For a complete 
discussion we should somehow incorporate the terms with tilde vertices as well or give ar- 
guments why this is not necessary. 



However, we will not discuss this and if ( 3.26 ) can be derived rigorously without per- 



forming the adiabatic limit too early - to do that we should first have an IR-finite theory at 



hand, and this is not the case now as already mentioned (see Sec. 3.4.2 ). So we just want 
to discuss some consequences of this formula to compare our results with the literature. 
First we see that (27r) 2 £ 'appears as a mass' in the propagator - this point of view was also 
fruitfully adopted in the (T — 0)-theory, e.g. to determine normalization constants in the 
causal approach (see |P2| , p. 213). But as we will see soon, S for T > has - beside the 
scalar component - terms proportional to j> and ji, which perhaps should be interpreted as 
a thermal change of momentum of the particle. This idea finally leads to the definition and 
consideration of thermal spinors and a thermal Dirac equation as introduced in J63| , |ic| . 

Here we just calculate £ and evaluate it especially for k B T <C m c -, which practically 
means that we can set the Fermi-Dirac distribution functions f± to zero (i.e. we disregard 
terms 0(e _/3m e- )). Then we derive a thermal mass, which is in agreement with the results 
in the literature (e.g. @, JTol). 



To achieve this goal we have to investigate the various terms of (3.23). There are four 



terms to be considered, the others can be derived from them (e.g. by using the relations 



from app.A.2) 



I : 


7" 




X) 


T>0 ^(y- 




II : 


7" 


T> %(?/- 


x) 




^)7m 


III : 






x) 


t>0 dUv- 




IV : 


7" 


T>Q sI 2 (y~ 


x) 


T>0 D F 12 (y- 





(3.27) 



A(x)B(x) — ► AB[p) = / d 4 kA(p - k)B(k), (3.28) 



First we turn to I. Using the formula for the Fourier transform of a product: 

1 

W) 

we get the following: 



d4 ^ jt + ^h 1 



(3.29) 



(2tt) 5 J {p-k) 2 -m 2 + iO ' e P\ku\_ 1 - 
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This has to be a covariant expression and therefore is of the following form: 

Co + Cp'J) + C u ji, with Co, p ,u Lorentz scalars, depending on p 2 and pu. (3.30) 



To evaluate ( 3.29| ) we use 7 M (^ — jt + m)^^ = — 2^+2^ + 4m, and the fact that Co origins from 
the part proportional to m, C p and C u from the one proportional to p^ + u^ (after omitting 
the 7^ from $ + ft). If we multiply these terms ~ (p^ + u M ) by p^ and u^, respectively, we 
get the following three equations: 



Co = 

C p {pu) + C u ■■ 
C p p 2 + C u (j>u) 







(2tt)5 , 




[d A k- 
i 




(27T) 5 . 


jd^k 



4m 



(p-kf 



Tfl 



i0 



-2pu + 2pk 



[p — k) 2 — m 2 + iO 
-2p 2 + 2pk 



S(k 2 
5{k 2 



1 



e /3\ku\ 
1 



(p — k) 2 — m 2 + iO 



-5(k' 



e 0\ku\ . 

1 



(3.31) 



The calculation of Co iPiU is done in App. |C.3| .4. The full result (in the approximation 
k B T <C m c - and including the factor e 2 ) on the mass shell p 2 = m 2 is 



C o =0 



p 48/3 2 tt 2 m 2 



(pu) 2 



C u 



A8(3 2 tt 2 m 2 — (pu) 2 



- 1 



pu 



pu 



2\J (pu) 2 — m 2 

™2 



In 



pu + \J (pu) 2 — m 2 



2^uf 



In 



pu - 

pu + y/ (pu) 2 



(3.32) 



pu 



V (pu) 2 



The terms II — IV of (3.27) are calculated in the same way, making an ansatz Cq + C p j/> + C u ji 
and determining these coefficients as above. However, in the approximation k B T -C m c - 
only the term I is relevant, therefore we will not investigate the other terms anymore. 

The result ( p. 32] ), after inclusion of the factor (27r) 2 from (27r) 2 E in (3.26), is in agreement 
with [JhJ (see their reference [16] for their calculations). How to extract a possible mass 
correction due to the temperature from this result is not that clear as often claimed. One 
way would be to proceed in the usual way (e.g. Q, |L0]]) by determination of the pole of the 
total propagator (3.26) for \p\ — > 0. This leads to the result known from literature, 



Amp 



12m 2 (3 2 



(3.33) 



Another viewpoint could be the more direct one of taking the term Co as a mass correc- 
tion and the 'slash'-terms Cpfi+C u fi as some sort of momentum change; then there would be 
no mass correction since Co = (see (3.32)). We emphasize that the contribution Cpfi + Cufk 
cannot be absorbed in normalization constants since they are not of a polynomial structure. 
We will not start to discuss this here any further since we take the viewpoint that we should 
first have a consistent thermal field theory for all temperatures. The determination of a 
mass and momentum correction certainly involves a discussion of the Dirac equation and 
the spinors in the thermal setting as it was done in ]63[ ]. As already mentioned it is not 
necessary to be concerned with this in the approximation k B T <C m e - - but it may has to 
be pursued to formulate a consistent thermal QFT for all temperatures (cf. Sec. ||). 



3.1.4 Matrix Elements for Bremsstrahlung 



To discuss the IR- limit in the cross sections to fourth order (see Sec. 3.4) we need the contri- 
butions from the one photon bremsstrahlung processes. Because they are given essentially 
by Compton-graphs we calculate the matrix elements already in this section on second order 
terms. The process is described by Compton graphs with an external classical 'non-tilde' 
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photon leg (a C-number potential) and an absorbed or emitted soft photon. According to the 
theory we have to build states using the emission operators in the thermal Fock space. We 
start with the most simple case and discuss general coherent and incoherent superpositions 
at the end of this paper. Thus we have the following matrix elements: 
Emission: 



SJ? = (bMP) + ^(k)^(k,0) + n\T^ oinptm \b a (q,P)+n), (3.34) 



Absorption: 



S$ s = (6 s (p,/3)+6^(-fc)^(fc,/3)+0|T 2 Compton |6 ff (g,/3)+f7). (3.35) 



We emphasize that there is no sum over the indices /z, therefore we write both of them as 
superscripts. A remark on e: Since physical photon states are transversal it is convenient 
to introduce real polarization vectors e M with e M = (0, e), k ■ e(k) = 0, e 2 = 1. Photon 
states are then given by |e A1 a A1 (/c) + f2). For convenience we choose e(— k) for the absorption. 
Further remarks on this can be found in plj around (3.5.9) and (2.11.32). To describe the 
emission of a photon by the final state^] (a M (fc, l3) + fl/3 | with initial state l^g) is quite natural. 
The other graph has to be considered since it also gives a contribution to the investigated 
process where we start with \£lp) as the initial state and end up with one ordinary or one 
tilde photon. These are the only two possibilities to have (0^| with a one particle change 
in the photon sector as the final state. Thus it is quite natural to interpret this second 
contribution as an absorption from the photon bath since this is the only physical process 
besides emission, which takes place and changes the heath bath by one particle in the photon 
sector. It was already argued in pEfl , p. 120, that the second process describes the generation 
of a 'hole', i.e. absorption of a photon. 

Now we consider the following contributions to BS: in x 2 we have a classical (thus non- 
tilde) external potential A v {x-i) c ^ and an incoming field 'J', in x\ an outgoing field and a 
field or iff and A M . Or we have A v (x 2 ) ext and an outgoing field ^ in x 2 , an incoming field 



and A^ or if and A^ in x-y. Thus there are four terms from T 2 ( 3.14 ) that contribute: 



Afi + A^ 1: + A21, Cfjj + and C 21 + B^. In momentum space, the contributions with 
X\ and X2 interchanged give the same result, but this resulting factor two is canceled by the 
factor 5T = 5 from the S- matrix expansion (^J). This gives 



S e f ™= Id 4 



x 1 d A x 2 g{x l )g{x 2 ) ^ l / 5 ip) u {-le 2 ^ ^S( x {x 2 - Xl )>f) ■ 
(2tt)2 



(2tt)s (2tt)* 
the three other terms. 



Mq)e- iqxi — t / dk S(k 2 )e(k )^/l + f(k)e lkxi e^k)Al xt (x 2 ) + 



(3.36) 



Here we used the covariant expression from 2.58 for A M (x) in order to keep track in a rigorous 
way of the conditions k 2 = and fco = The absorption part gives the same except that 
we have e~ lkxi instead of e lkxi and \ff instead of y/1 + f. Now we use the following Fourier 
transformed expressions: 



g(*i) = 7A2 / d%g{h)e-^\ S(x 2 - Xl ) = ^ J d 4 QS(Q)e^ x ^ , 



A(X 2 ) = l^rMr), 



(3.37) 



36 In the following discussion we omit the irrelevant polarization vectors and electron emission operators 
in the states. 
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then the ^-integrations give (27r) 4 <5(— k\ + Q — q + /c)(27r) 4 <5(— fc2 + P — Q — r), which allows 
to calculate the Q- and r- integrals (in the other three terms we proceed in an analogous 
way). In each term one ^-integration can be done trivially leading to a factor (2ir) 2 (cf. 
end of Sec. 2.1.1). With analogous calculations for the absorption we thus get the following 
matrix elements for emission and absorption: 



BS.cm 



(2rr)i 

% CXt / 



d%g{h) dk 6(k 2 )G(k )u s (p) 



'(p-q + k- k t ) ^(q + h - k)i(k) ^(1 - /-(p))(l - /-(?)) a/1 + W) + 

+ HMfc) -k 1 + k)A cxt ( P -q + k- h) y/(i - /_(p))(i - f-(q))y/i + nk) - 



u<y{q), 



i{k) PS F 21 (p -k 1+ k)A CK \p -q + k- fc 1 ) A //-(p)(l-/-(9))v7(*) 



A cxt ( P -q + k-ki) ?sUq + h- ky(k)V^- f-(p))f-(q)Vm 

a 2 



Sf 



BS,abs 



(2n)i 



d 4 fcig(fc!) / dk Q 5{k 2 )Q{k a )u s {p) 



+ (-i)A cxt (P -q-k-kt) ^(q + h + k)ft-k)y/(l - /_(p))(l - /-(«)) V7(*) + 

+ (-M-k) psUp -h- k)A cxt ( P - q -k- kV(i-/_(p))(i -/_(«)) v7(*) - 
- #-fc) ^(p - fcx - fc) A oxt (p - g - fc - Ai) A //-(p)(i-/-(5))v / i + /(fc) - 



A ext b - g - fc - h) es? 2 (q + k 1 + ky(-k)y/(i - JlyfiTMy/i + f(k) 



u a {q). 



(3.38) 



Thus we see that the absorption is the same as the emission with k replaced by — k (except 
in9(fco)), l + /by/. 

Inserting the ex pressi ons for the propagators from App. A. 3 in ( 3.38| ) - the definition of 
< ... > is given in ( fo.86 ) - using p 1 ^ = —fy + 2pe, etc. and the Dirac equation we get the 
following 2 times 6 terms: 



S 



BS.cm 
fi 



(27T)i 



d«hg{ki) / dk S(k z )e(k )u s (p)A ext (p -q + k- kx)u a {q) 



(2pe - ih + #)(-i)V(l - /-(p))(l - /-(g)) V 1 + /(*) 



(2tt) 2 ((p - fei + k) 2 - m 2 + iO) 



2pe-, 



W(i~f-(pW'f-(q))V^+W)- 



2ir 

6{( P -h + fc) 2 - to 2 ) < 0/ ± (p - ki + k) > + 
f {2pe - 



h + #)(-l)v//-b)(l - f_(q))^f{k)( — ) • 

• ^(p-fci + fcf-TO 2 ) <e v //±(p-fc 1 + fc)> + 



(3.39) 



(2ge + ft / - ^ (-*) V(i-/-(p))(i-/-(g))Vi + /(*;) 



(27r) 2 ((q + /ci -fc) 2 -m 2 +i0) 



2qe + fat- 
2n 



V(l-/_(p))(l_/_(g))VTT7(fe) 



■<J((g + *i - fc) 2 - m 2 ) < 6/±(g + fci - fc) > + 



+ (2qe + hi- - f-(p))f-(q)Vm(^-) ■ 

Z7T 

■ 6{(q + ki - k) 2 - TO 2 ) <e^f ± {q + k 1 -k)> 
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The matrix-element for the absorption S^ ' 3, s is the same, but with k «-> — k (except in 
0(*o)) and / «-» 1 + /. 

3.2 Calculation of some Third Order Graphs 

In this section we calculate the third order graphs that play or could play a role in the 
discussion of the IR-behaviour of the theory in cross sections to fourth order. Thus we 
will calculate T3 for the Vertex and the third order self energy graph. These have to be 
considered in investigating the IR-problem since they will turn out to be infinite in the 
adiabatic limit g — > 1. We also discuss the third order vacuum polarization and the graphs 
with two bremsstrahlung-photons - they are finite in the limit g — > 1, therefore we will not 
consider them further. To caution the reader we want to emphasize that in the following 
calculations we do not consider (T = 0)- and (T>0)-parts separately. Therefore, all formulae 
have to be taken literally for the (T > 0)-parts only, the (T = 0)-parts (and (T = 0)-parts 
of subgraphs with u> > 0, e.g. in the third order SE) are understood to have been treated 
correctly, i.e. by splitting the m non -trivially. We do not write down this explicitly (cf. also 



the remark at the end of Sec. 3.1.2 



3.2.1 The Vertex Graph 

D 3 is calculated in appendix p.l| , its singular order as well. We will only consider the 
terms with a non-tilde field A v in x 3 (since this will be replaced by a classical external 
potential in our calculations for the IR-behaviour) and or $ in ii, ^ or ^> in x-x. The 
terms with permuted arguments give the same result in momentum space (as can be shown 
by calculating just as we do for the terms we are concerned with here). Thus this gives 
an additio nal f actor 3! = 6 from the various permutations. But from the expansion of the 



^-matrix, (2.1), there is another global factor, ^ = i, which cancels the one above. 
Here we give the result for and som e hin ts for the calculation: 



We take the first part of D 3 (see App. DA) and change its form in a suitable wa; 



D% = p S- ^5+ ?D F - S iet S + f3 D+ - (>S- ^ av "D- - 
- S+ f>S- t3 D F + ^ av S~ D~ + S+ <% ret ?D+ = 

/3^»av 0gtTet 0J^-\~ _j_ /3gret fljjav ^^ av Pj^) rc ^ _i_ ^ 

■ /3gret (3ga,v &J^~\~ /3grct fijj av _j_ /Sg'av /^^ret 



We take the retarded part of it (by trivial splitting since to < —1, see App. D.l ) which gives: 
Rj = - S~ p S av D mt + ^S* 1 ' ' p S av P D + - S ret S~ D av . (3.41) 

Then we have 

T 3 = R 3 -R' 3 =- S- S av D rct + S ret S av f3 D+ - ^S 10 * f3 S~ f3 D av 

+ l3 S-(x 1 - x 3 ) S+ D F - ^ ret f3 S+ f3 D+ - S~ ^ av D~ = (3.42) 

= _ PgF PgF f3 D F 

To obtain this result we have added S~ 5 av D F - S~ S av D F and S~ S F D F - S~ S F ' D F , skilfully 
splitted S ret S F D F in S Tet S av D+ etc. and used the fact that supp(S Tet S Tet D av ) = {0}. 

37 The arguments of the products in the following formula are (xi — x-j), (x^ — X2), {xi — £2), always in 
this order. 
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Considering all factors we thus get: 



T 3 V = - e 3 : V( Xl )r fSf^xi - x 3 )>f "Sfi(a* ~ a*) ^(zi - x 2 ) 1 ^(x 2 ) : 4,(13) 

tffojy "S&fa - x 3 )Y ^(x, - a*) P D F 12 {x x - x 2 ) 1 ^{x 2 ) : 4(z 3 ) 

*( Xl )r "Sfifri - x 3 )Y PSY 2 (X 3 - x 2 ) D F 2 { Xl - x 2 ) 1 ^[x 2 ) : A v (x 3 ) 

*( Xl )^ ^( Xl - x 3 ) 7 " ?S F 2 (x 3 ~ x 2 ) PD* 2 ( Xl - x 2 ) 7 ^(x 2 ) : A v (x 3 ). 

(3.43) 

3.2.2 Third Order Self Energy 

There are two graphs that contribute to the third order SE: 'A', with non-tilde fields A v 
and ^ in x 3 and \& or 4 r in x\, and a self energy insertion on the incoming leg; and 'B', with 
non-tilde fields A v and ^ in x 3 and ^ or \f in x\ , and a self energy insertion on the outgoing 
leg. We first consider graph 'B'. For some more remarks on the choice of the external legs 
have a look at the begin of Sec. [00] . The calculation of £>3.ord.SE and w 3.ord.SE = _ QO ig 
done in appendix p.2| . Here we just give the result for y3-°rd.SE_ ^ e s -p]ft £)| ordSE trivially 
and build T = R-R / : 



T. 



3.ord.SE,B 



= e 3 (z 3 ) ( - 7" ^S F x {x 3 - x 2 )Y % F (x 2 - Xl )^ ^( Xl - x 2 ) 

- 7" % F 2 (x 3 - x 2 ) 7 ^ ^f 2 (x 2 - xx)^ D F 2 ( Xl - x 2 )^( Xl ) : A v (x 3 ) + 

+e 3 :*(aj 3 ) ( - 7" "Sfifo - x 2 )>f p S F 2 {x 2 - ^ F (x! - x 2 ) 

+ 7" ^ F 2 (X3 - x 2 )Y S F 2 (x 2 - a; 1 ) 7 „ ^ F (^i " ^2)^1) : A, (a*). 



(3.44) 



Graph A is calculated in the same way: 



T. 



3.ord.SE,A 



= e 3 ( - 7 M % F i(^i ~ 2:2)7^ % F i(^2 - 0:3)7" ^Fifo - asi) 

- 7 M "Sf 2 (a* - x 2 )7m ^(x, - z 3 )7* "Df 2 (a: a - a*))*^ : 4(z 3 ) + 
+e 3 ( + 7 M ^2 F i(^i - a; 2 )7„ % F i(*2 - x 3 ) 7 " D F 2 (x 2 - x x ) 

- 7 M Sl 2 {x x - x 2 )7m ^2 F i(^2 - x 3 ) 7 l/ ^22(^2 - xi)j *(x 3 ): A v (x 3 ). 



(3.45) 



3.2.3 Third order Vacuum Polarization 



£)3.ord.VP 

for T 3 OTd VP : 



and w 3 ord VP are calculated as in Sec. 3.2.2. This leads to the following expression 



^{x 2 )^^(x 2 ) D F 1/12/etc Xx 2 - x 3 )n 11/21 ^(x 3 - x^A^ixx). 



(3.46) 



We now insert the correct expressions for the vacuum polarization n and for D F . Disregard- 
ing some factors that are inessential with respect to convergence properties, we thus arrive 
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at 



d 4 k 2 d%g(k 2 )g(h) ( + v /(l-/_(p))(l-/_( (Z )) ^Df x (p - q - fe)n£(p - q - k 2 - k 3 ) 

+ v/(i-/_( P ))(i -/_(?)) ^i 2 (p - <z - fc 2 )n^(p - 9 - ** - fc 8 ) 

+ V.f-(p)f-(<l) P D^(p - q - feajfiff (p - 9 - fca - fca) 

+ Vf-(p)f-(q) dUp - q ~ k 2 )ii%(p -q-k 2 - k 3 )j . 



n is given by T 2 VP . This leads tcQ 



(3.47) 



d 4 k 2 d 4 k 3 g(k 2 )g(k 3 ) / d 4 fc (~ fc 2 + - fc + const.) 

, ' 

'7-matrix structure 7 



(2tt) 2 (jp-q- k 2 ) 2 + iO 2tt 



+ 7r s ((p-v- k 2) 2 )f(p-<i- k 2) 



1 



fe 2 _ m 2 + j 2tt 



— 5{k 2 -m 2 ) 



(Q - k) 2 - m 2 + i0 ' 2tt 
f-(Q ~ k)0(Q o - fc ) +f+(Q- k)d(-Q + k ) 



f_(-k)9(-k Q ) + f + (-k)e(k ) 
+ 7^5((Q ~ k) 2 - m 2 ) • 

y(l-/_(p))(l -/_(<?)) + 



+ ^*((P - 9 ~ k 2 ) 2 )y/f(p -q-k 2 )(l + f(p-q- k 2 )) ■ 



2tt 



<5(/c 2 - m 2 ) 



9(-ko)y/f-(-k)(l - /-(-*0) " f(Ao) - /+(-&)) 



2tt 



^((Q-A:) 2 -m 2 ) 



- fe )V/-(Q-A)(i-/-(Q-fc)) 



-Qo + *o)V/+(Q-*)(i-/+(0-*)) 



V(i-/_(p))(i -/_(?)) + 



+ ^((p - 9 - k 2 ) 2 )^f(p -q-k 2 )(l + f(p-q- k 2 )) ■ 



1 



k 2 - m 2 + iO 2ir 



+ —5(k 2 - m 2 ) 



f_(-k)e(-k ) + f+(-k)6(k ) 



1 



{Q-k) 2 ~m 2 + iO 2tt 
/_(Q - £O0(Q O - fco) + /+(Q - k)6(-Q + *o) 



38 Here we set 



p — q — k2 — ks. 
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2n 



8{k 2 - m 2 ) 



-^^ + i-, ((p _ g _ fc2)2)/( ,_ ? _ fc2) ). 

0(-koWf-(-k)(l - f-(-k)) - 9(k )Vf+(-k)(l - /+(-*)) 



2vr 



5{{Q - kf - m 2 ) 



e(Q - ko)y/f-(Q-k)(l-f-(Q-k)) 



0(-Qo + k a )^/f + (Q-k)(l-f + (Q-k)) 



V/_(p)/_(g)l . (3.48) 



Performing the multiplications and investigating the adiabatic limit we see that it exists and 
is finite (pay attention to the fact that the whole expression as a distribution need not make 
sense evaluated at specified points p, q. It has to be smeared out in these variables, therefore 
p = q causes no problems) . Thus the third order VP is not involved in the discussion of the 
IR-problem in the cross sections to fourth order (we will see later (see after ( |3.88| )) that a 
third order graph contributes with twice its real part to the fourth order cross section). 

3.2.4 Graphs with Two Bremsstrahlung-Photons 

This contribution has th e following form (here we choose special external legs again; cf. the 
discussion in Sec. 3.2.1 ): 



(3.49) 



d A x 1 d A x 2 d A x^{x 1 )fi.^ t {x l ) l3 Sf 1 (x 1 -x 2 )Mx2) ■ 
■ f) Sf 1 (x2-x 3 )Mx3)y{x 3 )g(x 1 )g(x 2 )g(x 3 ). 

We form matrix elements with emitted/absorbed photons (described by (a^ (k , (3) + VL p\ and 



(a M (fc, /3) + Qp\, see Sec. 3.1.4[ ) of momentum k and k', which are equal to 

jd 4 k 1 d A k 2 d 4 k 3 y/l-f-(p)u s {p)A CXt (p-q-k-k'-k 1 -k 2 - k 3 ) 
i(k) 



■ ^(q + k + k' + k 2 + k 3 ) y/l + f(k) 



0S{ 1 (q+k'+k 3 )y/l+f(k> 



2\k' 



■Ua(q)g(ki)g(k 2 )g(k 3 ). 

(3.50) 



Here we can have y/fjjk) instead of a/1 + f(k), or 12-, 21- and 22-propagators instead of 
11-propagators. In addition, not both propagators have to depend on q + k+ one or both 
of them can as well depend on p + k+ .... But all these modifications do not change anything 
in the following considerations on the existence of the adiabatic limit in these contributions. 
In the further calculations we use 



+ m)fu ~ u^f^u, 



(3.51) 



and the corresponding for (j> + m) instead of one or both of the factors (^ + m) . Here we set 
k = k! which is valid in the contributions from two bremsstrahlungs photons to the cross 



section to fourth order, cf. Sec. 3.4.3. 

If we calculate now the cross section and integrate over smalj^J k and k', it is clear that 



we have to investigate integrals of the following form (here p 

39 Since possible divergencies origin from small values of k and k' . 



q and 'P' denotes 



3G 



the 'principal value' of the integral): 



/ 



/•Wo ru ^ pi pi 

d 4 kid 4 k 2 g(k\)g(k 2 ) / / drdr'r 2 r I I dxdx 
Jo Jo J-lJ-l 

_Lp \ (3.52) 

rr' (q r' + q r + q{k\ + k 2 ) - \q\r'x' - \q\rx)(q r - \q\rx + q(ki + k 2 )) 



This is finite and so are the corresponding integrals with one or two (5-distributions from 
S F . Thus the graphs with two Bremsstrahlung-photons need not be considered in order to 
investigate the IR-behaviour in cross sections to fourth order. 



3.3 Thermal Corrections to 
3.3.1 Calculation 



In section 3.2.1 we calculated the vertex graph, i.e. T$ , eq. (3.43). One part of it will 
be important for the discussion of the IR-problem in the cross sections to fourth order (see 
Sec. 3.4). The other parts will give some thermal corrections to the magnetic moment of 
the electron. In this section we will be concerned with the calculation of these corrections. 
Before we explain how the vertex contributes to the electron magnetic moment, we first 
investigate its structure more closely to be able to identify the before- mentioned parts. To 
do that we consider the vertex in momentum space, i.e. its Fourier tran sform . Since we will 
be interested in electron scattering to determine fx e - (see below, after ( 3.70|) ) we study the 
matrix element 

S fi = (b s (p,(3) + ^\T^K(q,P) + n p ) (3.53) 

and write it down explicitly with Fourier transformed quantities. This leads to 

1 



d xid x 2 d x 3 g(xi)g(x 2 )g{x 3 )(-e 



(27T) 



Vl-/_(p)u(p)e^ 



• r "SFl(*l - X 3 )-f ^fl(*3 - X2h» DUx X - X 2 ) ■ 

1 -y/l-f^q)u(q)e-^A^ + 



the other three terms of T 3 = 



d xi 



2.3 



(2tt) 6 



rf 4 fci,2 : 33(fci,2,3)e 



— ikiXi—ik2X2—ik^X3 



(2tt) 3 



■r 



(2tt) 



va - f-(p))(i - f-(q)Mpy pxi - iqx2 ■ 

J 

d 4 rA c ^{r)e- irX3 + 



(27T) 



d 4 Q l3 Sf 1 {Q)e- lQ(xi - X3 ' > Y 



(27T) 
1 

W) 2 

+ the other three terms of T, 



7p ' 



d 4 k ^(^e-^-^uiq) 



V 



(3.54) 



Having performed the x- integrations we get three ^-distributions, each with a factor (27r) 4 . 
The ^-integration is trivial and gives a factor (2tt) 2 . In the other three terms we proceed 
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in the same way, which finally gives 



(2tt) 3 



d i k 1 d A k 2 g(k 1 )g(k 2 )u s (p) 



+ v/a-Z-CpWl -/-(<?)) / d 4 k^ ^i(p - fci - fc) 7 " ^(5 + - k) D F n (k) 



- V/_(p)(l-/-(g)) y dV ^i(P - *i - *h" XiOz + *a - *) ^DS>(A) 
+ (-i)V(l-/-(p))/-(g) / d 4 fc7 M ^i F i(P - h - k)Y P S F 12 {q + k 2 -k) ?D F 12 (k) 



Vf-(p)f-(q) I ^i(p - *i - ^OV ^2(9 + fc 2 - k) f3 D F 22 (k) 

■ l»K Xt (P -Q-h- k 2 )u a {q). 



(3.55) 



In this expression, the factors 7^ @S F (p — k\ — k)^f" l3 S F (q + k 2 — k)^^ lead to the follow- 
ing '7-matrix structure' of the vertex cont ribution (see the explicit representation of the 
propagators in momentum space, App. A. 3.): 



7" ^{p -h- k)Y ^(q + k 2 - k)^ ~ 7^ - jti - jt + m)-f {4 + + m) 7/i . 

(3.56) 

Sin ce we are inter ested in the adiabatic limit of these quantities (cf. the discussion in Sec. 
3.4 and after eq. ( 3.58 )) and this can be perfo rmed trivially, we can omit ki, i = 1, 2, in the 
numerator (cf. as well the end of Sec. |2.1.l[) and get (using the definition j) := p\j x and 
7^7" +7^ = 2g^ v , 7 A 7^7a = -27^ etc.) 



Y L {j> - J6 + m)Y(4 -jt + m) 7Al = 
{-2ii'j) + 4m<f + 4mp' y - 2m V) + ffi'ft + Hl v $ - %mk v ) - 2$-?% 



(3.57) 



Because of u(p) on the left and u{q) on the right of the whole expression ( |3.55[) this gives 
(using —2<fp( v j) — 2jrffc — 4p v i — 4q v ft + Apq-f" and the Dirac equation) 



u(p) Apqj v + A{ft(p v + q") 



mk v - (pk + qk)j v ) - 2ft-y v ft u(q). 



(3.58) 



The first of these terms, independent of k, is the one we will discuss in section 3.4 where 
we investigate the IR-limit. It is the only term of ( |3.58 ) that leads to IR-divergencies. Here 
we will be concerned with the finite part of the vertex, for which the adiabatic limit can be 
performed trivially without problems. The goal of the following calculations is to determine 
the thermal corrections to the magnetic moment of the electron, fi c - . We will not do that 
in the most general setting, w e will restrict ourselves to 'small' temperatures: k B T <C m c - . 
Fro m th e four terms of T$ , ( 3.43 ), then only the first survives, and only the third term 
of ( 3.58 ), proportional to k x k p , contributes. In this approximation, the first term of ( |3.58] ) 
does not contain any finite parts and gets cancelled completely by the contribution of the 
bremsstrahlung (see Sec. |3.4| ), and the second term (proportional to fc A ) vanishes identically 
as can be seen by symmetry arguments^. 
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dder the transformation k — > —k in (3.55) after inserting the explicit form of the propagators from 
i and considering only the terms ~ fc A in the numerator. 
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Thus we have to calculate the following part of the matrix element^ ( [3.55 ) 



(27T) 



This can be written as 



d 4 k(k^k v ) T —^ — — • — — - — ^(k 2 )- 1 



(3.59) 



{-2pk + iQ) {-2qk + iQ) y > e mA -\ 



(27T) 

with 



' u s (p)(l^ U lxK(l)ATn X (p,q,u), (3.60) 



T^(p,q,u) :=-2 /d 4 fc<5(fc 2 )-^ -. 

v 7 7 V Vl fc "l-1 (-2pfc + i0)(-2qfc + i0) 



(3.61) 



This integral has to be calculated. Because of the Lorentz covariance of this expression and 
the symmetry with respect to interchanging p <-> q it must be of the form 

T^(p, q, u) = AgT+BWlT + q»q v ) + C{p»q v + q»p v ) + 

+ D(// + // + + + G(//). ' ' 

Now we take the trace in this expression or multiply it with p^p v , p^qu, P^Uv or u^Uy, 
respectively, and thus get the following system of five equations: 

= 4^ + 2m 2 B + 2pqC + 2{pu + qu)D + G (3.63) 

h := / ^MC ^ih.!) ' ^k + iO) - + (™ 4 + (W) 2 ) 5 + 2m ^ C + 

+ (2m 2 (pu) + 2{pq){pu)^D + {pufG 

(3.64) 

h:=~\ ! d 4 kS(k 2 ) I = (pq)A + 2m 2 {pq)B + (m 4 + ( M ) 2 )c + 

_> ./ e/w-i v y (3 65) 

m 2 (qu) + m 2 {pu) + (pq)(qu) + (pq)(pu))D + (pu)(qu)G 



14 :- 



rf 4 H(fc 2 ) (e/3|M _^ • ( _ 2gfc+ . Q) 

=(pu)A + (m 2 (pu) + (pq)(qu)\B + (m 2 (<?u) + {pq){pu)j C + ( 3 - 66 ) 
+ (m 2 + (pu) 2 + (ptjf) + {pu){qu)jD + {pu)G 

1 /Un^ 1 ( fcM ) 2 



d kS(k 

2 J v ; e /3|fe«l_i (pjb-i0)(gfe-i0) ( 3 67 ) 



((pu) 2 + (<jM) 2 ) J B + 2(pM)(g M )C+ (2(pu) + 2(gu))L> + G. 



41 So this is the temperature dependent part of the first term of (|3.55|) in the approximation kgT -C m 
pe 

get a factor (27r) 2 from each fci-intcgration 



In addition we performed the adiabatic limit: we set fc; = in the arguments of Ty" and A® xt , g = 1 and 
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We will solve this linear system later. 

Now we explain why it is the vertex that determines the magnetic moment. First we 
define A Q := Ty l '7 M 7 Q 7i,. By the definition of the 7-matrices we have then 

A" = 2T£V„7p ~ ?T7m7,7 q - (3.68) 

With the Gordon decomposition^], 

u s (p)(jf + q' l )u a {q) = 2mu s (p)y i u a (q) - i{p v - q^a^u^q), (3.69) 



we get from ( p.62|) , after some manipulations and taking the u on the left and the u on the 
right end into account, 

A a (p, q, u) =7" (-2A+ 2m 2 B + Am 2 C - 2(pq)C - 2(pu)D - 2(qu)D - G) 

^ ' (3. TO) 

- 2im{B + C)(p p - q^)<r af * + 2u a (2mD + jiG) + 2{p a + q a )fiD. 



To extract information on /i c - from ( 3.70 ) we proceed as described for the case T = in 
@], p. 236. We consider electron scattering on an external static C-number potential at low 
energies. We calculate in the rest frame of the heat bath: u = (1,0) and use the standard 
representation for spinors with the same normalization convention as in |22| : 



E + m ( X- 



^(£> = \/^H^-y I' (3J1) 



where E := yjp 1 + m 2 . This is most convenient for the non-relativistic case |p| <C m e - we 
are interested in. Using this approximation and 

n _ A \ j _ ( a 2 \ klrn fa m 



where cr, are the Pauli-matrices, we find to leading order in — 



'3 

,.n 



Usipj^u^q) = [xt,Xv), 
u^fj^u^q) = ^- (xt, {p + q + iSy. {p - q)) 3 Xa) , (3.73) 
u s (p)<J^{Pn - q^)ua(q)A u = -(xt^Xa)((p- q)xA). 

These expressions will now be used to calculate the matrix element describing the before- 
mentioned electron scattering w e are inte reste d in. In the low temperature approximation 
this matrix element is given by (3.53) or (3.59), respectively, and just this term contributes 
as explained above. To compare with the value of /i - without any corrections, we consider 
the first order matrix element, which is given by 

(b s (p,P) + ^\TiK(q,(3)+n ) (3.74) 



in addition, where the T\ here consists of the non-tilde part of ( 2.68 ) only, since A° xt has to 
be a classical potential, i.e. of non-tilde type. In momentum space this is given by 

l e^-u s {p) 1 a u a {q)AT(p-q)- (3.75) 

42 Use a* 1 " := ^MSt"!. 
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Altogether, the first and third order matrix elements ( |3.74| ) and (3.59) are equal to (use 



Sfi ■= ie.—u s {p) 



=2 , 

—3 7 Q + 2m 2 B + 4m 2 C* - 2(pq)C - 2{pu)D - 2{qu)D - G 



(27T) 



2im(B + C) (p„ - % ) a<*» - — -g 2u a (2mD + jiG)- -— 3 2 (p a + q a )f,D 



(2tt) 3 



(3.76) 



Now we substitute ( |3.73| ) in and getg 



S 



fi 



ie 
2^ 



1 «i,ext „2 

' 2m (27r) J 



2 ' 2 (2mD + G)iS xt - t|4?^(p + q°)A c x 



2m(27r) 3 
2e 2 



IU(p + g)A 



(27T) 



(2tt) 3 



2m 2m(27r) 



(2tt) 3 
i2me 2 

W) 



D{p + q)A 
3 (B + C) 



where 



(3.77) 



:= -2A + 2m 2 B + Am 2 C - 2(pq)C - 2{pu)D - 2{qu)D - G. 



(3.78) 



The second summand in Q3.77 ) is proportional to the magnetic fieldf^ B(p — q) — i{p— q) x 
A{p — q) 7 therefore we can identify the magnetic moment^ 



e 
2m 



2 C 2 



1 - 



(2tt)= 



■W 



4m e 
"(2^ 



-(B + C) 



Using ( 3.63| ) and (3/75 ) this gives 



2m 



(2tt) : 



A2A) 



(3.79) 



(3.80) 



Thus we only have to know A, the other parameters are irrelevant for the magnetic mome nt 
i n th is approximation. Nevertheless we cannot avoi d so lving the whole system ( 3.65 ) - 
( 3.67 ). The solution will be calculated in appendix D.3. We only need the solution in 
the case of p — q, therefore we investigate the solution of the system for q := p — rj with 
q 2 = (p — rj) 2 = m 2 and rj 2 — if <C m 2 and then consider the limit rj — + 0. Pay attention 

43 As mentioned above: we calculate in the system u = (1, 0). 
44 The following formula is the Fourier transform of B = rotA. 

45 These are the first order term and the thermal correction; the QED corrections at T = are not included 
here. - To distinguish factors that do not belong to ji c - we use the fact that is the magnetic moment of 
the electron determined in the framework of relativistic QM To set p = q to identify the magnetic moment 
stems from a less intuitive discussion of expressions like (3.77) as it is done in [B4J, p. 347. There it is shown 
that the quasi-static limit p = q allows the discussion and identification of ji e - quite similar to the one in 
the framework of relativistic QM. 
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to the facht that due to the above definitions we have 770 = po — Po y 1 — (^p? — |?)- We 
emphasize that setting p = q in the linear system and then solving it is in general not 
allowed, but here it leads to the same result. We will discuss that and what could go wrong 



this way at the end of Sec. 3.3.2 



Now we present the result for the magnetic moment of the electron, including the QED- 
correction for T = and the corrections due to T>0, each up to third order in e (see App. 



D4 eq. ( PJ6Q ): 



M "--2^i 1 + 8^~36^rf>- !:iXI1 



3.3.2 Discussion of the Result 

In this section we want to discuss the result ( |3.81| ) for /i c - a little bit and to compare it with 
the existing results in literature. We are aware of nine papers wherein /i c - is calculated^: 
§ and |uj - |7). The results of these papers disagree partially; we can find three different 
values, which we report here, in units of 'magnetons' jp-, where in some of the references 
temperature dependent values for e and m are used. But in the approximation k B T -C m e - , 
which we are considering here this temperature-dependence can be omitted (we refer to the 
cited literature for a more exact treating, here we are just interested in giving an idea of 
what the results in literature are like and how they differ from our) . Now we give the values 
from the literature^]: 



18/3 2 ™ 2 " 
e 2 

36/3 2 m 2 ' 
19e 2 



reported by |10| - |13|, |§ and |17[, 

reported by § and @, (3.82) 

reported by pa . 



So all the results are of the same magnitude with respective factors 1, 2, ~ |. We agree with 
the second result and in absolute numbers it reads approximatively (for T = 293 K, h and 
c not set equal to 1 anymore) 

M e - = -0.62- 1(T 17 • ^- . (3.83) 
2m 



For comparison we give the value of the (T = 0)-corrections up to the same order: ^ ■ = 
1.16- 1CT 3 -^. 

Among all these calculations there is just one that is formulated covariantly, pi, - as is 
ours. We agree with them that all the differing results quoted arise because it is quite tricky 
to identify the magnetic moment of the electron in a non-covariant setting. 

To caution the reader we emphasize that there could be a difference between our result 
and the one from pj] that could be traced back to their setting p — q too early (cf. the 



remark after (D.14): Because the determinant of the linear system is zero for p — q it is not 
allowed to set p — q at the beginning. The procedures of taking the limit q ~ p — r/, 77 — ^ 
and solving the equation system do not commute!). If we set p = q^rj = QdX the very 

46 There are other papers that give values for /^ c - , e.g. Jli| ] or |pc|] , but for the calculation they refer to 
one of this nine papers. One paper, Ejj , we could not find. 

47 [|l I] has the wrong sign due to some error; it is corrected in [[L 3| . And | |lq | argues that a QFT calculation 
should reproduce a correction of — 3 g^ m ^ , which was obtained by some other method, an effective Hamilto- 

nian method, incorporating some intuitive ideas on the nature of the electron and not using covariant QED. 

2 

He then identifies the value of — 5j — ^ as part of it and gives some arguments why the calculations that 



lead to the first result in (3.85) just can account for T of the desired value 
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beginning of our calculation, i.e. in (D.l), we get five equations. Two of them are identical. 
We throw away one of them, arguing with ( [3.62 ) that B and C must not be distinguished 
for p = q and therefore just four variables, A, B, D, G, have to be determined. Solving this 
new linear system we end up with 



.4 



9/3 2 m 2 ' 

which fortunately is the same result as we have - but this will not be true in general 



(3.84) 



3.4 The Infrared Problem in thermal QED in Cross Sections to 
Fourth Order 

In this section we investigate the adiabatic limit in the cross section for scattering of an elec- 
tron on an external classical C-number potential to fourth order. To do that we have to con- 
sider the contributions from the vertex, the third order self energy and the bremsstrahlung 
with one BS-photon. There are no other second order contributions and the other third 
order graphs that could contribute are finite (see Sec. 3.2.3 and 3.2.4). We calculate the 
cross sections for these processes, which will contain IR-divergent parts in the limit g — > 1. 
Then we investigate if the various divergent parts perhaps cancel. This is not the case for 
arbitrary temperatures, whereas it is in the approximation k B T <C m e - . We caution the 
reader that we often write down the expressions for the (T = 0)-parts as if they had been 
calculated using ordinary Feynman rules. This is not the case, they have to be treated 
correctly according to the causal theory as described in |2.1.2 . But we have chosen to write 
them down in this simple manner to keep track of them without complicating the notation. 
This is possible since the T-independent parts can be treated separately (cf. end of Sec. 
3.1.1) and their cancellation was prooved in p2|. The same remark applies to some terms 



of the temperature dependent part of the third order SE including sub-graphs of singular 
order u > 0. 

Thus we have to calculate the cross sections for these various processes, i.e. 



da_ 



4. ord. 



C-SE 

°/i,3.ord. 



Sfi, 3. ord. 1 2 



4. ord. 



5 em 
, /i,2.ord. 



2 + l^abs 



7 /i,2.ord. 



— |*S/t,l.ord 

nem / nem \* _i_ cabs / r»abs \ * 

— D /i,2.ordA D /i,2.ordJ "+" J /i:,2.ord. l D /i,2.ord J 

+ S/M-ord.^^OTd .)* + £/i^.Sd.(£/i,l-ord.) 



To shorten the notation in the calculations we define the following quantities: 
<0^jm)>:= (ofaWf-ipW - TM) - 0(-paWf+(p)0- - /+(P))) ■. 



(3.85) 



< 0f±{p) > := 0(po)/-(p) +0(-pt»)/+(p) 



(3.86) 



Since we are interested in the adiabatic limit we set fcj = wherever this is already possible 



(according to the remark after (2.32)) 



3.4.1 Cross Sections of the Involved Processes 

Vertex: We consider the matrix element ( |3.53 ) 

S fl = (b s (p,P) + n \TV\b a (q,l3) + ^). 



(3.87) 
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We proceed as in Sec. 3.3, but now we are interested in the first, i.e. fc-independent part in 
(3.58). Since the vertex is of third order it contributes to the cross section ( 3.85| ) with 



S'/i,l.ord.(5'/i)* + <S/i(S7i ! l.ord.)* ■ 



(3.88) 



From ( 3.58 ) we see that the interesting part of the vertex is proportional to u(p)^ u u(q) , i.e. 



proportional to Sf 



i,l.ord. 



Using Sj. =: hS 



'/i.l.ord. 



and by inserting the explicit expressions 
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for the propagators into (3.55), we thus have 



da_ 

^; 4 . ord . 



|S/i,i.ord.| 2 2Re(/i) 



= — u 3 (p)fi {p-q)u a {q)u a {q)A ip ~ q)u s (p) — y/l - f-ip)^/ 1 ~ /-(«) 



: ext,* 



•Re 



d 4 /c 1 d 4 fc 2 5(fc 1 )3(/c 2 )8( M ) / cffc 



{2-Kf 



1 



1 



(2tt) 2 (p-ki- k) 2 - m 2 + iO 



+ S({p - k, - k)' - m 2 ) <Qf ± (p-k 1 ~k)> 

Z7T 



(2tt) 2 (<? + k 2 - k) 2 - m 2 + iO 



+ S((q + k 2 -k) 2 -m 2 ) <ef ± (q + k 2 -k)> 

Z7T 



1 -1 



(2tt) 2 k 2 + iO 2tt 



+ -s(k 2 )f(k) -/-(<?)) 



i[ (-1)— <5((p-fc 1 -fc) 2 -m 2 ) <e v //±(p-fc 1 -fe)> 
1 1 



(2tt) 2 (q + k 2 - k) 2 - m 2 + iO 



+ —S((q + k 2 - k) 2 - m 2 ) <e/±(g + fc 2 -fc)> 

Z7T 



<5(fc 2 )V/(fc)(l + /(fc)) V/-(p)(1 -/-(g)) + 



(2?r) 2 (p-h- k) 2 -m 2 + iO 



+ Sdp-k.-ky-m 2 ) <Qf ± (p-k 1 ~k)> 

Z7T 

(-l)^S((q + k 2 -k) 2 -m 2 ) <ey/f ± (q + k 2 -k)>y 



■ (-6(k 2 Wf(k)(l+f(k))j^(l - f-(p))f-(q) - 
((-l)^((p-ki-k) 2 -m 2 ) <OVf±(p-ki-k)>y 

■ ((-i)^((9 + fc 2 -fc) 2 -™ 2 ) <&Vf±(q + k2-k)>y 

-1 1 



(2?r) 2 -fc 2 + iO 2tt 



5(fc 2 )/(fc) K//-(p)/-(g) 



(3.89) 
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Performing the multiplications this gives the following 14 terms 



16 : 

1 : 

2 : 

4 : 

3 : 

5 : 

6 : 



'da 

on 



V 



+ 
+ 
+ 

+ 



4.ord. 

— i 



da\ lold - 8(pg)e 2 
dn ) (2tt) 2 



Re 



d kid k 2 g(k 1 )g(k 2 ) 



(i-/_(p))(i -/_(?)) 



+ 



+ 



(2tt) 6 ((p - h - fc) 2 -m 2 + iO)({q + k 2 - k) 2 - m 2 + iO)(k 2 + »0) 

_1 (l-/-(p))(l-/-( g )) 2 

(2tt) 5 ((p-fci-fc) 2 -TO 2 +iO)((g + fc 2 -fc) 2 -m 2 + zO) ^ m ' 

1 (l-/_(p))(l-/_(g)) 
(2tt) 5 ((p - h - k) 2 -m 2 + iO)(k 2 + iO) 

• 5((g + fc 2 - fc) 2 - ™ 2 ) < 9/±(g + fc 2 - fc) > + 
-j (l-/_(p))(l-/_(g)) 

(27r) 4 ((p-fci-fc) 2 -m 2 +i0) 

• <f((g + fc 2 - fc) 2 - m 2 )5(k 2 ) < Gf±{q + k 2 -k)> f(k) + 
1 (l-/_(p))(l-/_(g)) 

(2tt) 5 ((g + k 2 - fc) 2 - to 2 + i0)(fc 2 + iO) 

■ 5((p - h - fc) 2 - to 2 ) <e.f±(p-fci-fc)> + 
-j (l-/_(p))(l-/_(g)) 
(27r) 4 ((g + fc 2 -fc) 2 -m 2 +i0) 

• <y((p - fc! - fc) 2 - TO 2 )<5(fc 2 ) < 0/±(p - fc! - fc) > /(fc) + 

(2tt) 4 (fc 2 + iO) 

• 5((p - fci - fc) 2 - to 2 ) < 0/±(p - fci - fc) > • 

• 5{(q + k 2 - fc) 2 - m 2 ) < 9/±(g + fc 2 - fc) > + 
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1 (1 - /_(p))(l - f-(q))5({q + k 2 - fc) 2 - m 2 ) <0/±(g + k 2 - k)> 



(27T) 

• <J((p - fci - fc) 2 - m 2 )5{k 2 ) < 0/ ± (p - fci - fc) > /(fc) + 



+ 7 ^ (1 - / - ( ; )V /r 2 (p)(1 2 - / - ( f )) ^((P - *i - fc) 2 - ™ 2 ) • 

(2tt) 4 ((g + fc 2 - k) 2 - m 2 + iO) u ' ' 

■ <eVf±(p-ki-k)> <5(fc 2 )v//(fc)(i + .f(fc)) + 

1 (1 - - /-(p))*((p - ^ - fc) 2 - m 2 ) • 



(27T) 



• < Q\/f±(p-ki-k)> 6((q + k 2 - kf - m 2 ) • 

• <6/ ± (g + fc 2 - k)>S(k 2 )^f(k)(l + f(k)) + 

■ <e^/f ± (q + k 2 -k)> 5(k 2 Wf(k)(l + f(k)) + 



1 



r(i - .f-ip)Wa-f-(q))f-(q)s((p -h- k) 2 - m 2 ) 



(2tt) 3 

• <e.f ± (p - h - k)> S((q + k 2 - k) 2 - m 2 ) 



■ <QVf±(q + k2-k)> 6(k 2 )^f(k)(i + f(k)) + 



-i y/-(p)(i-/-(p))/-( g )(iT(g)) 

(2tt) 4 A: 2 - iO 

• <J((p - fci - fc) 2 - to 2 ) < eV/±(p-fci-fc)> • 

• <5((g + fc 2 - fc) 2 - m 2 ) < e v / /±(<Z + fc 2 -fc)> + 

+ 7^3 V/-(P)(1 - /-(p))/-(g)(l - /-(g)) • 



5((p-fc!-fc) 2 -m 2 ) <e v //±(p-fci-fc)>- 



<*((« + fc 2 - fc) 2 - m 2 ) <e v //±( (Z + fc 2 -fc)> <5(fc 2 )/(fc) 



(3.90) 
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Self Energy: We consider the matrix element 

s fl = {b s (p 7 P) + ^\Tl E K(q^) + ^), 



(3.91) 



with T 3 SE from Q3.44| ) and ( |3.45| ). Proceeding as we did with the vertex to obtain ( |3.55[ ) we 
get in the momentum space representation 



S 



3.ord.SE,A _ e 

f* ~ (2tt) 3 



d i kid A k 2 g{k 1 )g{k 2 )u s (p) 



+ V(1-/-(P))(1 -/-(?))■ 

d * k [ - y (p + fc - fc^ ^(fc)] ^(p - fca - fc 2 ) + 



+ V(i-/-(p))(i -/-(?))• 

- y ^f 2 (p + fc - fc 1 ) 7 „ ^ F (fc)] ^(p - A* - fca) + 



Vfc[ + 7 M ^(p + fc - fci)7„ ^izW] ^(p - fci - fc 2 ) + 



+ V/-(P)(1 -/-(«))■ 

- 7 M ^5 2 F 2 (p + fc - fc 1 ) 7 „ ^ F (fc)] ?SUp -h- k 2 ) 

■7"^(9)^(p-g-*i-*a) 



(3.92) 



and 



-r3.ord.SE,B _ e 

"(27T)3 



d fcirf k 2 g{ki)g(k 2 )u s {p)~f v 



+ V /(1-/-(P))(1 -/-(<?))■ 

d 4 fc[ - ^(g + fci + fc 2 ) 7 ^ ^S F (g + fc + fcxhp ^(fc)] + 

! V L--/_(p))(l -/_(<?))• 

d 4 fc[ - ^ F (g + fc a + fc 2 ) 7 ^ <£f 2 (g + k + h)^ D\ 2 {k)] + 



+ -/-(p))/-(g)- 

- ^S F (g + fc x + fc 2 ) 7 ^ ^(g + fc + fc 1 ) 7/1 ^ F 2 (fc)] 



+ H)v / (i-/-(p))/-(?)- 

+ ^ F (g + fc! + fc 2 ) 7 M ^(g + k + k,)^ ^(k)] 
•u CT (g)i° xt (p~g-fci -ha). 



(3.93) 



As in the vertex contribution (see after ( 3.88 )) Sff is proportional to S/jji.ord. and thus the 
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cross section is given by 

SE 



( da\ 2 / c 3.ord.SE.A . c 3.ord.SE,B\ 

V^A.ord. = ' /M - 0rd -' { " ^ > 



= ^s s (p)A (p-gK(gK(«)4 ' (p- g)« a (p) — V(l - /_(p))(1 - /_(?)) • 
•Re 8m 2 J tPkKpkigikJgikt) Jd 4 k- 

+ Vd -/-(P))d -/-(«))(-!) • + i 5(fc2)/(fc) ) ' 

— + -^—5({p + k — k\) 2 — to 2 ) < Qf±(p + k — k\) > 



(2tt) 2 (p + fc-fci) 2 -m 2 +i0 2tt 
1 1 



(2tt) 2 (p - fci - fc 2 ) 2 - m 2 + iO 



+ 



+ ^5(( P -k 1 -k 2 ) 2 -TO 2 )<e/±(p-fc 1 -fc 2 )> ) + 



+ V(i - /-(p))(i - /-(<z))(-i) (^(fc 2 )VMi + • 

• (^((P + *-*i) 2 -m 2 ) <©v / /±(P + fc-fc 1 )>) • 

• (^5{{p-ki-k 2 ) 2 -m 2 ) <9V/i(p-fci-y>) + 

+ V/-(P)(1 - /-(?)) (^((P + fc - fc i) 2 - <®V.f±(P + k-k 1 )>y 

+ ^((p -ki- h) 2 - to 2 ) < e/±(p - fcx - fc 2 ) > ) + 

+ iV/-(P)d - (^to + sH ' 

! + _! ( 5( (:P + a : -a :i ) 2 -to 2 ) <e/±( P + fc-fc 1 )> 



(2tt) 2 (p + fc-fci) 2 -to 2 -iO 2tt 
(^H(p-ki-k 2 ) 2 -m 2 ) <e v // ± (p-fc 1 -fc 2 )>) + 
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+ Vd -/-(P))d -/-( 9 ))(-D ( (^p^o + SF^*)) ' 
1 1 
(2tt) 2 (g + fci + fc 2 ) 2 - to 2 + iO + 

+ + fc 2 ) 2 - to 2 ) <6/±(g + fci + fc 2 )> 

Z7T 

1 1 



(2tt) 2 (q + k + ki) 2 - m 2 + iO 



+ ^5({q + k + fci) 2 - to 2 ) <6/±(<z + fc + fci)> ) + 
+ v/(l - /_(p))(l - f-(qM-l) [^5(k 2 Uf{k){\ + JW)) 
■ (^S((q + h+ k 2 ) 2 - m 2 ) < Q^f±{q + k 1 +k 2 )> 

^-5({q + k + h) 2 - m 2 ) <6V/±(? + Hfci)>j + (3.94) 

+ (-i)V/-(«)(i-/-(p))(-i) (^(fc 2 )V/(fc)(i + /W) 

1 1 



(2tt) 2 (g + fci + fc 2 ) 2 - to 2 + iO 



+ irS({q + h + k 2 f - to 2 ) < e/±(g + + fc 2 ) > 

Z7T 



• (^((g + fc + fc!) 2 -™ 2 ) <eV/±(? + fc + fci)> 

+ (<?)(! - /_(p)) (^((<Z + fci + fe) 2 - to 2 ) <0 A // ± ( 9 + A:i+A:2)>) • 



" 1 ^ + ^(* 2 )/W 



(2tt) 2 -fc 2 + iO 2tt " v 7 V(2tt) 2 {q + k + kif - m 2 - iO 
+ ^S((q + k + k 1 ) 2 -m 2 ) <e.f ± (q + k + k 1 )> 
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After having performed the multiplications, this gives the following 30 terms: 



o : 



a : 



b : 



c : 



d : 



e : 



f : 



g : 



da_ 



SE 

4.ord. 

— i 



da_\ 

dn) 



Lord. 



8m e 



2^2 



■Re 



/ 



d 4 k 1 d 4 k 2 g(k 1 )g(k 2 ) 



(2n) 2 

(l-/_(p))(l -/_(«)) 



+ 



(2tt) 6 ((p + k- fci) 2 -m 2 + i0)(k 2 + iO)((p -h- fc 2 ) 2 -m 2 + iO) 

1 (l-/_(p))(l-/_(g)) 
(2?r) 5 ((p + k - fci) 2 -m 2 + i0)(k 2 + iO) 

■ 5((p -h- k 2 ) 2 - m 2 ) < 0/ ± (p -k 1 -k 2 )> + 
-1 (l-/_(p))(l -/_(«)) 



+ 



(2tt) 5 ({p + k — fci) 2 - to 2 + z0)((p - fci - k 2 ) 2 - m 2 + iO) 



^ 2 )/(fc) 



+ 



i (l-/_(p))(l -/_(<?)) 



5(fc 2 )/(fc) 



5((p - fci - fc 2 ) 2 - to 2 ) < Of±(p -k 1 -k 2 )> + 



(2tt) 4 ((p + fc-fci) 2 -TO 2 + iO) 
,-k^-k^-m 2 ) 
1 (l-/_(p))(l-/_(g)) 
(2tt) 5 (fc 2 + i0)((p - fci - fc 2 ) 2 - to 2 + iO) 

•<5((;p + fc-fci) 2 -m 2 ) <e/±(p + fe-fei)> + 
i (l-/-(p))(l -/-(«)) 



+ 



+ 



(2tt) 4 (fc 2 + »0) 

• <e/ ± (p + fc-fei)><e/±(p-fei -fc 2 )> + 
-* (i-/_(p))(i -/_(«)) 



(2tt) 4 ((p - fci - fc 2 ) 2 - 



S((p + fc - fci) 2 - m 2 )(5((p - fci - fc 2 ) 2 - m 2 ) 
P 

S(k 2 )f(k) 



iO) 



(2tt)3 



• ^((p + fc-fc!) 2 -?^ 2 ) <e/ ± (p + fc- fc 1 )> + 
(i-/_( P ))(i-/_( g Mfc 2 )/(fc). 

• <f((p + fc - fcx) 2 - m 2 )(5((p - fci - fc 2 ) 2 - m 2 ) 

• < 9/±(p + A; -&i)><e/±(p- fci -**)> + 
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+ (^3 (1 - /-(P))(l - f-(q))S(k 2 )Vn^ +JW)- 



6({p + k- fci) 2 - m 2 )S({p - fci - fc 2 ) 2 - m 2 ^ 



<eVf±(p + k-k 1 )><e^/f ± ( P -k 1 -k 2 )> 



i (i - /-( g )V/-(p)(i - f-(p» S (k 2 W .f(kKi + / W) • 



(2tt) 4 (p - fci - fc 2 ) 2 - m 2 + iO 

■ sdp + k-hf -m 2 ) <e v / /±b + fc-fci)> + 

+ ^)s(l - f-(l)Wf-(p)(l - /-(p))<5(fc 2 ) v //(fc)(l + .f(fc)) • 

• <J((p - fci - fc 2 ) 2 - m 2 )5((p + k- fci) 2 - m 2 ) • 

• < e/±(p - fci - fe) >< e V/± (p + fc - to) > + 



, -l (l-/-(g))v / /-(p)(i-/-(p)) 

(2tt) 5 ((p + fc - fci) 2 - to 2 - iO)(-fc 2 + iO) 

• 5((p - h - k 2 f - m 2 ) <e v //±(p-fc 1 -fc 2 )> + 

(27r) 4 (p + k — ki) z — m z — iQ 

5((p - h - k 2 f - m 2 ) <e^f±( P ~k 1 -k 2 )> + 



i (i-f-(q)Wf-(p)(i-f-(p)) 

(2tt) 4 -k 2 + iO 

■ S((p + k- fci) 2 - m 2 )5((p - fci - fc 2 ) 2 - m 2 ) • 

• <e.f±(p + k-k 1 )><e^f ± ( P -k 1 -k 2 )> + 

i 

(2^) 3 



+ 77^3 a - fAq)Wf-(p)(i- f-(p)Mk 2 )f(k) ■ 



■ 8{{p - fci - fc 2 ) 2 - m 2 )5((p + fc - fci) 2 - m 2 ) 

• <e/±(p + fc-fc 1 )><e v // ± (p-fci -fc 2 )> - 
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, _=» (i-/-(p))(i-/_( g )) 

(2tt) 6 ((g + fc + /ci) 2 - m 2 + i0)(fc 2 + iO)((q + fa + k 2 ) 2 - m 2 + iO) 
1 (l-/_(p))(l-/_(g)) 



<5(fc 2 )/(fc) 



(2tt) 5 ((g + fc + fa) 2 - m 2 + iO)(k 2 + iO) 

■ 5({q + ki + k 2 ) 2 - m 2 ) <Of±(q + fa +k 2 )> + 
,_1 (l-/_(p))(l-/_(g)) 

(2tt) 5 ((g + fc + fci) 2 - m 2 + iO)((q + fa + k 2 ) 2 - m 2 + iO) 

-i (l-/_(p))(l-/_(g)) 
(2tt) 4 ((g + fc + /ci) 2 -m 2 + iO) 

• 5{k 2 )f(k)5((q + fa + k 2 f - m 2 ) < 9/±(g + fa + k 2 ) > + 
1 (l-/_(p))(l-/_(g)) 
(2tt)5 (fc 2 + iO)((g + fci + fc 2 ) 2 - m 2 + iO) 

•5((g + A: + fc 1 ) 2 -m 2 ) <6/±(g + fc + fci)> + 

+ ^ 1 ' f ^% Uq)) ^ + ^ + *i) a - + fa + fa) 2 - m 2 ) 

■ <®f±{q + k + fa)><Of±(q + fa + k 2 )> + 
-i (l-/_(p))(l-/_(g)) 



Ak 2 )f(k) ■ 



(2tt) 4 {(q + fa + k 2 ) 2 -m 2 + iO) 

■ 5({q + k + fa) 2 - m 2 ) <Qf ± (q + k + fa)> + 

+ j^(l-f-(p))(l-f-{q))Kk 2 )f{k)- 

■ S((q + k + fa) 2 - m 2 )6{{q + fa + k 2 ) 2 - m 2 ) • 
• <&f±(q + k + fa)><Qf±{q + fa +k 2 )> + 
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h' : 



k' 



-1 



+ 7ra(! - /-(P))(l - f-(q))S(k 2 )Vf(k)(l + f(k)) 



(27T) 



<5((g + fc + fci) 2 - m 2 )<5((g + k x + fc 2 ) 2 



<eVf±(Q + k + k 1 )><e^f ± (q + k 1 + k 2 )> + 

-8{k 2 )^f(k){l + f{k)) 



-i (i-/_( P )) > //-(g)(i-/-(g)) 



(2tt) 4 ((g + fci + fc 2 ) 2 - m 2 + 2 0) 

• <5((g + fc + fci) 2 - m 2 ) < e y /f±(q + k + k 1 )> S(k 2 ) + 

1 



-(1 - /_(p)) - /-(g)^(fc 2 ) V/(fc)(l + /(fc)) 

• <5((g + fe x + fc 2 ) 2 - m 2 )<5((g + fc + fci) 2 - m 2 ) ■ 



■ < e/±( g + fci + k 2 ) >< e v/±(g + fc + fci) > + 

-1 (l-/-(p)V/-(g)(l-/-(g)) 
(2tt) 5 ((g + k + fci) 2 - m 2 - i0)(-fc 2 + iO) 



•<5((g + fci+fc 2 ) 2 -m 2 ) <Q^f±(q + k 1 + fc 2 ) > + 



<5(fc 2 )/(fc) 



2 (l-/_(p)V/_( g )(l-/_( g )) 
(2tt) 4 ((g + fc + fci) 2 -m 2 -i0) 

• 5(( 9 + fci + fc 2 ) 2 -m 2 ) <eV/±(g + fci + fc2)> + 
i (l-/-(p))v//-(g)(l-/-(g)) 



(2tt) 4 -fc 2 + iO 

• <J((g + fc + fci) 2 - m 2 )<5((g + fci + fc 2 ) 2 - m 2 ) • 

■ <e/±(g + fc + fci)><ev//±(g + fci + fc2)> + 
+ ^3 (i - /-(p))V/-(g)(i-/-(g))<s(fe 2 )/W • 

■ 5((g + fci + fc 2 ) 2 - m 2 )8((q + fc + fci) 2 - m 2 ) • 

• <e/±(g + fc + fci)><eV/±(g + fci + fc2)> | 



(3.95) 



Bremsstrahlung: The bremsstrahlung is of second order, therefore (see (3.85)) it con- 
tributes with 



BS,em/ abs / qBS,ciii/ abs\ * 
fi 



Q DO, LULL/ fclUS / qD 



(3.96) 



to the cross section. The matrix elements are known from ( 3.39 ). This gives for the cross 
sections 



48 The factor 1 = J d 3 kS(k — fc') is inserted in order to have k' independently from k in the second factor. 
We omit the terms proportional to fa in the numerators since we are interested in the adiabatic limit and 



these terms could at most give finite contributions (as will be 



terms proportional to pe or qe. These are investigated in Sec. 3.4.2). On similar grounds we can omit the 
terms proportional to ft. 



clea 



from the divergence properties of the 
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/ o \ BS,cm , 4 » „ „ 

Mj} k \ = (2^)9 J d4k ^9(kiW(k2) J dk a 5(k 2 )e(k ) J dk' Q S(k' 2 )0(k' a ) 



BS.em 

(k) 
4.ord. 

d 3 k'5(k - k>)\ u s {p)^ xt {p-q + k)u^q)u a (q)^*(p-q + k')u s {p) 



(peV(l-/-(p))(l-/4g))y/r + 7(fc) 

((p-fci+fc) 2 -TO 2 +iO) 

+ i2n(pe) V(l - /-(?)) (1 - /-(g)) Vl + /(*) • 

• S((p - fci + fc) 2 - m 2 ) < 0/±(p - fe! + k) > + 

+ i27r(pe) V/_(p)(l - /-(g))v7(fc)<5((p - fci + fc) 2 - ™ 2 ) < e A //±(p-fci + fc)> + 



( ge y(i-/_( P ))(i-/_( g )) v /r+7w 

((g + fci - fc) 2 -m 2 + iO) 

+ i27r(ge)V(l - /-(P))(l - /-(g)) V / l+7(fc) • 

• + fcx - fc) 2 - m 2 ) <e/±(g + fc!-fc)> + 

+ 12^(56)^(1 - /_(p))/_(g) v7W<^((9 + *i - fc)' - ™ 2 ) <©V / /±('Z + fci-fc)> 



(p e )V(l-/_(p))(l-/_( g )) v /l + /(A/) 
((p - fc 2 + fc') 2 -m 2 - iO) 

- i27r(pe) V(l - /-(P))(l - /-(?)) a/1 + /(fc') ' 

•<5((p-fc 2 + fc') 2 -m 2 ) <e/±(p-fc 2 + fc')> - 

- i2*{pey v //-(rt(i-/-(g))v / 7(fc 7 ) • 

• 5((p - fc 2 + fc') 2 - m 2 ) < 6 v //±(p-fc 2 + fc')> + 

(g e )V(i-/-(p))(i-/-(g))v / i + /(fc') 



+ 



((g + fc 2 - fc') 2 - to 2 - iO) 



iM^T v / (i-/-(rt)(i-/-(g))v / i + /(fc') • 

• 5((g + fc 2 - fc') 2 - to 2 ) < 6/±(g + fc 2 - fc') > 



- i27r(ge)V(l-/-(p))/-(g)V7(feO 



• 5((g + fc 2 - fc') 2 - to 2 ) < ev//±(g + fc 2 -fc')> 



(3.97) 



We get the terms for the absorption from the ones for the emission by changing fc <-» — fc 
(not in O(fco)) and / <-> 1 + /. Performing the fco-, fco- and fc'-integrations we get a factor 
-Ar, and we see that fc' = fc. In order to compare this with the contributions from the vertex 

2\k\ 

and the self energy we integrate over small photon energies |fc| < loq <C rn e -, where too can 
be interpreted as some finite energy resolution of the detector. Since we are interested in 
the adiabatic limit we can replace < Qf±(p — fci + fc) > by < Qf±(p) >, etc. Furthermore, 
we set p 2 = m 2 = q 2 . Performing the various multiplications gives the following 2 times 36 
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terms^J: 

-m 2 (l-/_(p))(l-/_(g))(l + /(fc)) 



A 



(2p(k - ki) + (k — k x ) 2 + i0)(2p(k - ha) + (k - k 2 ) 2 - iO) 



C : + 2 P (k-k 2 + (k-k 2 ) 2 )-^ 1 - /-(P))^ - + **» 



■5(2p(k-k 1 ) + (k-k 1 Y) <8/±(p)> + 
-2-Kim 2 



2p(fc - fca) + (k - fc 2 ) 2 - 2 



v7-(p)(1-/-(p))(1 - f-(q)Wf(k)(l + f(k)) 



■ 5(2p(k - ki) + (k- fcO 2 ) < 6 y/fjfi > + 
v . -pg(l-/_(p))(l-/_( g ))(l + /(fc)) 

' ■ (2?(*i - fc) + (fci - fc) 2 + i0)(2p(fc - k 2 ) + (k- k 2 ) 2 - *0) 

Vn : + ^ ,r 2 T'^ 2 .-» (! - /-(?))(! - /-(?))(! + 



2p{k- 


k 2 ) + (k - k 2 ) 2 - iO 




■ 5(2q(k 1 -k) + ( 




—2nipq 


2p{k- 


hi) + (k- k 2 ) 2 - zO 




• S(2q(h -k) + ( 




2irim 2 



!-i : + szti its , n ,^ 2 , n V/-(g)(i-/-(g))(i - f-(p))y/Hk)(l + f(k)). 



• 6(2p(k - k 2 ) + (k - k 2 ) 2 ) <0/±(p)> + 
D : + 47r 2 (-m 2 )(l - /_(p))(l - /_(?)) (1 + f(k))S(2p(k - h) + (k - k x ) 2 ) 

■ d(2 P (k - fc 2 ) + (fc - fc 2 ) 2 ) <e/±(p)><e/±(p) > + 

0: + 47r 2 (-m 2 )V/-(p)(l - /-(p))(l - + /(*0) ' 

■ $(2p(fc - fci) + (k - k 1 f)6(2p(k - k 2 ) + {k- k 2 f) ■ 



<©V/ ± (p)><e/ ± (p)> + 



49 Here we summed over four polarizations of the photons, including the unphysical scalar and longitudinal 
ones (cf. S p. 245 and eq. (3.12.42)), and we used £n=l e^T^ei™' - e^*^ = -g^u. 
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• 6(2 P (k - k 2 ) + (k - fc 2 ) 2 ) <e/±(p)> + 

VIII : + 4^ 2 (- P<Z )(l - /_(p))(l - + f(k))S(2q(k 1 -k) + (h - kf) ■ 

■ s(2 P (k - k 2 ) + (k - k 2 f) <e.f±(q)><e.f ± (p)> + 

1 II : + 47r 2 (-pg)v//_( g )(l - /_(g))(l - /_(p)) A //(fe)(l + f(k)) ■ 

■ 5(2q(k 1 -k) + (fci - k) 2 )S(2p(k - k 2 ) + (k- k 2 ) 2 ) ■ 

■<evCT)><e/±(p)> + 
7: + 2p(fc _ fcl) + (fc _ fcl)2 + , V/-(P)d - » - + 7(*)) • 

• $(2p(fc - fc 2 ) + (fc - fc 2 ) 2 ) < 6v</^(rt > + 

S: + 47r 2 (-mV/-(p)(l - /-(P))(l - f-b)Wf(k)(l + JW) ■ 

■ 5(2p(k - h) + (k- k l f)8{2p{k - k 2 ) + (k- k 2 f) ■ 

■ <e/±(p)><e A //±(p)> + 

e : + 4^ 2 (-m 2 ).f_(p)(l - f_(q))f(k)5(2p(k - h) + (k - k,) 2 ) ■ 

■ s(2 P (k - fc 2 ) + (k - k 2 f) <e A //±(p)><e A //±(p)> + 

i-ni: + 2g(fcl _ fc) 2 + 7- fc)2 + io ^Z-Wd - /-(p))d - /-(«W7W(i + 7W) • 

• 5(2 P (fe - fc 2 ) + (fc - fc 2 ) 2 ) < e > + 

l-IV : + 47T 2 (- M ) v //_(ri(l-/_(p))(l - .Uq)Wf(k)(l + f(k)) ■ 

■ 6(2q(k 1 - k) + (fci - k) 2 ) ■ 5(2p(k - k 2 ) + (k - fc 2 ) 2 ) • 

<e.f±(q)><Q^Mp)> + 

l-V : + 4TT 2 (-pq)V.f-(p)(l - /-(p))/-(g)(l - f-(q))f(k) ■ 

■ 6(2q(k 1 -k) + (fci - k) 2 )S(2p(k - k 2 ) + (k- k 2 ) 2 ) ■ 

■<®^Mq)><®^Mp)> + 
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-pg(l-/_(p))(l-/_(g))(l + /(fc)) 
(2p(jfe - fci + (fc - fci) 2 ) + i0)(2q(k 2 - fc + (fc 2 - k) 2 ) - iO) 

(1 -/_(p))(l -/-(?))(! + 



2g(fc 2 - fc) + (fc 2 - fc) 2 - M 

■8{2p{k-k 1 ) + {k-k 1 f) <Of±(p)> + 
— 2iripq 



2q(k 2 - k) + (fc 2 - k) 2 - iO 



W/-(p)(l-/-(p))(l f-(q))y/f(k)(l + f(k)) 



■ d(2 P (k - k t ) + (k- fci) 2 ) < eVMp) > + 

-m 2 (l-,f_(p))(W_( g ))(l + ,/-(fc)) 
(2q(k 1 - fc) + (fci - k) 2 + i0){2q{k 2 -k) + (k 2 - k) 2 - iO) 

■ S(2q(h - *) + (*i - *) 2 ) <e/±(«) > + 

■ «2 9 (t, - k) + (ti - 1) 2 ) <e-SJM) > + 

27rip5 



■(l-/_(p))(l-/_(g))(l + /(*)) 



2p(k - fci) + (fc - fci) 2 + iO 

• <5(2g(fc 2 - fc) + (fc 2 - fc) 2 ) <ef±(q)> + 

+ 47T 2 (- M )(1 - /_(p))(l - + f(k))5(2p(k - fci) + (fc - fci) 2 ) 

• <5(2g(fc 2 - fc) + (fc 2 - fc) 2 ) <e/±(p)><e/±(g)> + 
+ 47r 2 (- Pg ) v //_(p)(l-/_(p))(l - Uq))y/f(k)(l + f(k))- 

■ S(2p(k - fci) + (fc - fci) 2 )(5(2g(fc 2 - fc) + (fc 2 - fc) 2 ) • 

■<ev^W><e/±(«!)>i 
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• <5(2g(fc 2 -k) + (fc 2 -k) 2 )< e/± (g) > + 
H : + 47r 2 (-m 2 )(l - /_(p))(l - + f(k))6(2q(k 1 - k) + (h - kf) 

■ S(2q(k 2 -k) + (k 2 - kf) <ef ± (q)><Of ± (q)> + 



+ 47r 2 (-m 2 V/-(<7)(W-(9))(l - /-(p))V/(fc)(l + /(*))• 

• 5(2g(fci -k) + (fci - fc) 2 )<5(2q(fc 2 - k) + (fc 2 - A:) 2 ) 



<©VM9)><e/±(?)> 



i-vii: + 2p(fc _ fcl) + ( r_ fcl)2 + i0 ^-w - - /-(p))V7wa + 7W) ■ 

• ,5(2 g (fc 2 - fc) + (fc 2 - fc) 2 ) <ev^(^)> + 

l-VHI: + 4^ 2 (- TO )v//-(<z)(l - /-(<Z))(1 - f-(p))Vf(k)(l + 7W) ' 

• 5(2p(k - h) + {k- k 1 f)d(2q(k 2 -k) + (k 2 - kf) • 



■ <ef±{p)><QVMF)> + 

l-IX : + 47r 2 (-p<zV/-b)(l - /-(p))/-(g)(l - f-(q))f(k) ■ 

■ S(2p(k - h) + (k- k 1 f)d(2q(k 2 -k) + (ha - kf) 



■ <®Vf^><evMoj> + 

27TITO 2 



2q(fa -k) + (fci - k) 2 + i0 



V/_( 9 )(l-/_(g))(l - /-(P))V/(A)(1 +/(*)) 



<5(2g(fc 2 - fc) + (fc 2 - fc) 2 ) <6V/±(g)> + 



k: + 4^ 2 (-m 2 ) v //_( <Z )(l - - /_ (p)) ^/(fc) (1 + /(fc)) • 

• *(2g(Ai - fc) + (fci - kf)5{2q(k 2 - fc) + (fc 2 - fc) 2 ) 



■<e/±(g)><ev^(g)> + 

+ 4^ 2 (-m 2 )(l - f-(p))f-{q)f(k)8(2q(k 1 - fc) + (fcx - fc) 2 ) 



• <5(2 g (fc 2 - fc) + (fc 2 - fc) 2 ) <eVMi)><QVMi)> j + 

+dito(fc <->-&,/<-► 1 + /). 
3.4.2 Investigation of Possible Cancellations 



(3.98) 



Now we compare the cross sections we calculated in the foregoing section. We investigate if 
there is a net cancellation of all the divergent terms or not. The result will be that there is 
no complete cancellation, certain contributions from the self energy and the bremsstrahlung 
do not cancel. The vertex on the other hand causes no problems. Furthermore, we will see 
that there is a complete cancellation if we have small temperatures k B T <C m c - . 

Concretely we will consider the following things: We will integrate the BS-contribution 
from |fc| = to too -C 1, which accounts for soft photons that remain undetected (here ujq can 
be interpreted as some finite energy resolution of the detector). This also fixes the Lorentz 
system: it is the laboratory system, i.e. the rest system of the bath, i.e. u = (1,0). The 
fc- integration in the V- and SE-contributions will then be considered for values |fc| £ [0,wo] 
as well, since possible divergencies origin from fc = only. Thus we can treat fco = |fc| as 
small and omit it where this is possible. Since we are interested in the adiabatic limit, we 



will set e = <-> fci = where this is allowed (cf. the remark after (2.32)). Thus we can 
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for example replace < Qy/ f±{q + ki + k 2 ) > by < 9 y/ f± (q) >, < Qf±{p + k — k 2 ) > by 
< Qf±(p) >, etc. In addition we choose the test function g to be symmetric and real, which 
is no big restriction since we are interested in the limit g — > 1. Now we list how several 
terms compensate and which terms remain in the end. 

Vertex and pq-terms of BS: We often use the fact that the two BS-terms that are 
considered together with one V-term are the complex conjugate of each other and thus the 
sum is twice the real part. 

-: V 1 + BS 1 + BS V gives the terms proportional to pq of the (T = 0)-BS multiplied with 
(1 — — /-(<?)) and compensates V 16 (cf. [p2| , Sec. 3.11). This can be seen as 

follows: V 1 has an integrand proportional to (since we are interested in the real part 
we consider the principal value part only) 

-f 

(3.99) 



(2p(-h - k) + (-ki - k) 2 + i0)(2q(k 2 -k) + {k 2 - k) 2 + iO) ' 

The BS-contributions BS 1 + BS V plus the corresponding terms from absorption give 
after some transformation of variables an integrand proportional to 

- 0(fco) f (3.100) 



(2p(k 1 + k) + (fci + k) 2 + i0)(2q(k 2 -k) + (k 2 - k) 2 + iO) ' 

Now we expand the factors in the denominators (using a formula like 'ytpj = 1 — 
x' adapted to this case) and compare the resulting terms. We see that the leading 
divergencies ad up to the (T = 0)-part of the bremsstrahlung and thus compensate 
V 16 . The next to leading order term is of the form 

p /( fc X fc + fc i) 2 I ( 3ioi) 

{2p(k + k 1 )) 2 2q(k 2 -kY [ ' 

The transformation k — > — k, ki — > —ki then generates a change in the global sign 
and thus the integral is identical zero. The other terms are finite and vanish in the 
adiabatic limit. 

V 4 + BS n + BS vn : Proceeding as above we get for the vertex an integrand proportional 
to 

+6(2p(-k - h) + (-k - k 1 ) 2 )8{2q{k 2 -k) + (k 2 - k) 2 )f(k) (3.102) 

and for the bremsstrahlung 

-S(2p(k - fci) + {k- k 1 f)5{2q{k 2 -k) + (k 2 - k) 2 )(e{k Q ) + f{k)). (3.103) 

With the help of S(k 2 ) it can be seen that the terms \k 2 — k) 2, etc. in the other 
^-distributions can be neglected. Then it follows that the divergent terms (the terms 
proportional to f(k)) cancel between BS and V and the /-independent term from BS 
is finite. 

V 5 + BS m + BS VI : analogous to the previous calculation 

V 7 + BS IV + BS vm : dito (Because we consider the real part only). 

V^ + BS^ + BS 1 ^ 11 : dito. 
v ii +bs i-ii + bs i-viii. ditQ 

V 12 +BS 1 " III + BS 1 " X : dito. 
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V 13 + BS 1-IV + BS 1-VI. ditQ 

V 15 + BS 1_V + BS 1_IX : dito. 

V 14 : as above since we are interested in the real part only. 
V 6 : dito. 

V 2 : The relevant structure is 

1 



((2p(-fci -k) + (-fci - fc) 2 ) + iO)(k 2 + iO) 



5{2q{k 2 -k) + (k 2 - k) ). (3.104) 



We use the ^-distribution to perform the fco-integration and therefore decompose it 

according to %(x)) = T,k,y(*k)=° ifeff ' We see that this gives W\^ k ° ~ 2q °^ + 
S(k )]. Using this in the complete expression for V 2 leads to a leading order divergency 
that vanish es because of antisymmetry under the transformation of variables ki — > — fci 
(cf. ( |3.10l| )). The next to leading order term is finite. 

-: V 3 analogous to V 2 . 

Self energy and m- terms of BS: We use the same methods as for the vertex; in addition 
we use the following partial fraction decomposition: 



1 1 



(2p(fc - fci) + (fc - fci) 2 + iO) (2p(fc - fc 2 ) + (k- fc 2 ) 2 - iO) 
! A B 



(3.105) 



(2p(fc - fci) + (fc - fci) 2 + iO) (2p(fc - ha) + (fc - fc 2 ) 2 - »0) ' 
where we demand A = — B, which gives 

A = 2p(k 1 - fc 2 ) + (fci - fc 2 )(2fc - fci - fc 2 ) - iO" (3 ' 106) 

-: SE b + BS A is discussed as the terms V 1 + in addition using (ft.l05| ). The leading 
order gives the terms proportional to m 2 of the (T = 0)-BS multiplied with (1 — 
/-(p))(l - /-(g)) and compensates SE° (cf. [||, Sec. 3.11). The other terms are 
finite or can be shown to vanish because of asymmetry under certain transformations 
of variables (cf. V 1 + ...). 

-: SE b + BS E + SE° : compensation in an analogous way. 
-: SE C + SE f : the relevant term is 

-4S(2p(k - fci) + (fc - k 1 ) 2 )8{2p(-k 1 - fc 2 ) + (-fci - k 2 ) 2 )8(k 2 )f, for (3.107) 
BS B + BS C it is 

2<5(2p(fc - fci) + (fc - k 1 f)5(2p(k - fc 2 ) + (fc - fc 2 ) 2 )(6(fc ) + /). (3.108) 

Using 8(k 2 ) we see that the terms ~ (fc — fci) 2 etc. in the argument of the other 
^-distributions are negligible. Thus there is no cancellation between the contributions 
proportional to / because of the factors '4' and '2', respectively in front of these 
expressions. The term proportional to 9(fco) is finite. This factor '2' which inhibits 
a cancellation shows up here since we cannot use the afore-mentioned partial fraction 
decomposition which precisely led to the compensation in the case of SE b + .... 

-: SE d ' + SE C ' + BS F + BS G : is investigated analogously. 
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-: BS D + SE S : dito. 

-: BS H + SE g ': dito. 

-: SE 1 + SE m + BS" + BS 7 : dito. 

-: SE 1 ' +SE 1 ' +BS" +BS': dito. 

-: SE j + SE n + BS P + BS S : dito. 

-: SE j ' + SE n ' + BS K + BS e : dito. 

-: SE e ,SE f ,SE 1 ,SE m ': These terms have the same structure as the finite contribution 
to SE C + ... . Thus they are finite as well. 

-: SE a ,SE a ',SE k ,SE k ',SE d ,SE c ': Here we have the case of a subgraph with singular 
order uj > 0. The self energy part has to be splitted using the causal method. The 
result is given in p2| , eq. (3.7.41). If we have a close look at that for the momentum 
p — eki with e — > 0, we see that the terms that could diverge vanish because they 
are antisymmetric under the transformation of the integration variable ki — > — ki (cf. 
V 1 ...). 

-: SE h + BS e : These terms contain three ^-distributions that are treated as in V + ... 
. More important, however, are the factors that multiply them. For the self energy 
it is 2y/ f(l + /)(1 — /_(p))(l — f-(q)), whereas for the bremsstrahlung it is /(l — 
f-(q))f-(p) (the term proportional to 6(fco) is finite). Thus there is no cancellation 
because of the factor '2' and the different structure in f±. 

-: SE h ' +BS A : analogously. 

We want to add some remarks on the crucial factor 'two' in the calculations above: each 
vertex contribution comes with two corresponding BS-contributions, thus resulting in a net 
cancellation. For each self energy contribution, on the other hand, there is just one BS- 
contribution. For SE b + BS A + SE° the partial fraction decomposition (3.105) leads to an 



additional factor two in BS A , thus leading to a net cancellation. The same considerations 
apply to SE b + BS E + SE° . However, there is no partial fraction decomposition ( |3.105 ) 



involved in the (T > 0)-contributions since one or two of these factors are replaced by a 
^-distribution. Thus there is no additional factor two for these terms and in consequence no 
cancellation either. 



3.4.3 Investigation of Possible Cancellations with Mixed States 

In this section we want to analyse possible cancellations using more general asymptotic 
states. It will turn out that not even the use of general mixed states or coherent superposi- 
tions leads to a cancellation of the IR-divergencies. 

First, we consider incoherent mixtures of states with one electron, on electron and one 
photon, ordinary and 'tilde', one electron and two photons, ordinary and 'tilde'. Higher 
numbers of photons do not have to be considered since we are interested in the cross section 
to fourth order only. From more than two photons some are merely spectators and not 
involved in the processes. 

We will not use a density matrix formulation for the problem. Thus we will have to 
calculate the cross sections corresponding to the various parts of the initial and final state, 
take the weightened mean over the initial contributions and the sum over the final ones. We 
write for the initial state (we have suppressed to write down the vacuum on which the 
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emission operators act) 

^initial) :=+ / d i q\b + (q))(b+---\c (q) 



MS 



V / A / d 3 fci|& + (g)a+(fc 1 ))(---| Cl ( g ,fc 1 ). i + 

i JMS J<u 

V / A / d 3 fci|& + (g)5+(fci))(---|£i( g ,fci) l + 

V / A / d'fci / d 3 fc 2 |6+(( Z K+(fc 1 )a+(fc 2 ))(. ■ -102(9, fci,*2)« + 

- J MS J<cj -/<oj 

V / A/ ^ 3 fci / rf 3 fc 2 |& + (g)a+(fci)5+(fc 2 ))(.--|5 2 (g,fc 1 ,fc 2 ) 4 , + 

4 • JMS J<w J<w 

V / A / d'fci / d 3 fc 2 |6+( <Z )a+(fc 1 )a+(fc 2 ))(.--|5 2 (g,fc 1 ,fc 2 ) y , 

■ ■ JMS J<Wn J<Wn 



(3.109) 



The g-integration is over the mass shell (MS), the fci i2 -integration over |fci. 2 | < loq and the 
sum runs over four polarizations^"] i, j — 1, ...,4. The /3-dependence in the various thermal 
emission operators is suppressed in the notation. Since the Cj's etc. represent probability 
densities we have cq, ci, c 2 > and the normalization condition 



/ d 4 q\c (q)+y2 d 3 fcici(g,fci) 2 + V / d 3 fci£i(cj, fci) 4 

JMS 1 , J<u>o a J<ui 



/ d 3 h d 3 k 2 c 2 (q, fei,fc 2 )y ; +^2 / ^ 3fc i / d 3 k 2 c 2 {q,k x ,k 



j j- J <l*JQ J <UJ y J <UJ Q J <U(, 



id 



<u>a J <u> 



X <i 3 fci / d 3 k 2 c 2 (q, ki, k 2 )ij 



Vij + 



(3.110) 



Furthcron c 2 and c 2 are symmetric in the arguments fci and fc 2 and the indices 

The S'-matrix terms involved are the same as in the sections |3.4.1 and 3.4.2| . We group 
them according to the number of bremsstrahhmg photons (BSP) involved: zero BSP : the 
first order contribution and third order vertex and self energy [^J One BSP : the second 
order BS contributions. Two BSP: the three third order terms with two BS photons. 

Performing the sum and integration over the possible final states we get 36 contributions 
to the cross section, each involving a weighting factor Co, Ci, ... and integrals over the electron 
and photon momenta. They are listed in appendi x p]. 

Now we investigate the cross section given by (E.2). We first consider the case of sharp 
momenta for the incoming and outgoing electron, i.e. there is no integration over these 
momenta. Then^J the coefficient of |(6 + |0|fo + )| 2 equals 1 by the normalization condition 

( pioD - 

Using the ^-distributions, the fact that / ss 1 + / ~ /jjibj f° r the leading contributions to 
the IR-divergencies and the structure of applied to |a + ), etc. we see that the contributions 
proportional to the product of zero and two BSP factors are of the following form: 

50 We decided for shortness and completeness to indicate the four involved photon polarizations in this 
way directly in the emission operators. 

51 The third order vacuum polarisation could at most give next to leading order in divergencies (cf. Sec. 

ra. 

"In the following we write |0|, 1 1| and |2| in the matrix elements for the zero, one and two BSP-contribution 
to the S'-matrix. 
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d 3 k 1 (0)(b + a + {k 1 )\2\b + a + (k 1 )) • [ terms involving c u ...}. (3.111) 

Since the c's obey the normalization condition and are positive, no integral of the form 
J dkc(k) can show IR-divergencies. Furthermore, the products (0) (2) are finite as well (see 



Sec. 3.2.4), thus the whole contribution involving these products of zero and two BSP 
contributions is finite. 

Therefore, we only have to consider the zero BSP contributions, which are the same as 
in the calculation with pure states as already mentioned and the one BSP contributions we 
will examine now. 

Using the relations between (6+o + |l|6+), (6+a+|l|6+), (6+|l|6+o+) and (b+\l\b+a+) 
(k <-» — k and / 1 + / replacements are involved) and the fact that / 1 + / ~ ^tct for the 
leading contributions to the IR-divergencies we can write the leading one BPS contribution 
as follows: 



1BPS = 



|(6+a+(fc 1 )|l|6+)| 2 [c ( g )+^ j d 3 l lCl (q,h) m +J2 I dHxCi{q,h) 
+ cx(g, fci)j + y^ / d 3 l 2 c 2 (q, h,l 2 )in +E d3 ^C2(q,kx,h)in 

n n 

|(6+5+(fc 1 )|l|6+)| 2 [ Co (g)+^ J dH lCl {qM) m + Yl J d *hh{q,h) 

m m 

+ ci(q, fci)j + y^ d 3 l 2 C2{q,k 1 ,l 2 )in +Y d 3 l 2 c 2 (q,k l7 l 2 ) in 



(3.112) 

We look at the p 2 , q 2 - an d pg-c ontributions in the cross section factors of this expression 
separately as we did in Sec. 3.4.2 . Since the zero BPS contribution are the same as before, 
we see that each pg-contribution should finally give the same value as before, thus cancelling 
the vertex contributions. The p 2 , ^-contributions, however, should undergo a change which 
amounts to a factor 2 for the (T > 0)-parts, accounts especially for the SE-graphs SE h,h 
(BS A , BS C have to change in this way) and does not spoil t he ca ncellation of the (T = 0)-parts 
(i.e. BS A ,BS E should give the same as before); cf. Sec. 3.4.2. 

To investigate this further we write the one BSP cross section in the following form (cf. 
eq. (3.97)). Unimportant global factors are omitted and the '+' in the argument of F 
indicates the 0(+fco) involved): 



d 3 k\{b+(p)af(k)\l\b + (q))\ 2 = 

i J 

= ]T J d 4 h J d 4 fc 2 g(fc 1 )r(fc 2 ) J d 3 k- 

i 

■ [F( P , q, kx, k, B(k))i + F{q,p, -k lt -k, Q(k))i 
F(p, q, k 2 , k, Q(k))* + F(q,p, -k 2 , -k, 9(fe)X 



(3.113) 



The analogous term with d + instead of a + reads the same after having replaced k with — k 
(not in the 0-function involved). We have from (3.112) 
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1BPS = + ^] J d^hj d A k 2 g{k 1 )g*{k 2 ) J d 3 k< 

Z 

+ (f(p, q, fci, fc, +)iF(q,p, -k 2 , -k, +)* + F(q,p, -k 1 , -fc, +) l F(p, q, fc 2 , fc, +)*J 

• Q(q,uo) + Q'(q,vo,k)i + 

+ (P{p, q, fci, -fc, +)iF{q,p, -fc 2 , fc, +)? + -k u fc, +) i i f, 0, g, fc 2 , -fc, +)*) 

• w ) + Q'{q,ujQ,k)i | + 

+ Z! f d% J d i k 2 g(k 1 )g*(k 2 ) J d 3 ki 
+ (^(p, 9, fci, fc, 9, fc2, fc, +)* + F(qr,p, -fci, -fc, +)iF(g,p, -fc 2 , -fc, +)*J 

• Q(q,uo) + Q'(q,uo,k)i + 

+ [F(p, q, fci, -fc, +)iF(p, g, fc 2 , -fc, +)* + F(q,p, -fci, fc, +)iF(q,p, -fc 2 , fc, +)*J 

• Q(q,uj ) + Q'(q,ui ,k)i 



(3.114) 

The Q's abbreviate the terms involving the c's. Consider the leading contributions to the IR- 
divergencies in GG* , i.e. the terms '(fc — fc^) 2 ' can be omitted. Then each 'pg'-contribution 
equals the positive or negative of the corresponding l p 2 , q 2 '-contribution for p — > q. Using 
the continuity of G in p and q and G(p,p, ...) ^ we seef^] that it is im possib le to find Q's 
(and hence c's) that meet the afore-mentioned conditions (cf. after eq. ( 0.112 )) and lead to 
a cancellation of all IR-divergencies. 

The same consideration applies to the case of sharply peaked momenta centered around 
q and p of the in- and outgoing electrons, respectively, even if we allow small deviations from 
p 2 = m 2 = q 2 . The case of arbitrary coherent superposition of states can be investigated 
using similar arguments and does not lead to a cancellation as well. 



3.4.4 Discussion of the Results 

In this section we make several remarks on the results of the last one and discuss the results 
established in literature in relation to ours. In section ^ we will present some ideas how the 
problems that arise here in the IR-behaviour possibly could be cured. 

- First we want to emphasize that there is no cancellation of IR-divergencies in QED at 
finite temperature in the cross sections to fourth order if we do not restrict ourselves to 
small temperatures (see the second point). The vertex contributions to the divergency 
cancel with some part of the bremsstrahlung, whereas the self energy contributions 
do not. Especially the terms 'SE h ' and 'SE h ' cause problems since they are of a 
structure that does not arise in the bremsstrahlung (in 'SE h ' and 'SE h ' we have non- 
tilde external legs leading to a factor (1 — f_(p))(l — /_(?)) and an inner tilde vertex 
leading to 12-propagators. The BS terms either have a factor (1 — f_(p))(l — f-(q)) 
and no 12-propagator or a 12-propagator and the wrong factor involving /_). The 
other SE-terms that do not compensate differ by a factor two from corresponding 



53 For p = q, perform the transformation fci ,2 
see that we can write "(GG** + GG*)[2Q + Q' + 



—fci, 2 in the second and fourth summand. Then we 
for the sum of the first two and last two summands, 



respectively. 
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BS-terms. Thus we disagree with the existing literature that claims IR-finiteness 
for all temperatures ([Q, g), etc.) and have doubts if proofs of IR-finiteness to all 
orders ([Q) could rigorously be performed. We will discuss the shortcomings in the 
literature in point three. In addition to the external states considered usually we have 
investigated mixed states as well. We find the same result: there is no cancellation of 
IR-divergencies. 

Secondly we sec that the restriction to small temperatures k B T <C m c - finally leads 
to IR-finite values. In this case only V 1 , V 16 , SE b , SE°, SE b ', SE°', BS 1 , BS V , BS A 
and BS E contribute. In contrast to the literature we have this cancellation in a very 
straightforward way without the necessity to introduce temperature dependent spinors 
or a thermal Dirac equation as it is introduced in ([|S3|) and used by many other 
authors. Thus we see that the consideration of the Fermi-Dirac-distribution f±, which 
in itself is IR-finite causes more than just minor irrelevant changes with respect to the 
low temperature case (as claimed in Q], Q): the finiteness present in the latter gets 
completely lost. 

Pay attention to the fact that in this approximation the contributions with |fc| > luq 
from V 1 , BS 1 and BS V cancel exactly. That is why there is no contribution to the 
magnetic moment from this part of the vertex. 

In this third point we discuss several papers and the reasons for which they end up 
with finite results and why they are wrong or do not contradict our results. 



A first group of papers ([lfj, |L3|], |33| and j35| - |6j) presents calculations in the low 
temperature approximation k B T <C m c - only. As mentioned above we agree with the 
IR-finiteness found by them in this approximation; but as our calculations show, one 
cannot extract from this any results concerning arbitrary temperatures, and to use 
this finiteness as a hint for the existence of IR-cancellations to all orders at T > is 
too optimistic. 

A second group of papers (jij - 0) is concerned with arbitrary temperatures but 
neglects the thermal modification of the fermion propagators via f± on more or less 
clear reasons. Thus the fact that they find a complete cancellation of IR-divergencies 
does not contradict our results since we have this as well if we omit the terms with 
factors f±. 

The papers (Q, || and ||) leave out the doubling of fields in their cancellation and 
thus especially do not have to face the contributions SE h , SE h , which cause the major 
problems in our calculations. Especially in there are some terms missing: without 
reason, they claim to have no modification of the bremsstrahlung by the Fermi-Dirac- 
distribution and in their contribution to the vertex and the self energy they only have 
the terms V 2 , V 3 , SE a , SE d , SE a and SE C , which do not cause any problem anyway. 

In the paper |Q] there are some shortcomings in the considerations after their equation 
(32): In their general investigation of the n-th order they consider SjtS for each vertex 
and set g = 1 from the very beginning. But in the third order we are interested in, 
one has to investigate u^SftSj^u ■ b and the limit g — > 1 has to be performed at the 
very end. The net effect of the former treatment with respect to the latter is again a 
the omission of contributions as SE h , SE h . 

The paper |5| also neglects the doubling of fields. 

Besides that there are some papers that consider other problems than we did, |08) and 
f69| , to mention just two, which consider <p 3 and ip 4 theories, respectively, which have 
quite a different structure from QED. 

Furthermore it is clear from the discussion that there are a lot of finite contributions 
from these various terms considered above. As in the (T = 0)-case they possibly depend 
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on the test function g (see p^ |, Sec. 3.12). We do not consider this here any further 
since it is of no importance as long as the whole theory stays divergent in the adiabatic 
limit as it does now. 



4 Conclusion 

In this work we investigated QED at finite temperatures in perturbation theory. Using the 
causal approach to QFT we are not confronted with any UV-divergency and have a method 
at hand to rigorously discuss the IR-behaviour of the theory by investigating the adiabatic 
limit. 

We have calculated several second and third order graphs and derived corrections to 
m c - and [i c - due to the effects of the temperature. Our results agree with the ones in the 
literature. 

The main result, however, is the fact that the theory is not IR-convergent. Already in the 
cross sections to fourth order there is no complete cancellation of the various IR-divergent 
terms in the case of arbitrary temperatures. This result remains valid even if we set the 
chemical potential fi equal zero with the consequence that /_ = /+ or if we use more general 
mixed states or coherent superpositions as initial and final ones. For small temperatures the 
cancellation can be established in a straightforward way without the necessity to introduce 
thermal spinors. 

Now we give some ideas, which perhaps could be pursued to construct an IR-finite theory 
for QED at T>0. 

A first idea would be to use the normalization freedom of the causal approach due to 
the distribution-splitting with singular order u> > to get rid of the divergent terms in the 
adiabatic limit. But as a close inspection of these terms readily shows, this is not possible 
since these terms do not have the structure of a polynomial - not even in the adiabatic limit. 



Another idea is to change the decomposition of the field operators (2.57), (2.63) into 
positive and negative frequency parts or to define the thermal absorption and emission 
operators in a different way. But it seems quite impossible to do that in a reasonable way 
without violating the causality of the support of the Z)2-distribution. Since this property 
is crucial for the causal construction of the S'-matrix, the implementation of this approach 
seems impossible. 

The third - and perhaps most promising - idea is to change the involved scattering states 
in a suitable way. On grounds of the observation that the divergent vertex contributions 
cancel and only the self energy (and the corresponding bremsstrahlung terms) cause diffi- 
culties we could try to absorb these SE-contributions in the fermionic states by summation 



(proceeding as for the thermal mass corrections, see Sec. 3.1.3 ). This reasoning leads to 
a modified Dirac equation and temperature dependent spinors as it was done in |63| for 
the low temperature limit. Difficulties surely will arise through the two types of vertices 
that should be considered in a complete summation scheme and through the mathematical 
subtleties that origin in the structure of a thermal Dirac equation. 

Compare to the case of charged particles in the (T — 0)-theory. To have a unique finite 
adiabatic limit the unphysical degrees of freedom of the photon have to be considered in the 
bremsstrahlung (p^j, Sec. 3.12.). This can be seen as some manifestation that the physical 
scattering states are not free electrons but have always a Coulomb-field attached to them. 
The presence of a thermal bath perhaps influences somehow the effects of these unphysical 
degrees of freedom. 

Another point that is not yet clear is the process of 'thermalization': intuitively it should 
be the case that a particle immersed in a heath-bath looses momentum by repeated scattering 
processes; this is not the case in our framework. But to take intuition as a guide is not always 
the best thing to do. 

Some other ideas we did not investigate at all are the following: Implementation of the 
temperature in a different way, i.e. not as a scalar (cf. footnote g(|), different treatment of 
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the fermionic part of the theory, i.e. not through the implementation of the Fermi-Dirac- 
distribution in the theory in a way so similar to the photons or description of the theory 
using infraparticles (see e.g p6| and references therein) to mention just three. 

Finally, we mention that in the opinion of some authors (e.g. ]70|) a consistent pertur- 
bation theory for a QFT including temperature is not possible at all or at least not in such 
a direct way as known from the (T = 0)-case. We do not adopt this view, but the physical 
content and foundations of the theory have to be investigated and thought over very care- 
fully in order to set up a sensible thermal QED. Thus it could be that the theory is not good 
for arbitrary temperatures because the fermionic contribution from the heat-bath (virtual 
pair generation) is not incorporated correctly by treating it more or less analogously to the 
photon-part. 
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APPENDIX 



A Contractions, Propagators etc. 

A.l Contractions, Causal Functions and Feynman Propagators at 
T > in a>space 

In the following let A~ , "P - etc., and A + , etc., resp. denote the parts of A, ^ etc. (see 
(2.57), (2.63)) with absorption and emission operators, respectively. 
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Contractions: 



(2tt) 3 7 2\k\ 
-9p.v f d 3 k 



[4")( a; ),4+)(y)] = 
[4+)(x),i(-)(y)] 



(2tt) 3 
g^% p D^\x-y) 
-g^ f d 3 k 



(27r) 3 7 2|fc| 



+ e 



ik(x- 



:g^i^D { +\x-y) 

W / S(/(^- ife(x - !/) + (i + /WK fe(3: - y)N 
"Dfr'O* - id = ^n } (y - x) 

g^i f3 D[-\x -y) = -g^i PD ( +\y - x) 
g^i ^(x -y) = -g^i f>D%\y - x) 
g^i PD ( -\x -y) = -g^i D ( +\y - x) 
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{*(-)(x),*i +) (y)} = -i-3 / - ^ U-(p)(l-/-(p))(^ + mW 



ip(x-y) 



+ iV/+(p)(l - /+(P))G* - "») o6 e 



ip(x- 



=:i f3 S ( + ] {x-y) ab 



d 3 p 



iV/-(p)(l-/-(p))(^ + m) ai 



be 



^ + m 2 . 

- iy/f+(jp)(l-f+(pW - m) ab e^ 

{^(^^^(y)}-^ / 7= | f/-(p)(^ + m)„ t e-^-») + 

+ (l-/ + (p))G*-m) o6 e*<*-«'^ =: 
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+ (l-f+(pM-m) ba e- i * x - y)> j =: 

=:\ f3 s[-\y-x) ba 
{^-)( x )M+\y)} = j^ J 2 ^ P + m2 ( i V/-(p)(i - f-ipW + 



ip{x-y) 



- Wf+(p)( l - f+{p)){i> - m) ha e 



ip{x-y) 



--:i S[ 2 \y-x) ba 
d 3 p 



{&-)( X ),^\y)} = J— f (zV/-(p)(l -/-(?))(?> + m) 6a e'^-") + 

y^j J 2y<p 2 + m 2 \ 

+ V/ + (p)(l-/+(p))(jJ - m) fea e-^^-^ 



=: -i^fo-a;)*, 
{*(-)( 2; ),*[ +) (y)} = 7 ^ / - 7= = ((l-/-(p))(^ + m) bo e^) + 

+ f + (p)W-m) ba e- i X*-v) S J = . 



'22 VW x Jba 



Causal Functions: 
D = D(-) +£>(+>: 



?D 12 (x-y) = 0= P D 21 {x-y) 
D 22 (x-y) = - T =°D(x-y) 



e -ik(x-y) _ e ik{x-y) = T=0 £) ( a , _ y ) 



S=SH + g( + ) ; 

^(z - y) ah = / - 7^== (05 + mUe-^-ri + {j> - m) o6 e*<*-»> 



(2tt)3 7 2^+™ 2 
= T=0 S(x - y) ab 
S 12 (x - y)ab = = S 21 (x-y) ab 
% 2 (x-y) ab = T=0 S(x-y) ab 
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Feynman Propagators: 

D F _ £jret _ £,(-) _ jrjav _|_ £)(+). 



- i/) = ^3 y ^ e -*fc(»-") + /(A)(e- ifc <»-*> + e ik ^~ x ^ + T=0 D* ct (x - 
= ^o DF{x _ y) + _±_ J 0^ /(fc)(e -*(-,) +ett (-y)^ 

X 2 ^-2/) = -^2 ) (^-y) 

^ 1 (a:-i/) = -^r ) (*-») 

^40r-t/) = ^3 y 0^ fe (y-x) +/(A; )( e -^-x) +e ifc( y -x)^ _ T=0^rct (a ,_ 

= -( ^ F (, - „))• + ^JL_ I (/(fc)( C - tt <-«> + e*<— >)) 

gF _ g.(-) _ gret _ _g(+) _ S av : 

+ /_(p)C|* + m) ab e-^-y> - f + {p)(j> - mUe^-v^J - T =°S Ict (x - y)ab 

= T -°S F (x-yU + 
+ / - f P ( /_(p)0 + m) ob e-^(-^ - / + (p)0 - m) o6 e^-") 

(Z7Tj 7 2^p2 + TO 2 ^ 

% F 2 (x-y) afc = ^(x-yU 
% F 1 ( a; -y) a 6= ^(ar-yU 

% F (z - y ) afc = * / f P f (jJ + m)^-») + 
(27r) J 7 2VF^ + m 2 y 

+ / + (p)(^-m) ab e^-^-/_(p)(^ + m) afce -^-^^ - T =V ot (* - j/U 

= ( T =°S F (y-xUr + 

+ 7^3 / - J / + (p)(ii-m) a ^-'»-/-(p)(/ + mUe-^-») 

(27T)' 3 7 2^^ + m 2 \ 



A. 2 Relations 



"DtfW = - "Dtfi-z), ?D[- 2 \z) = ?d[ + 2 \-z) 
<>Dk\z) = D&(-z), <>d£\z) = D[V(-z) 
?d£\z)= D[+\z)= PD[t\-z) 

^ F (z)= ?D F n (-z), Dl 2 {z)= Dl 2 {-z) 
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T>0 b? 1 {z) = T>0 C 2 F 2 (,)= T>0 D[V(z)= T >°D[l\-z) 
^\z)=^\z), ^\z)=^\z) 

es[- 2 \z) = -es[t\z) 
% F 2 (z)= esb\z), = <fef a (z) 

T >°S?i(z) = T> "Sl 2 {z) = - T >°s[t\z) = T >°S ( h \z) 

A. 3 Covariant Notation and Momentum Space Representation 
Covariant Notation 

Pd£\x -y) = -^3 J d 4 fc«5(fc 2 )e-^-^(fc ) + f(k)) 

W+\x -y) = ^J d 4 fc«5(fc 2 )e-^)V/(fc)(l + /(fc)) 
< 3 D F u (x-y)= T=°D F (x-y) + -^_ j^e^m 



(>s£\x -y) = j±pj d%b{p 2 m 2 )8(p ) ■ 



^3 J d 4 pS(p 2 - m 2 ){p' + m)e 



ip(x-y) _ 



(27T) 

'f(y'o) - f-(impo)-f+(i>)e(-i>u) 

-ip(x-y) _ 



-y) = Jd i p5(p 2 -m 2 )(p' + m)e- 



■ (e(po)Vf- (p) (i - /- W) - (?(-po) v7f (p) (i - /+ (p)) 

• ( - 0(-po) + /_(p)0(po) + U(p)0(-po) 
^(x -y)= T =°S F (x -y)+*J d 4 P 6(p 2 m 2 ){j> + mje"*^-") • 

'/-(pMpo) + /+(p)0(-Po) 
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Momentum Space Representation 



-s(k")(e(ko) + f(k)) 



2n 



S(k 2 )y/f(k)(l + f(k)) 



1 



-1 



(2tt) 2 k 2 + iO 2ir 



6(k 2 )f(k) 



W(k) + -5{k 2 )f{k) 

Z7T 



T=0 D F 
-1 I 



(2tt) 2 -fc 2 + iO 2tt 



5(k 2 )f(k) 



T=0 S F (fcr + _L J(fe2)/(fc) 



^(p) = - m 2 ){i> + m) (d(p ) - f-(p)d(p ) - /+(p)0(-po) 

^ } (P) = - m 2 ){j> + m )U Po )^Up){l Up)) - 6(-p Wf + (p)(l f + (p)) 



^P(p) 
^i(p) 



— <5(p 2 - m 2 )^ + m) f - 6(-po) + /_ (p)6(p ) + f+(p)9(-p ) 



(2?r) 2 p 2 - m 2 + iO 2tt 



T=0 ^ F (p) + ^(p 2 - m 2 W + m) ( f-(p)6(p ) + f + (p)6(- Po ) 



B Three 'Intuitive' Observations in TFD 

In this appendix we present some unrigorous considerations, which may help to develop 
some intuition for TFD: 

Particle Number Operator: Calculate the expectation value of the ordinary (T = 0)- 
particle number operator for scalar bosons in the thermal vacuum: 



(n \a(k) + a{k)\np) 



(2.64) 



(v/TT7 <*(*, P) + - v 7 / a(*, &)) (^/^+7 a(*, /?) - v? a(*, /?) H 

(- v ^fa(fc,/3))(-V75(fc,/3) + )|^) = /W. 



<1> 



(B.l) 



= (0 



Here we made use of (flp\a(k, (3) + = and a(k, f3)\flg) — 0. So the thermal vacuum seems 
to mimic a black body heath bad as we expect it to do. 

Electron in a Heath Bath: Consider the following calculation (for the formulae with 



- ,G 



we refer to [E4| ) : 



b(p,(3)+\n p ) = e- lG b(p)+e lG \n p ) = e- lG b(p) + e lG e- lG \0) = 
= e-' G b(p) + \0) = b(p) + e- iG \0) = b(p)+\n p ). 



(B.2) 
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So it could be said that b(p,(3) + generates an electron with momentum p in the thermal 
setting. But this should not be taken as a sound physical interpretation of the formulae. 



A third observation: set T — in the formulae of subsection 2.2.2; then we get two copies 
of ordinary QED at T = 0. 

C Second Order Calculations 
C.l The £> 2 -Distribution 

D 2 (x,y) = -e 2 ^ b Y cd { 



:* Q (x)* 6 (x)* c 
{*<->(*),*<+> 

:^ a (x)^ h (x)^ c 
:^ a (x)^ b (x)^ c 

{¥-Hx),^ 

:* Q (x)* fc (:E)* c 

{*~\x),¥+^ 



y)} :^ a (x)t> d (y 
y)} :f a (x)f «,(!/ 
y)} :^ Q (rr)^ d (2/ 
y)} :* a (x)* d (?; 



\y)} j ( - 



(y)} :* 6 (x)* c (y): + 



(y)} : J 4 M (a ; )I,( y ):+[4-)^),i(+)(y)] 



(y)})(:A^x)A v (y):+[A(-\x),AM(y)] 



(x),* d + \y)} : § 6 (z)# c (j/): + 



(y)})( :A fl (x)A v (y):+[Ai-Hx),Al + \y)] 



e2 lablcd{ dito with x y 



+R' 



(C.l) 
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+ ^ :* o (a0*6(z)* c (l/)*«l(l/): +7 ^ ~ :* 6 (x)§ c (y): + 
+ - i P S [ +\x-y) bc :* a (x)* d (y): + 

+ 1 "Sfr>(i, - x) da j ^(x - y) b ^ ( ^(a^y): *£><+> (a - y) 

-^&a(x)*b{x)* c (y)h{v)-- +(-*) ^(l/-*)*. :* b (x)* c (y): + 
+ i^+)(a ; -y) bc :* (a:)f d (y): + 

+ H) ^ } (?y - *)dai S[t\x - y) b )j ( :^0r)i,a/): +$„,,» ^(a: - y)' 
- ^ :*„(a;)f 6 (ic)*c(l/)*d(l/): "sfc \v - x) da & h (x)V c (y): + 

+ i ^(y - x) da {-i) f>s£\x - y) b )j ( :4(a;)^(y): +y^i ^(a; - y)) + 

+ ^ :* a (x)^(x)* c (y)^(y): +1 S^\y-x) da :i b (x)* c (y): + 

+ i /3 4+ ) (x-y) bc :* a (x)|r d (y): + 
1 

+ \ S&{y - x) da l - ( 3 S ( +\x - y) b )j ( ^(^(y): (a; - y)) j 

f — e 2 7^ b 7^ d |dito with a; <-> y (and : a <-» c,b <-> <-» za)| J . 



+-R2 

(C.2) 

The four summands of the A^-part of this expression shall be denoted by A a , B a , C a ,D a , 
the respective four summands in the first factor of A a , D a by 1, 2, 3, 4, the respective two 
in the second factor of A a ,...,D a by 1,2. For the i?2~P ar t we choose the same notation 
with a superscript 13 instead of a . After carrying out the multiplications in this expression 
we then can refer to the several parts by A^jD^ etc. By Am, Dm we will denote 

Act 1 A.P p)a 1 
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After having done this calculation, we arrive at the following expression for D 2 : 

p D 2 {x,y) ='A + B + C + D' = 

4?i : -e 2 7M &a(x)* b (x)* c (y)* d (y):-A li (x)A v ( V ): 

AI2 ■ -e 2 7^7cd &a(x)* b (x)* c (y)* d (y): g^i ^d{+ ] (x - y) 

< : - (-e 2 7M, :*c(2/)* d (y)*a(x)* 6 (x): - x)) 
A a 21 : -e 2 ^^ &b{x)9 c {y): \ ^\y - x) da A^(x)A v (y): 

< ■■ (-e 2 7« :* d (y)*a(x): \ ^\x y) bc A v (y)A^x):) 
A% 2 : -e 2 l^ ab l u cd ^a(x)^ d (y): \ S[t\x - y) bc A^A^y): 

42 ■ (-e 2 7M, ■■VciyWx)-- \ P s[t\y - x) da A v (y)A,(x):) 

^1 : -e^Jcd &b{x)9 c {y): \ fj S ( h \y x) da9llv i ?d[+\x y) 

4i ■■ (-e 2 « 6 &d(v)*a{x): \ <>stf(x y) bc g, v i f>D<&\y x)) 

A% 2 ■■ -e 2 l^ cd :*a(x)* d (y): \ S[t\x y) bc g^i "£><+>(* y) 

A§2 ■■ (-e 2 7« ^c(yWx): \ ^{y - x) da9liv i ^[+\y - x)) 

At, : -e 2 r ab l v cd \ P S[:\y x) da \ <&<+>(* - y) bc A,{x)A v {y): 

< : - (-e 2 7M,T f3 S { n \x y) bc l - ^\y - x) da A„{y)A,{x):) 
A%i ■ -z 2 l" ab i: d \ ^u\y x) da \ S[t\x - y) bc g^i "£><+>(* y) 

A f L ■ (-e 2 7«7 S { u\x y) bc ) W+\y - x) da g, v i ^+\y x))+ 
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+e 2 ^ ab Y cd ^ a {x)^ b {x)^ c {y)^d{y)y.A li {x)A v { y y. 

+z 2 la b l V cd :*a(aO*b(aO*c(i/)^(y): <W - 2/) 

+ (-e 2 7 c ^ b ^(^(^(^(x): (y - x)) 

+e 2 r ab l v cd :Mx)*c(y): (-*) "S^fo " ^(*)^(»): 

+ ("^M, :*d(l/)*a(x): (-i) ^'(x - y) 6c :^(y)^(x):) 
+e^ b l v cd &a(x)h{v): i S[?(x - y) bc :A^x)A v {y): 

+ (-e 2 7^7 a M 6 ^c{y)h{x): i S[t\y - x) da :A v {y)A,{x):) 
+^ b lcd ^b(x)^ c (y): (-*) "S^fo - x) Aa9llv i PD { +\x - y) 

+ W d :* d (2/)*a(x): (-i) ^-'(x - 2/) bc5 ^ - x)) 

+^7^ :t«Wl d (!/): * ^+ } (x - yj^i ?d[+\x - y) 

+ (-e 2 7M, &c(vA(x): i S[? \y - x) da g^i - x)) 

+e 2 rM-i) s£\y - x) da i Ps[tHx - y) bc :A^x)A u {y): 

+ - v)bJ P S[t\y x)da :A v {y)A,(x):) 

+e 2 r ab l v cd {-i) s£\y - x) da i Os[tHx - y) bc g^i d[+\x - y) 

+ - I/W S[p(y - x) da g^i (>D[+\y x 

+e 2 ^ cd ^a(x)^ b (x)^ c (y)^ d (y)::A^x)A,( y ): 

- (e S )X :*c(y)* d (y)*a(x)* b (x)::i 1 ,(y) J 4 A1 (x):) 

" W :*c(l/)*d(l/)*a(*)*b(aO: - x)) 

+e 2 7 Q M 6 7cd :^W*c(l/): i ^2 - ^O^M 

~ (eSM, ^(f)*-^ 1 ^(x - !/)6c :l,(y)^(x):) 
+^7^ :*a(*)*d(l/): (-*) ^^(x - y) bc ^(x)^): 

" (eSM, :*c(2/)* b (x): (-i) ^(j, - x) da :i,(j/)^(x):) 

+e 2 7^7c d :*6(a;)*c(l/): i ^2 ^ - *)<ia<W - V) 

~ (e 2 7M, :*d(y)*a(x): i ^(x - "D^fo - x)) 

+^7^7^ :*a(x)* d (y): (-*) ^+ } (x - y) bc g^i (>d[+\x - y) 

~ (eSM, :*c(l/)*6(«): H) ^(l/ - *)<ia<W ^(l/ - *)) 

+e 2 7» ^2 } (y - »)«J-(-0 ^(x - v)bc M^My)'- 

~ (eSM ^(x - 2/)6c(-i) "S^fo - :i,(2/)^(x):) 

+e 2 7» ^2 } (y - x) da (-i) s£\x - y) bc g^i ^d[+\x - y) 

- (eSM ^2 } (* - l/)ic(-i) - *)d«9v>i D[V(y x)) 



78 



2 ^ b Y cd ^a(x)M^c(y)^d(y)::A^x)A u (y): 

2 < b lcd :* a (x)* b (x)* c (y)f d (y): g^i ^\x - y) 

- ("^ :* c (y)* d (y)* a (x)* b (x): <^ - x)) 
Hrfcd Mx)*M: \ fj S { 2 2 ] (y ~ x) da :A ll {x)A v {y): 

{-eHdl^b Mv)*«(*)- \ ?S { -\x - y) bc :A v (y)A^x):) 
2 M :*a(x)* d (y): \ \x y) bc M^Mv)'- 

(-e 2 lc d i: b ^c(y)k(x): \ ?Si?(y-x) da :A v (y)A,(x):) 
tabled ^b(x)*c(y>. \ P S^\y - x) da g^i Pd£\x - y) 

(-eHd^b ^d(y)*a(x): \ ^~\x - y) bc g^i *£><+> (j, x)) 
2 ^cd ^a{x)^ d {y): \ Sit\x y) bc g^i (>D<&\x - y) 

- (-<?l v c dll b ■■*c(yA(x): \ f3 s£\y - x) da g^i ^+\y - x)) 
2 lablcd\ P S£( V - x) da l - y) bc :A,(x)A u (y):) 

(-e 3 ^ "S^ \x - y) bc ± f3 S { +\y x) da :A l ,(y)A,(x):) 
2 <blcd\ P S^(y x) da l - Ps&\x - y) bc g^i D<£\x y) 

- i-e 2 lc d lZb- ^(x - y) bc - f>S&\y x) da g, v i f>D%\y x)). 



79 



C.2 Proof of the Causal Support of D 2 

Pay attention to :* d * a := - :^ a ^d-- 

An -0 '■ causal 

Ai 2 : D^\x -y)- D[+\y - x) = D^\x - y) + D { ~ ] {x - y) = 
= D n (x -y)= T=0 D(x - y) : causal 
A 2 i + A 3 i :1. + 2.term : S[i\y - x) da + S { +\y - x) da = 

= p S n (y - x) da = T=a S da (y - x) : causal 
2. + l.term : ^s[i\x — y)bc + S[i\x — y)b c ■ analogously — > causal 
A22 + A 32 :1. + 2.tcrm : S[i\y - x) da D{+\x - y) + S[\\y - x) da D { +\y - x) = 
= Ps[i\y - x) da t>D£\x -y)+ s[l\y - x) da D^\x - y) + 

+ f>s[t\v - *)da f3 D { +\y -x)- S[t\y - x)da f3 D [ +\x - y) = 
= Sn{y-x) d jD^\x-y)- 

- f>s[t\v - *)da (- D[t\y -x)+ d[+\x - y)) : causal 

V v ' 

=D(x-y) 

2. + l.term : S[t\x - y) bc D[+\x - y) + S[i\x - y) bc D { + ] \y - x) : 
analogously — > causal 

An : s[i\y - x) da s[t\x - y) bc - s\i\x - y) bc s\t\y - x) da = 

= s[i\y - x) da S[+\ X - y)bc + s[i\y - x) da s\i\x - y) bc - 

- s[i \x - y) hc s[+\y - x) da - S\i\y - x) da s[i \x - y) bc = 
= S[i\y-x) da Sn(x-y) bc - 

- s[i \x - y) bc [ s[t\y - x) da + s[i\y - x) da ] = 
= S[i\y-x) da Sn(x-y) bc - 

- 0s u\ x - V)bc Su{y - x) da : causal 
A 42 : s[-\y - x) da S { + ] \x - y) bc d[+\x - y) - 

- s{-\x - y) bc S ( +\y - x) da D { + ] (y - x) = 
= f>s[? (y - x) da S^ (x - y) bc (x - y) + 

+ S[-\y - x) da s[-\x - y) bc D[+\x - y) - 

- S[i\x - y) bc S\t\y - x) da D^\y - x) - 

- S[i\y - x) da s[i\x - y) bc d[+\x - y) 
the first two terms are proportional to Sn(x — y) bc , 

the other terms are treated as the first one in A22 + A 32 — > causal 
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Bu + Cn :0 : causal 

B12 + C12 :1. + 2.tcrm : p D { + ] (x - y) - P D { + ] {y - x) = 

2. + l.tcrm : - } (y - x) + fa ~ 3/) = 

B21 + C31 :1. + 2.tcrm : - S^\y - x) da - S { 2 f> \y - x) da = 

2. + l.tcrm : S^\x - y) bc + s£\x - y) bc = 
B31 + C21 :1. + 2.tcrm : S{? (x - y) bc + S { 12 \x - y) bc = 

2. + l.tcrm : - S[? (y - x) da - S[ 2 \y - x) da = 
B22 + C32 :1. + 2.tcrm : - s£\y - x) da D[+\x - y) - S { +\y - x) da D { +\y - x) = 
= - D { +\x-y) S 2l {y-x) da = Q 
2. + l.tcrm : S^\x - y) bc D { +\y - x) + "S<+ ; \x - y) bc ^ (x - y) = 
= D[+\x-y) S 21 (x-y) bc = O 
B 32 + C22 :1. + 2.tcrm : s[? (x - y) bc D { £ (x - y) + ^ (x - y) bc D { £ (y - x) = 
= D<£\y-x) S 12 {x-y) bc = Q 
2. + l.tcrm : - s[? (y - x) da f3 D { + ] (y - x) - s[ 2 \y - x) da f3 D { + ] (x - y) = 
= - D[+\x-y)% 2 (y-x) da = O 
B41 + C41 :1. + 2.tcrm : - S^\y - x) da S{?(x - y) bc + 

+ s£\y-x) da S[ 2 \x-y) bc = 
= - S { +\x-y) bc S 21 {y-x) da = Q 
2. + l.tcrm : S^\x - y) bc S{?(y - x) da - S { +\x - y) bc S[ 2 \y - x) da = 
= - s£\x-y) bc S 12 (y-x) da = O 
B 42 + C 42 :1. + 2.tcrm : - ^ (y - x) da s[?(x - y) bc D ( +\x - y) + 

+ f>s£\v - x) da s[ 2 \x - y) bc D^\y - x) = 
= f>sk\y - x) da s[+\x - y) bc [- D { +\x - y) + D^\y - x)} = 
2. + l.tcrm : S^\x - y) bc S[?(y - x) da D { +\y - x) - 

- s£\x - y) bc S[ 2 \y - x) da d[+\x - y) = 
= "S<+> fa - V)bc (y -x) da [ D[t } (y-x)- D^ (x -y)} = 
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Dn :0 : causal 

D 12 : D<£\x - y) - (y - x) = D^\y - x) + D[-\y - x) = 

= D n (y -x)= T=0 D(y - x) : causal 
D 21 + D 31 :1. + 2.tcrm : - S^\y - x) da - S ( +\y - x) da = 
= - 0S 22(V -x)da = - T=0 S da (y - x) : causal 
2. + l.term : S 22 \x — y)bc + S 22 \x — y)b c : analogously — ► causal 
D22 + D 32 :1. + 2.tcrm : - S { 22 \y - x) da D<£> {x - y) - S { +\y - x) da D { +\y - x) = 
= - P S { 2 -\y - x) da d£\x - y) - S { +\y - x) da D<£\x - y) - 

- S ( +\y - x) da D^\y - x) + S { 2 ?(y - x) da d£\x - y) = 
= - S 22 {y-x) da D<£\x-y) + 

+ S { +\y - x) da [- D { +\y - x) + D { +\x - y)] : causal 

= T=«D(y-x) 

2. + l.term : S { 22 \x - y) bc D^\y - x) + 

+ S 22 \x — y)bc D 22 \x — y) : analogously — > causal 
D41 : - S { 22 \y - x) da - y) bc + S ( +\y - x) da S { 22 \x - y) bc = 

= - f3 S [ 22 \y - *)da Sit\x - y) bc - S { +\y - x) da S { +\x - y) bc + 
+ S ( +\y-x) da S { 22 -\x-y)b c + 
+ siP(y-x) da sit\x-y) bc = 
= - S 22 (y-x) da S ( 2 ?(x-y) bc + 

+ S ( 2 P(y - x) da [ S { 22 \x - y) bc + S { +\x - y)bc\ = 
= - S 22 {y-x) da s£\x-y)bc + 

+ siP{y~x) da S 22 (x-y) be ■ causal 
D 42 : - S { 22 \y - x) da S { +\x - y) bc D<£\x - y) + 

+ S { +\y x) da si~\x y) bc D%\y x) = 
= - S { 22 \y - x) da s£\x - y) bc D ( + } (x - y) - 

- s£\y - x) da S { 2 p(x - y) bc Di+\x -y) + 
+ S ( 2 p (y - x) da S& (x - y) bc f Hf (y-x) + 
+ Si+\y - x) da Si+\x - y) bc D ( +\x - y) 

the first two terms are proportional to S 22 (y — x) da , 

the other terms are treated as the first one in D 22 + D 32 — ► causal 

Thus we proved that D 2 has causal support, and the requirements for the inductive causal 
construction are met. In addition we can see that the sum of the B- and C-contributions to 
D 2 equals zero. 
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C.3 Some Calculations for T 2 
1. Calculation of u> VP 

The (T > 0)-parts of the vacuum polarization have u) < — 1, thus they can be split trivially. 
This will be shown in the following. 



/i 41 -f- a 41 . 



^(y - x)-f ^(x -V)- f) S { A\v - P S[~\x - y) 
Ps[l\y - x)-f ?S lx {x -y)- S n (y - x)^ ^{x - y) 
T =°s[-\y-x)Y T =%(^y)- 

- ™S n (y-x)>f T =°sl-\x-y) 
T >°s[-\y-x)r T =°S n (x-y)- 

S n (y-x)r T >°S^\x-y) 



T=0r 



: [(T = 0)-part] r 



T >°S[~\y-x)r T =°S n (x~y) 
T =°S n (y-x)^ T >°Stf(x-y) 



(C.4) 



We consider the singular order of the (T > 0)-parts in momentum space. According to 
( 2.36| ) we thus have to investigate the following expression (we suppress some unimportant 
factors) : 



iimp(«5)(n(f),0( P ))^(n o (p),0(p)), 



with 



d 4 qtr 



5{q 2 ~ to 2 ) 



+ injl^ii —j) + m) 
■ [f-(q)6(q ) + U(q)6(-q )} S((q - pf - m'Jaignfa, - Po)- 
Now we calculate the go-integral using 



5{q 2 - m 2 ) = 



+ ra l \ J 
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which gives us (we choose polar coordinates, x := cos(0)) 



lim p{8) J d p<fi(p) I drr j dx j d<fi 



1 r 2ir 



-1 Jo 



S S H y 6 Wr 2 + m 2 e^M-/^ + l 



2 



• sign(\/ f 2 + «t 2 — f-)tt^( — ?\/?+ mi + T ra '' 

o 2 o o 2o z 



^ l l 



<5 5" y <5 J Vr 2 +m 2 e^«l+/^ + 1 

c ^l - 2 P^ 1 *^ AT^ 2 t \p\ , P 2 

sign(— V? + m z F)7;°\~rv r + m + -i~ rx + 

o 2 o o 2o z 



where we have go = +Vr 2 + m 2 in the first summand and qo = — \/r 2 + m 2 in the second. 
With the remaining ^-distribution we calculate the x-integral: 



ki^J d4p J Q 



OO p27T 



dr I d/ ip) 



\l>\ Vr 2 ~+ 



( 



p\ %/r 2 +m 2 ti -(ui sin <p+u 2 cos <k)^\J -p 2 r 2 -m 2 p^42 fyPo ">A 2 +m 2 -^- - ^ (p %/r 2 +m 2 — §y ) | — 



e 

analogous terms 



-l 



In the first summand we have x = Pa ^^^ m - — • Because of — 1 < x < 1, this equation 
requires r to fulfil some additional condition. This claims r ~ i for 5 — > 0, and the following 
considerations are thus valid also for u = (1, 0). After investigation of the integrand we can 
verify that all assumptions of the theorem of Lebesgue are met and it thus applies^]. Thus 
we can perform the limit in the integrand and it follows that the whole expression converges 
to zero, for all p(d) ~ 5~ m . m £ N (because lim5^o(Tsr^gr — ) = Vp ^ 0)! Thus we have 

eT+l 

tj VP =—00, and the quasi asymptotics equals zero. 

2. On the Vacuum of Thermal QED in 2. Order 

We have to show that the vacuum expectation value of the vacuum graph is finite in the 
adiabatic limit. T hus w e consider the Fourier transform of Jd 4 xd 4 yT 2 G (x 1 y)g(x)g(y) (use 



T VG - SSD, cf. (p0[)), which gives 

- J dVr 2 VG Mg (p).9o(-p) 

For the B + C-terms this equals^ 
i 



d P9o(p)go(-p) / d qid q 2 tr 



7 M (^i +m)7 M (^ 2 +m) 



5{q\ - m 2 )5(ql - m 2 ) 



■ [0(91,0)^ - 6{-qi,o)Vl+\ [6{q 2fl )V^ - e(-q 2fl )V 7 2+\6{{q 1 - q 2 - ep) 2 ) y/f(l + /). 

Evaluation of the trace gives 8(2m 2 — qiq 2 ), and with S(q 2 — m 2 ) the qix,- and (^.o-hrtegrals 
are calculated, then we insert 1 = Jd 4 kd(qi — q 2 — ep— k) and evaluate the /c -integral with 

54 Some problems possibly could be caused by p £ {p £ M | p 2 = 0}. But since this is a set of measure 
zero with respect to the measure d 4 p this is not the case. 

55 Here ^T/2± := J f±(q 1/2 )(l - /±(</i/ 2 ))- 
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5(k 2 ) and the <i 3 (^-integral using 6(q\ —cfe — eff— k). Thus we arrive at terms of the following 
form: 



1 



d 4 pg Q (p)go{-p) / d 3 k / d i q 1 



[2m 2 - \/\qi\ 2 + m 2 y (q\ - ep- k) 2 + m 2 + q\ - eqip - q ± k\ 
1 =VT^v / 2^V/(fc)(l + /(fc)) 



■\/|(fi| 2 + m 2 y (qi — ep — k) 2 + m 2 



2 j_ m 2 _ ^1 ^ _ ( p — k) 2 + m 2 — epa — |fc|). 



Here the limit e — > can be performed without problems, the whole expression (without the 
factor \) converges to a ^-independent constant, which can be cancellated by a suitable 
normalization ternQ Thus we established the existence of the adiabatic limit of the vacuum 
expectation value for the B + C-part of the vacuum graph. 
As for the A + D-term: 



J d 4 q 1 d 4 q 2 tr[....]S(q 2 1 - m 2 )5(q 2 - m 2 ) • 

• %i,o) - /-(<7i)%i,o) - /+(gi)0(-?i,o) 



6(-92,o) + /-(<? 2 )%2,o) + /+(92)0(-<Z2,o) 



S((qi - qi -p?) 



%i,o - 92,0 - po) + f(qi -q2-p) 



Its (T = 0)-part is discussed in p2| , p. 263 (it equals d(p) as it is defined there) and does not 
cause any problems because of the factor 6(p 2 — Am 2 ). The temperature dependent part has 
the same '^-structure' as the term discussed above and thus converges in the adiabatic limit 
to a constant, which can be cancelled by suitable normalization as well. 

3. Calculation of w SE 

The (T > 0)-parts of the self energy have a singular order of to — — oo. This will be shown 
in the following. 



^22 + ^32 : 



f>Stf (y - x) (x-y)+ s[t } (y - x) ^ (y - x) 

%,{y - x) PD[+\x - y) + s£\y - x) ^D n (y - x) 
T =°S n (y - x) T =°D[j\x -y)+ T =°s[j\y - x) T =°D n (y - x) 
T =°S n (y - x) T >°D[t\x -y)+ T > S[f\y - x) T =°D n ( y - x) 



(T = 0)-part 



T>0c(+) 



T=0 



>S n (y-x) T >°D[i\x-y) 



<S£>(y-x)*=»D n (y-x) 

-i ret 



+ 



56 See the discussion in the main text after eq. (3.20) 
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We consider the singular order of the temperature dependent parts in momentum space, ac- 
cording to (2.36) we thus have to calculate the following limit (we disregard some unessential 
factors): 



p(8)(±(^)J(p)) = (t (p)J(p)}, 

S^Q 



with 



£(p) := J d\ T >°s{t\p - q) T - D n (q) ~ 

- Jd A q8((p - q) 2 - m 2 )7 !/ (^ ~i+ m)^ v ■ 

■ [ - f-(p-q)0(po - qo) - f+(p-q)9(qo - Po)] S(q 2 )sign(q ). 
We calculate the go-integral using 8(q 2 ), which gives terms of the following form: 

fmp(8) JdYm Jd 3 q ^-^+2m)f_(P s - q )e(^-q )5(^-2^\q}+2^\q)x~m% 

2 

With the remaining ^-distribution we calculate the x- integral and get (we use S(fa — 2^& \q\ + 
t-g-\q\x-m ) - 2\pM d \ x 2\pm >>■ 

_ / 0\rua-(ui sini^fua cos cf>)^ \j ~Ap 2 r 2 +2p 2 rn. 2 -5 2 m 4 +2^-r~4p Sm 2 r-^u- 3 5^ (-T+m^Zpor)!-/^^-^ -1 

Proceeding as above for the vacuum polarization we finally get oj se = — oo. 



4. Some Calculations for the Thermal Mass 

Co: We choose u = (1,0), p || fc3-axis. The fco-integration is performed using the 8- 
distribution. Then we write the integral in polar coordinates, the (^-integration is trivial, 
giving a factor 2n, and the ^-integration is straightforward and yields an integrand propor- 
tional to 




p2 _ m 2 _ 2p Q r + 2\p\rx + iO p 2 — m 2 + 2por + 2|p|rx + iO 



1 1 (p 2 - m 2 - 2p r + 2\p\r + i0\ 1 ^ fp 2 - m 2 + 2p r + 2\p\r + i0\ 
2\p\r \p 2 — m 2 — 2por — 2\p\r + iO J 2\p\r \p 2 — m 2 + 2p^r — 2\p\r + iO J 



This is well defined for all values of r and x as can be seen readily. If we now have p 2 — m 2 it 
can be verified using ln(a + ib) = ln(V a 2 + b 2 ) + iarg(a + ib) that the real and the imaginary 
parts are zero each. 

C Vi u- Because the parts proportional to (pu) and p 2 , respectively, in the integrals for 
C P ,C U in (3.31) have the same structure as Co, which was defined by the corresponding 
integral with m 2 in the numerator, these parts can be omitted for p 2 = m 2 . Solving the 
remaining integrals as above and additionally using J °° dr e ^_ x = gjjy, equation ( 3.32 ) can 
be verified. 



D Third Order Calculations 

For D3 the general calculation scheme is as follows: we write down the various involved 
products of Tfe's according to the causal approach. Then we perform the contractions which 
lead to contributions with external legs we are interested in and consider only these terms. 



<s(> 



D.l Calculations for the Vertex Graph 
1. Calculation of 



R! 3 {x u x 2 ,x 3 ) =J2T(Y,x 3 )T(X) = 
Pi 



= T 2 (xi,x 3 )Ti(x 2 ) +T 2 (x 2 ,x 3 )Ti(xi) + T 1 (x 3 )T 2 (x 1 ,x 2 ), 



A' 3 (x x ,x 2 ,x 3 ) =J2t(X)T(Y,x 3 ) 

Pi 



Ti(x 2 )T 2 (xi,x 3 ) + T 1 (xi)T 2 (x 2 ,x 3 ) + T 2 (x ll x 2 )T 1 (x 3 ) . 
-> , ' v v ' >■ , ' 

"11 A™ Avt* 



This is calculated using the following formulae (see (2.7)) 



Ti(x) = -Tx(x), T 2 (xi,x 2 ) = -T 2 (xi,x 2 ) + T 1 (x 1 )T 1 (x 2 ) + T 1 (x 2 )T 1 (x 1 ). 

The three terms of the second equation define the three parts of R' 33 : R 33 x + R' 33 2 + R' 33 3 , 
analogously for A' 33 . We take here just the terms with external fields we are interested in 
(cf. Sec. |3.2.1 ): In x 3 there shall be an 'ordinary field' A c l ^- t (x 3 ) 1 no 'tilde field'; in x\ shall 



be ^(xi) or ^(xi), in x 2 ^(x 2 ) or ^(a^). 
R' 31 (xi,x 2 ,x 3 ) = - T 2 (xi,x 3 )Ti(x 2 ) 



■ *( Xl )^ % F !(x! -x 3 ) (-l)7"*(i3) :: A^xx)A v (x 3 ) : ■ 

■ I : §(x 2 )7 Q *(a; 2 ) : A a {x 2 )- : ^>(x 2 )^(x 2 ) : A a {x 2 ) 
e 3 : ffojy P S ¥ 21 { Xl - x 3 )(-l) 7 v V(x 3 ) :: Ap(x x )A u {x z ) : ■ 

*(x 2 )7 Q *(a; 2 ) : A a (x 2 )- : *(x 2 h a *(x 2 ) : A a (x 2 ) 



After performing the contractions, we have 

R' 31 = - e 3 : VfaW ^(a* - 13)7" 7 ^(xa - x 2 )^(x 2 ) : i *£><+> (a* - x 2 )A v (x 3 ) 
+ e 3 : tffcijy "Sfifri - x 3 )-fi S[P(x 3 - x 2 )^(x 2 ) : i l3 D { +\x 1 - x 2 )A v (x 3 ) 
+ e 3 : *( Xl )r "S&fri - x 3 )Y\ P S{+\x 3 - x 2 )j^(x 2 ) : 1 fJ D { +\ Xl - x 2 )A v (x 3 ) 
*{ Xl )Y PSl x { Xl - x z )-f% p s{t\x 3 - x 2 ) 7 ^(x 2 ) : * D { +\ Xl - x 2 )A v {x 3 ). 



e 3 
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Now to R' 32 : 

R 32 (xi,x 2 ,x 3 ) = - T 2 {x 2 ,x 3 )T 1 (xi) = 



mCompton 

— 1 2 

-e 3 : *(x 3 ) 7 1/ % F i(x 3 - x 2 )-f*(x 2 ) :: A^{x 2 )A v {x 3 ) : 
: *(x 1 ) 7 a *( a ; 1 ) : A a ( Xl )- : *(x 1 ) 7 «f (xi) : A a ( Xl ) 

- e 3 : *(x 3 ) 7 " P S ¥ 12 {x 3 - x 2 ) 7 ^(x 2 ) :: i^a^a*) : 
: ^h^On) : A a (xi)- : *(x 1 ) 7 «f(x 1 ) : A^) 



After performing the contractions (pay attention to :^(x2)---^(xi):— > — :<I'(xi)....<I'(x 2 ): !), 
this gives 



#32 - 

+ e 3 : ^{ Xi )y\ S[-\x 1 - x 3 )Y S^{x 3 - x 2 )j^(x 2 ) : i D { +\x 2 - Xl )A v (x 3 ) 

- e 3 : *(xi)y(-») "S^O&i - ^3)7" ^fi(^3 - £2)7^2) : 1 d[+ } (x 2 - si)^(x 3 ) 
+ e 3 : §(xi) 7 ^ "sfr^i - £3)7" ^f 2 (x 3 - x 2 ) 7a J (x 2 ) : 1 D { +\x 2 - Xl )A v (x 3 ) 

- e 3 : *{ Xl )^{-i) S^\ Xl - x 3 )Y Sf 2 (x 3 - x 2 )^(x 2 ) : 1 D<£\x 2 - Xl )A v {x 3 ). 

Now we turn to R' 33 : 

R' 33 .i =-T 1 (x 3 )T 2 (x 1 ,x 2 ) = 

S v ' 

rr7M0llcr 

= - e 3 | : *(x 3 ) 7 Q *(x 3 ) : A Q (x 3 )- : *(x 3 ) 7 Q *(x 3 ) : i Q (x 3 )| • 

• : *(x 1 ) 7 q *(x 1 )*(x 2 ) 7q *(x 2 ) : ^(x! - x 2 ) 
+ e 3 | : *(x 3 ) 7 Q vl/(x 3 ) : A a (x 3 )- : *(x 3 ) 7 Q *(x 3 ) : i Q (x 3 ) J • 

• : *(x 1 ) 7 q *(x 1 )#(x 2 ) 7q I'(x 2 ) : D* 2 (x 2 - Xl ) 
+ e 3 | : *(x 3 ) 7 Q *(x 3 ) : A a (x 3 )- : *(x 3 ) 7 Q *(x 3 ) : i Q (x 3 ) J • 

• : *(x 1 ) 7 q *(x 1 )*(x 2 ) 7q *(x 2 ) : D F 12 {x 2 - an) 
- e 3 | : *(x 3 ) 7 Q *(x 3 ) : A a (x 3 )- : #(x 3 ) 7 Q *(x 3 ) : i Q (x 3 ) J • 

• : #(x 1 ) 7 Q l'(x 1 )#(x 2 )7 Q l'(x 2 ) : D F 22 { Xl - x 2 ). 
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Considering only the interesting terms here as well we get 

^33,1 = 

+ e 3 : ^zih" \ ^\ Xl - 13)7" j S{+\x 3 - x 2 ) 1 ^(x 2 ) : ^(x, - x 2 )A v {x 3 ) 

- e 3 : ^ih" j ^(xj - ^3)7^ (S3 - x 2 ) 7aj *(x 2 ) : "£)f 2 (x 2 - x{)A v (x 3 ) 

- e 3 : *(zih"(-i) S^\ Xl - x 3 )-y v ] S[t\x 3 - x 1 ) 1 ^{x 2 ) : D\ 2 {x 2 - x x )A v [x 3 ) 
+ e 3 : *( Xl )^(-i) r3 S { 2 -\ Xl - x 3 )Yi S[t\x 3 - x 2 )^(x 2 ) : D^ 2 ( Xl - x 2 )A v {x 3 ). 

Now to the next term: 



#33,2 =Ti(x 3 )Ti(a;i)T 1 (a:2) =ie| : xl/^)^*^) : A Q (x 3 )- : foh^fo) : A Q (x 3 ) 

■iej : *(a:i)7 a *(a;i) : A Q (an)- : *(x 1 ) 7 a *(x 1 ) : 
■ie| : *(x 2 ) 7 Q *(x 2 ) : A Q (x 2 )- : *(x 2 ) 7 «i(x 2 ) : i Q (x 2 ) J . 

Contracting as above and disregarding uninteresting external legs we have 

-^33,2 = 

+ ie 3 : *( Xl )<fj S[~\ Xl - x 3 ) 7 ^ S[+\x 3 - x 2 )^(x 2 ) : i d[+\ Xi - x 2 )A v (x 3 ) 

- ie 3 : *( Xl )r] f3 s[-\ Xl - x 3 ) 7 "i S{+\x 3 - x 2 )j^(x 2 ) : i f3 D { + } (x 2 - Xl )A v (x 3 ) 

- ie 3 : *( Xl )r(-i) S { 2 -\x 1 - x 3 ) 7 ^ S { +\x 3 - x 2 )^(x 2 ) : i D[+\x 2 - Xl )A v (x 3 ) 
+ ie 3 : *{xi)-f{-%) f3 S { 2 -\ Xl - X3) 7 "i S[t\x 3 - x 2 )j^(x 2 ) : i D { +\ Xl - x 2 )A v {x 3 ). 



For the last term of R' 33 we have 



#33,3 = T 1 (x 3 )T 1 (x 2 )T 1 (x 1 ) =ie{ : *(x 3 ) 7 Q *(x 3 ) : A a (x 3 )- : *(x 3 ) 7 a *(x 3 ) : i Q (x 3 )| 

■ie{ : *(x 2 ) 7 Q *(x 2 ) : A Q (x 2 )- : (ar 2 ) 7 Q tf (a: 2 ) : i Q (x 2 ) 
■ie( : *(xi)7 a *(a;i) : : <f (£i) 7 a <l(xi) : A a (xi) 



Again contracting and considering interesting external legs only (Pay again attention to 
:*(ar 2 )....*(ori):-> - :*(xi)....*(x 2 ):!) this leads to 

-R'33,3 

■ ze 3 : ^zih" j ^r } (*i - *3)7"t ^(xg - x 2 ) 7m v]/(x 2 ) : * ^(^ - ^i)^(x 3 ) 

- *e 3 : *( Xl )r(-i) S ( 2 i\ Xl - x 3 )-f\ S[t\x 3 - x 2 ) 7aj vI/(x 2 ) : 1 D[+\x 2 - Xl )A v (x 3 ) 

- ie 3 : *(x!) 7 ^ S { -\ Xl - x 3 )7"i ^+ } (x 3 - x 2 ) 7ai *(x 2 ) : i ^+ J (x 2 - Xl )A v (x 3 ) 

+ ie 3 : *(a; 1 )7''(-*) S^\ Xl - x 3 )^ i S[?(x 3 - x 2 ) 7 ^ (x 2 ) : 1 D { +\x 2 - Xl )A v {x 3 ). 
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Now we calculate A' 3 



3 in quite the same way: 



A' 31 (xi,x 2 ,x 3 ) = - T 1 (x 2 )T 2 (xi,X3) = 

s .. ' 



iCompton 
2 



= + e 3 | : *(x 2 ) 7 «*(x 2 ) : A a (x 2 )- : *(x 2 ) 7 tt *(x 2 ) : i Q (x 2 ) j ■ 

• : tfOnh" f) S! 1 (x 1 -x 3 )(-l) 1 ^(x 3 ) :: ^(xx^xg) : - 
- e 3 | : *(x 2 h a y(x 2 ) : A Q (x 2 )- : *(x 2 ) 7 «f(x 2 ) : i Q (x 2 ) J • 

• : Si 1 (x 1 -x 3 )(-l)^(x 3 ) :: A^ Xl )A v (x 3 ) : . 



Doing the contractions this results in 



31 



+ e 3 : 9( Xl )-f "Sf^xi - £3)7^ "S^fo - ^2)7^2) : * - x^fe) 

- e 3 : *(a;i)Y« ^(an - x 3 ) 7 ^ ^ 2 } (x 3 - x 2 ) 7m I<(x 2 ) : i D 2 + ] (x 2 - Xi)^(x 3 ) 

- e 3 : *(x!)r ^f 2 (x! - x 3 )Y\ P S[-\x 3 - x 2 )j^(x 2 ) : 1 D[+\x 2 - x 1 )A l/ (x 3 ) 
+ e 3 : ^{xx)-f f3 Sf 2 ( Xl - x 3 )Yi f3 S { 12 \x 3 - x 2 ) 7ai *(x 2 ) : 1 D { 2 \ ] {x 2 - x x )A v {x 3 ). 

The next term is 

A 32 (xi,x 2 ,x 3 ) = - Ti(x\)T 2 (x 2 ,x 3 ) = 



- e 3 : *(xx) 7 ^ S[+\ Xl - x 3 )Y "Sfi(z 3 - x 2 )j^(x 2 ) : 1 d[+\ Xi - x 2 )A v {x 3 ) 

+ e 3 : VixiWi-i) S { 2 t\xi - x 3 )Y ^(a* - x 2 ) 7ai *(x 2 ) : * ^(xi - x 2 )A v {x 3 ) 

- e 3 : *( Xl )r] S{+\ Xl - x 3 )r % F 2 (x 3 - x 2 ) 7At f (x 2 ) : * - x 2 )A,(x 3 ) 

+ e 3 : ^(xih^-i) ^(x! - x 3 )Y ^f 2 (x 3 - x 2 ) 7a J(x 2 ) : i D ( +\ Xl - x 2 )A u {x 3 ) 



iCompti 




: *(x 1 ) 7 Q *(x 1 ) : A a ( Xl )- : *(x 1 ) 7 Q *(x 1 ) : ^(x^ 



• : §(x 3 )^ ^n(x 3 ~ x 2 )>f*(x2) ■■: A^{x 2 )A v {x 3 ) : 
e 3 | : *(a;i)7 a *(a:i) : A a ( Xl )- : *(a:i)7 a *(a:i) : ia(zi)} ■ 

• *(x 3 ) 7 1/ /3 5f 2 (x 3 - x 2 ) 7 ^(x 2 ) :: i M (x 2 )A„(x 3 ) : , 



which gives after contracting: 
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Now to A 33 : 

A 33 1 = _ T 2 (x 1 ,x 2 )T 1 (x 3 ) = 



rriM0llcr 

— 1 2 

-e 3 : *(x 1 ) 7 a *(x 1 )*(x 2 ) 7ct *(x 2 ) : D^( Xl - x 2 ) 



: *(x 3 ) 7 Q *(x 3 ) : A Q (x 3 )- : *(x 3 ) 7 Q *(x 3 ) : A a (x 3 ) 
+ e 3 : ^{x 1 ) 1 a ^{x 1 )^{x 2 ) 1 J{x 2 ) : ^(x, - Xl ) ■ 

■ | : *(x 3 )7 Q *(x 3 ) : A a (x 3 )- : *(x 3 ) 7 Q *(x 3 ) : ia(x 3 )| 
+ e 3 : *(xi) 7 Q *(xi)§(x 2 ) 7ct *(x 2 ) : ^(x, - Xi) ■ 

■ | : *(x 3 ) 7 Q *(x 3 ) : A a (x 3 )- : *(x 3 ) 7 Q l-(x3) : A a (x 3 )} 

- e 3 : ^{xi^ix^Xiha^Xi) : "D&fci - a*) ■ 

• | : *(x3) 7 Q <J/(x 3 ) : A Q (x 3 )- : *(x 3 ) 7 Q i(x 3 ) : i Q (x 3 )}. 

Disregarding terms with uninteresting external legs and performing the contractions we 
arrive at 

4' 

+ e 3 : -ar 3 ) 7 ^ ^ (x 3 - x 2 ) 7ai *(x 2 ) : ^(xi - x 2 )A„(x 3 ) 

- e 3 : ^(xi - x 3 ) 7 »i S { 12 \x 3 - x 2 ) 7a J(x 2 ) : ^f 2 (x 2 - Xl )A v (x 3 ) 

- e 3 : - 0:3)7" j } (^3 - ^2)7^(^2) : P Dj 2 (x 2 - Xl )A v {x 3 ) 
+ e 3 : *(xi) 7 "(-i) ^(^1 -a= 3 )7^ ^ (a* - ^2)7^(^2) : - x 2 )A,(x 3 ). 

The next term is 

^33,2 = T^T^T^) =ie| : *(x 1 )7 Q *(^i) : 4(h)- : * (an^f (a*) : i a (xi)j • 

■ie{ : *(x 2 ) 7 Q *(x 2 ) : A a (x 2 )- : *(x 2 ) 7 Q *(x 2 ) : i Q (x 2 ) 



ie\ : *(x 3 ) 7 Q *(x 3 ) : A Q (x 3 )- : *(x 3 ) 7 Q *(x 3 ) : i Q (x 3 ) 



Once more we perform the procedure of contracting and disregarding uninteresting terms 
with the result 

A 1 

+ ie 3 : §(xx) 7 ^ f3 S ( +\x 1 - x 3 )y\ ^ (x 3 - x 2 )j^(x 2 ) : i P D { +\ Xl - x 2 )A v {x 3 ) 

- ie 3 : ^ih" j ^(X! - x 3 ) 7 ^ ^ 2 } (x 3 - x 2 ) 7ai 1>(x 2 ) : * ^+ } (x 2 - Xl )A v (x 3 ) 

- ie 3 : ttOnh^-i) ^(xi - x 3 ) 7 ^ "S^fo - ^2)7^(^2) : * - X!)4(x 3 ) 
+ ze 3 : ^ih^-i) ^(^i - ^3)7"* ^12 } (^3 - x 2 ) 7a J(x 2 ) : i "D&^i - x 2 )A„(x 3 ). 
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Finally 



^33,3 = T 1 (x 2 )T 1 (x 1 )T 1 (x 3 ) =ie\ : *(x 2 ) 7 Q *(x 2 ) : A a (x 2 )- : §(x 2 ) 7 Q §(x 2 ) : A a (x 2 ) 



iej : *(x 1 ) 7 Q *(x 1 ) : A a ( Xl )- : *(x 1 ) 7 Q *(x 1 ) : A a ( Xl ) 
iej : *(x 3 ) 7 Q vl/(x3) : A Q (ar 3 )- : §(£3)7^3) : A a (x 3 ) 



We calculate as above ( :&(x2)....&(xi):—* — :*(xi)....4'(x 2 ):!) This gives 



4' 



33,3 



+ ie 3 : *(xx) 7 ^ ^(ari - x 3 ) 7 v - S^\x 3 - a*^* (* 2 ) : i (} D { +\x 2 - Xl )A v (x 3 ) 
i i 

- ie 3 : *( Xl )r(-i) P S ( +\x 1 - x 3 )^- ^ (x 3 - x 2 )j^(x 2 ) : i D^ (x 2 - Xl )A v (x 3 ) 



— ie 



+ ie 3 : *(a; 1 )7' , (-i) S i +\x 1 - x 3 )>fi f3 S { u \x 3 - x 2 ) 7 ^(x 2 ) : * ^+ J (x 2 - x 1 )A v {x 3 ). 

Now we calculate D = R' — A' and order it according to the external legs. We thus arrive 
at the following expression (the arguments of 'SSD' are {x\ — x 3 ), (x 3 — x 2 ), {x\ — x 2 ), 
always in this order) 



Dj = +e 3 



ii / "-"n -^n 



ii 7 



+ e 3 : ^(xi)^. 



+ e 3 : *(xi) 7 ^- 



+ e 3 : *(a;i) 7 ' i - 



ii 

; /3 

'ii 

- <^+V ^f 2 + 

+ % a rr ^ ^ 2 } + 
-^rV^ )/, -Df 2 +' 

+^ 1 +y^ 2 ) ^f 2 - 

-psi^r f3 sif fJ D[V . 



J^{x 2 ) : A v {x 3 ) + 



7m *(x 2 ) : A v {x 3 ) + 



l^{x 2 ) : A v {x 3 ) + 



v "5f? D^ } 



21 

_ flc(-) , 

J 21 7 "12 -^22 

- "s&V estf vd\ 2 ^ 

+ ^7" ^12 > ^2 } H 



7^*(x 2 ) : A„(x 3 ) 
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2. Calculation of u; v 



We have to investigate (see (2.36)) 



with 

A(p, g) := /" d 4 fc 7 ^ ^ av (p - fc) 7 " #? ret (g - fc) 7M ^W(jfe) + 

where the stands for terms without essential differences in their structure with respect 
to the determination of the singular order. We choose the Lorentz-system wherein u = (1, 0) 
and use 5(k 2 ) to perform the fc-integration. This gives integrals of the following form: 



lim P (S) J d 4 pd 4 q4>(p, q)Jd 3 k^- [ ~ i 



\k\ 

7- and slash-factors 

p 2 2pk 2 q 2 2qk 2 p q 1 

<H77 J m F" - m W F + M^T" - feo)- 



"5 2 S 6 2 S & o "' e P\k\ _ i' 

We calculate the |fc|-integral with one of the ^-distributions and the a;- integral with the 



other. Proceeding as in the determination of lo (see App. C.3) we derive lj = — oo 



D.2 Calculations for the Third Order Self Energy 
1. Calculation of D| ord SE 

The procedure is essentially the same as in the corresponding calculation for the vertex. 
Nevertheless, we write it down here in detail to give a reference of the various signs, factors 
i, etc., which are crucial for the discussion of a possible cancellation of the IR-divergencies. 
There are two graphs that have to be calculated, denoted by 'A' and 'B', respectively, (see 



Sec. 3.2.2). As in the calculation for the vertex we have: 



R' 3 (x 1 ,x 2 ,x 3 ) =YT(Y,x 3 )f(X) = -T 2 (xi > x 3 )T 1 {x2)-T2(x2 ) x 3 )T 1 (xi)- 

' ' * ' 

7?' 7?' 

^31 ^32 

- T 1 (x 3 )T 2 (xx,X2) + T 1 (x 3 )T 1 (xi)T 1 (x 2 ) + Ti(x 3 )T 1 (x 2 )T 1 {x 1 ), 



A' 3 (x-l,X2,x 3 ) = y^f(X)T(Y,x 3 ) = -Ti(x 2 )T 2 (xi,x 3 ) -T 2 (xi)T 2 (x2,x 3 ) - 
p . ' v . ' 

2 A' A' 

^31 ^32 

- T 2 (xi,X2)T 1 (x 3 )+Ti(xi)Ti(x2)Ti(x 3 ) +Ti(x2)Ti(a;i)ri(x3) . 
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First we consider graph B (disregarding terms with external legs different from ^(^3), A v {x 3 ) 
and &(xi) or 

R31 = - T 2 (xi 1 x 3 ,) Ti{x 2 ) = 

trivial graph 

+ ie 3 : ^(x 3 )Y^(x 3 )^{xx)j^(xx) :: A v {x 3 )A ll {xx) ■ ■ 

*(x 2 )j a V{x 2 ) : A a {x 2 )- : ^!{x 2 )^{x 2 ) : A a (x 2 ) 
ie 3 : *(x 3 )7 l/ *(x 3 )*(x 1 ) 7 ^*( a ; 1 ) :: A u {x 3 )A^{xx) : ■ 



1-' 

*(a: 2 ) 7 Q *(x 2 ) : A a (x 2 )- : *(x 2 ) 7 Q *(x 2 ) : A a (x 2 ) 



After having performed the contractions providing us with the external legs we are interested 
in, we have 

R 31 = 

- ie 3 : * (3:3)7" j s[t\x 3 - x 2 ) 7/1 ~ S { n\x 2 - x 1 ) 7 ^( Xl ) : A u {x 3 )i p d[+\ Xi - x 2 ) 
+ ie 3 : VfaWi s[?(x 3 - x 2 )^{-i) S^\x 2 - xx)l^{xx) : A l/ (x 3 )i d[+\xx - x 2 ) 
+ ie 3 : §(3:3)7" j S[+\x 3 - x 2 )^i fs s[ 2 \x 2 - xx)^{xx) : A v (x 3 )i d£\xx - x 2 ) 

- ie 3 : VfaWi P s[p(x 3 - x 2 )^ P S { 22 \x 2 - ii) 7 "*(a;i) : A„{x 3 )i d£\xx - x 2 ). 



The next term, R' 32 — —T 2 (x 2 ,x 3 )T 1 (xi), is calculated as li^j Brtex (see App. D.l): 

Compton 

R 3 2 =-e 3 : *(x 3 h» "Sf^a* " x 2 )r*(x 2 ) :: ^(a^A,^) : • 

*(x 1 ) 7 a *(a: 1 : A Q (x 1 )- : * (x 1 ) 7 Q! f (xj : i Q (xx) 



3 3 : §(x 3 ) 7 " Sf 2 (x 3 - x 2 )^(x 2 ) :: A^A^) : • 



*(x 1 ) 7 tt *( a ; 1 : A a ( Xl )- : *(x 1 ) 7 Q *(x 1 : A a ( Xl ) 



Performing the relevant contractions gives 



R ' 3 2 = 



- e 3 : §(x 3 ) 7 " Xi(^3 - x 2 ) 7 ^ S { +\x 2 - xi)^(x x ) : A^)* ^(^ - n) 

+ e 3 : §(z 3 ) 7 " "Sfifo - sa)V» S { +\x 2 - xxh^(xx) : A,(z 3 )i "M^fo - a*) 

§(x 3 ) 7 " ^(^3 - z 2 ) 7 M H) ^(^ - ^1)7^1) : A v (x 3 )i f3 D { +\x 2 - Xl ) 



e 3 



+ e 3 : §(x 3 ) 7 " Sf 2 (x 3 - z 2 ) 7 ^ ^(sa - x 1 ) 7 ^(x 1 ) : ^(13)1 ^(x 2 - Xl ). 
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The next term is 
#33,1 = - Ti(x 3 )T 2 (x 1 ,x 2 ) = 



SE 

„3 / 



- ie*^ : *(x 3 )7 Q *(x 3 ) : A a (x 3 )- : *(x 3 ) 7 a *(x 3 ) : A a (x 3 ) 

■ : §(x 2 ) 7 " S^{x 2 - *i) D^(x 2 - xjn* ( Xl ) : 

■ : *(x 2 ) 7 " S^(x 2 - D\ 2 {x 2 - x 1 ) 7 „*(x 1 ) : 

- ie 3 ^ : ^(x 3 ) 7 Q *(x 3 ) : A a (x 3 )- : (s 3 )7 a tf (a*) : i Q (x 3 )^ • 

• : *(x 2 ) 7 "(-l) Sf 2 (x 2 - xi) D\ 2 {x 2 - x 1 ) 7 ,*(x 1 ) : 

- ie 3 ( : *(x 3 ) 7 Q vI/( X 3) : 4(x 3 )- : *(x 3 ) 7 Q *(x 3 ) : A a (x 3 ) 

■ ■ * (a^V Sl 2 {x 2 - xi) Dl 2 {x 2 - Xl ) lv ^{ Xl ) : , 
which after performing the relevant contractions leads to 

#33,1 = 

- ie 3 : *(x 3 ) 7 ^ S[+\x 3 - x 2 )-f SUx 2 - x 1 ) 1 ^(x 1 ) : A v {x 3 ) D^(x 2 - Xl ) 
+ ie 3 : *{x 3 )-fi S[P(x 3 - x 2 )^ S^(x 2 - x^j^ix,) : A v {x 3 ) D\ 2 {x 2 - Xl ) 

- ie 3 : 9(x 3 )-fj S[+\x 3 - x 2 )Y Sf 2 (x 2 - x 1 ) 1 ^(x 1 ) : A„{x 3 ) D F l2 {x 2 - Xl ) 

- ie 3 : 9(x 3 )-fi S{+\x 3 - x 2 )^ ^ 2 F 2 (x 2 - x 1 ) 1 ^(x 1 ) : A v {x 3 ) D\ 2 {x 2 - Xl ). 
Now we come to 



#33,2 = T 1 (x 3 )T 1 (x 1 )T 1 (x 2 ) =»e| : §(x 3 ) 7 Q *(x 3 ) : A a (x 3 )- : *(x 3 ) 7 «*(x 3 ) : A a (x 3 ) 

■ie{ : *(x!) 7 a *(xi) : A a ( Xl )- : *( Xl )^(xi) ■ 4(zi)j 
■ie{ : *(x 2 ) 7 Q *(x 2 ) : A Q (x 2 )- : *(x 2 ) 7 «*(x 2 ) : i Q (x 2 ) J 



which gives after having done the contractions 

#33,2 = 

+ ie 3 : *(x 3 )7" j ^+ } (x 3 - x 2 )r] S { ^\x 2 - anh^fa) : 4,(3:3)1 Mt^i - **) 

- ie 3 : *(x 3 ) 7 ^ "sfcfta* - x 2 )^{-i) S { 2 ^\x 2 - xi) 7 ^(xi) : A v {x 3 )i f3 D { + } (xi - x 2 ) 

- ie 3 : §(3:3)7" j ^+ } (x 3 - x 2 )-fi f3 s[-\x 2 - xih^i) : A v {x 3 )i D ( 2 +\x 1 - x 2 ) 

+ ie 3 : *(x z )-fi S[t\x 3 - x 2 )-f\ S { 22 \x 2 - Xl )^( Xl ) : A„(x 3 )i D { + ] (*i - 12). 



i? 33 3 is calculated analogously, the result is the same as for R' 33 2 but with xi and x 2 ex- 
changed - before performing the contractions! After having done them, we have (disregarding 
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terms with uninteresting external legs): 

-^33,3 = 

- ie 3 : *(x 3 h v ^ P S[t\x 3 - x 2 )r \ S[t\x 2 - x 1)7^1) : A v {x 3 )i D[+\x 2 - Xl ) 

+ ie 3 : tffoh'j P S[t\x 3 - x 2 )ri s£\x 2 - Zi) 7 ^(zi) : A v (x 3 )i D[+\x 2 - Xl ) 
+ ie 3 : *(x 3 ) 7 ^ P S[t\x 3 - x 2 )^{-i) S [ +\x 2 - ^1)7^(^1) : A v {x 3 )i D { +\x 2 - Xl ) 

- ie 3 : *(x 3 ) 7 ^ "^(a* - x 2 ) 7 ^ s£\x 2 - x 1 ) 7 ^(x 1 ) : A„(x 3 )z D { +\x 2 - Xl ). 

Now we come to A' 3 : 

A' 31 = — Ti{x 2 )T 2 (x\,x 3 ) is calculated just as R' 31 , but with the order of the factors reversed. 
With the relevant contractions this gives 

A' - 

- ie 3 : §(x 3 ) 7 ^ f3 S { n \x 3 - x 2 )j^ S[+\x 2 - x 1 ) 7 ^( Xl ) : A 1/ (x 3 )i D { +\x 2 - Xl ) 
+ ie 3 : *(x 3 ) 7 ^ P S { -\x 3 - x 2 )^{-i) S { +\x 2 - £1)7^1) : A v {x 3 )i D 2 +\x 2 - Xl ) 
+ ie 3 : *(x 3 ) 7 ^ f3 S { n \x 3 - x 2 ) lfi i S { +\x 2 - x 1 ) 7 ^(x 1 ) : A v {x 3 )i d[+\x 2 - Xl ) 

- ie 3 : VfaWi S[ 2 \x 3 - x 2 )j^ S { 2 p(x 2 - x 1 ) 7 ^(x 1 ) : A v {x 3 )i D { +\x 2 - Xl ). 
The next term is A 32 = —Ti(xi) T 2 (x 2 , x 3 ). Performing the relevant contractions this leads 

S v ' 

Compton 

to (pay attention to the fact that :<]>.. - :*...*:): 

A' - 
^32 

■ e 3 : *(x 3 ) 7 " S^(x 3 - x 2 ) lt X fJ S { n \x 2 - x 1 ) 7 ^(x 1 ) : A„{x 3 )i d[+\ Xi - x 2 ) 

- e 3 : 9(x 3 )>f ^(xs - x 2 )^i S[ 2 \x 2 - x 1 ) 7 ^(x 1 ) : A v (x 3 )i D { +\ Xl - x 2 ) 

+ e 3 : H>(x 3 W Sf 2 (x 3 - x 2 )^(-i) S { ~\x 2 - x 1 ) 1 ^(x 1 ) : A v {x 3 )i d[+\ Xi - x 2 ) 

- e 3 : 9(x 3 )-f Sl 2 (x 3 - x 2 )ljr S ( 22 \x 2 - x 1 ) 1 ^(x 1 ) : A u {x 3 )i D 2 +\ Xl - x 2 ). 

Now we turn to 

A' 33 1 = —T 2 {x 2l x 3 ) Ti(x 3 ), which gives with the relevant contractions 

V v ' 

SE 

A' 

+ ie 3 : §(ar 3 ) 7 ^ S[^ (x 3 - x 2 )^ S^(x 2 - xi)-f*(xi) : A v {x 3 ) f3 D F n (x 2 ~ Xl ) 
- ie 3 : *(ar 3 )7^ S { 12 \x 3 - x 2 ) ll± S^(x 2 - x 1 ) 7 ^( Xl ) : A v (x 3 ) D\ 2 {x 2 - x{) 
+ ie 3 : *(ar 3 ) 7 ^ S[^\x 3 - x 2 )^ S ¥ 12 (x 2 - x 1 ) 7 ^(x 1 ) : A v {x 3 ) Df 2 (x 2 - Xl ) 
+ ie 3 : V(x 3 )Yi P s[ 2 \x 3 - x 2 )^ Si 2 (x 2 - x 1 ) 7 ^l-(x 1 ) : A v {x 3 ) D F 22 {x 2 - Xl ). 
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The next term is 

A' 33 , 2 — Ti(a;i)Ti(a:2)Ti(a;3) 

A 1 

^33,2 — 



(after having done the relevant contractions) 



- ie 3 : $(i s h"j S^\x 3 - x 2 ) 7 ^ S[^\x 2 - anh^fri) : A,(x 3 )» *£><+> (an - a*) 
+ ie 3 : *(x 3 ) 7 ^ S[-\x 3 - x^-i) S { ~\x 2 - Zih^i) : A,(x 3 )z "D&'fa - x 2 ) 
+ ie 3 : *(x 3 ) 7 ^ S^\x 3 - x 2 )^i S ( V2 \x 2 - x x )^{xx) : A v {x 3 )i D { +\ Xl - x 2 ) 

- ie 3 : VfaWi S[ 2 \x 3 - x 2 )^ ^ {x 2 - x 1 ) 7 ^(x 1 ) : A v {x 3 )i D { +\ Xl - x 2 ). 
Finally we have to calculate 

A333 = Ti(x 2 )T 1 {x\)Ti(x 3 ) =>■ (proceeding as above) 



4' 



33,3 



+ ie 3 : *{x 3 )Y\ S^\x 3 - z 2 ) 7 ^ f3 S { + ) (x 2 - ^h^O : A v (x 3 )i D { +\x 2 - Xl ) 

- ie 3 : *(x 3 ) 7 ^ S[-\x 3 - x 2 ) lfl i S[?(x 2 - x 1 ) 1 ^{x l ) : A v (x 3 )i D[+\x 2 - Xl ) 

- ie 3 : VfaWi fj s{-\x 3 - x 2 ) lfl {-i) S { +\x 2 - Xl )^{ Xl ) : A v {x 3 )i D ( +\x 2 - Xl ) 
+ ie 3 : ttfob"* S[-\x 3 - x 2 )j^ S { +\x 2 - x 1 ) 7 ^(x 1 ) : A v {x 3 )i D { +\x 2 - Xl ). 

Thus we know A' and R' from the third order self energy graph B. This gives the following 
result for D = R! — A' of graph B (the arguments of '55' are (x 3 — x 2 ), (x 2 — x\), always 
in this order): 



D 



3.ord.SE.B 



e 3 : V(x 3 ) 



-r s^ sit ] i, DiV 



D n 7 j ii7m -^11 



12 



-7^+V^+VM! 

-r s[p^ s { 2 V^ D 2 i 



+r s[i ) r s[Vi» D[t 



-), 



-7^5f 2 7^5<7 V ^ 

~r s[ 2 ] r ^ 2 F i7, ^f 2 
+r s[ 2 ) r s£h, D[t 



(xi - x 2 )- 
x 2 - x\)- 
x 2 - xi)- 
x 2 - xi)- 
x 2 - x\)+ 
(xi - x 2 ) + 
(xi - x 2 )- 

(x 2 - Xl)- 

(x 2 - xi)+ 
(x 2 - Xi)+ 
(x 2 -xi)- 
x\ - x 2 )- 
x\ - x 2 )- 
x 2 - x-t)- 
x 2 - xi)+ 
(xi - x 2 ) + 
(xi - x 2 )- 
(x 2 - Xi)- 



f p s (-) r s (+) 1 ^ D (+) {x2 _ xl) 



: A v (x 3 ) 
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+e 3 : #(ar 3 ) 



_ r /j 5 (-) 7m _ X2)+ - 

+7 ^ 5 (+) 7A1 ^(-) 7ai ^(+) (xi _ X2) _ 

-7" "s£V ^f 2 7M "Dfate - zi)+ 
+7" "S&V PS&ir Dl 2 {x 2 - Xl )+ 



"7 



_ 7- ^12 
+r/35 (+) r ^(+ 

7^ ^2> 7 



7 M ^t ) (^i-^)- 
7 M ^+ ) (x 1 -x 2 )- 
1 ^ r3 D[1 2 \x 2 -x 1 ) + 
1 ^D ( 2 1 2 \x 2 -x 1 ) + 
1 ^D ( i ^\x 2 ^x 1 y 
1 ^D i £{x 2 -x 1 )- 



-Y f> Sl ll ^s[ 2 ) 1 ^D ( 2 X\x 1 -x 2 )+ 
+r S F 2 7M _ X2 )_ 

_ 7 , ^^F (a . 2 _ + 

+7 , S (-) r S T 2lii D F 2{X2 _ Xi) + 

+r "s<:V "sfr^ (>Dlt\xi x 2 y 

_ r p s (~) r S (-)^ D (+) {X1 _ X2 )_ 

- 7 " ^rV - 
.+r^ ) 7"^ ) 7/4 ) (^-^). 

Now we simplify this expression; in the first part, the ten terms with a type f -vertex in x 2 
give after suitable modification (consider the support properties of the quantities involved) 



tf(sci) : A v (x 3 ). 



r-,3.ord.SE,B;2:2:typol _ 
3,** 


% a i> 3 


- Xi)^ p S u (x 2 - 


x 1 ) 1 ^D[-[\x 1 - 


- x 2 ) 


+ 


f3 S\f{xz 


-x 2 ) 1 ^s[-\x 2 




- x 2 ) 






- x 2 )r ^(x, - 


X1 ) ltl f>D$( xl - 


x 2 ) 




%i(xs 


-x 2 )^%f(x 2 - 


-x l ) 1 ^D { ~\x 1 


- x 2 ) = 


= + 




-x 2 )-f PSH{x 2 - 


x l ) 1 ^D i il \x l 


- x 2 ) 




S{f(x 3 


-x 2 ) 1 ^S i n \x 2 


-x^^DlKx, 


- x 2 ) 


+ 




-x 2 ) 1 ^S ( - 1 \x 2 


-x^^Dlfix, 


- x 2 ) 




%?(x 3 


-x 2 )^%f(x 2 - 


-x 1 ) 1 ^d[~\x 1 


- x 2 ). 



The terms with type 2 in x 2 give 



D. 



3.ord.SE,B;z 2 :type2 
3,** 



o. 



because both D 12 and S\ 2 are equal to zero. Likewise we have 

y-)3.ord.SE,B;a;2:typcl q an( ^ £j3.ord.SE,B;iE2:typc2 _ g 

3,** 3,** 



D 3 for graph A is calculated in the same way. 



2. Singular Order 

The structure of the third order SE graphs is S(x 3 — x 2 )S(x 2 — x{)D{x\ — x 2 ). This gives 
after Fourier transformation S(q) jd A kS(q + k)D(k), which means that only 

w 2.ord.SE ig 

relevant, because the first factor (S) has singular order ui < — 1 and the second one is 
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exactly the second order SE. Since the temperature dependent parts of the second order SE 
have singular order oj < — 1 the temperature dependent parts of the third order SE can be 
split trivially. 



D.3 Calculations for the Thermal Corrections to fj, e - 

In this appendix we present the calculations to solve the linear system (3.63) - ( p. 67 ), which 
reads in matrix form 

/o\ 

h 
h 
h 



2m 2 



n 4 + (pq) 2 
2m 2 (pq) 



2( P q) 
2m 2 (pq) 
n 4 + {pq) 2 



2(pu + qu) 
2(pu)(m? + pq) 
(m 2 +pq)(pu + qu) 



/4 

m 2 

pq 

pu m 2 (pu) + (pq)(qu) m 2 (qu) + (pq)(pu) m 2 + (pu) 2 + (pq) + (pu)(qu) 
\1 (pu) 2 + (qu) 2 2(pu)(qu) 2(pu + qu) 

(A\ 
B 
C 
D 
\G) 



1 \ 

(pu) 2 
(pu)(qu) 
pu 

1 / 



(D.l) 



Before solving that, we calculate the four integrals I2, ■■■,h, as defined in ( 3.63 ) - ( 3.67 ) 

1 pk 



d 4 kS(k 2 



( e /3|M _ 1) (-2qk + i0)' 



(D.2) 



We choose the Lorenz system with u = (1, 0) and q — (qo, 0, 0, |g|). Then we perform the fco- 
integration with the ^-distribution and choose polar coordinates for the spatial integration 
(x := cos(9)). This gives 



dr- 



d.r 



dip- 



1 



P\ku\ _ I 



Por — pir sin((/j)-\/l — x 2 — p2r cos(tp)Vl — x 2 — pzrx 
-2q r + 2\q\rx + iO 



+ 



~Po r ~ Pi r sm(ip)yl — x 2 — p2r cos(</3)Vl — x 2 — p^rx 



2q r + 2\q\rx + iO 



(D.3) 



Now we perform the (^-integration, which leads to a solvable r- integral (use J^°dr e p^_ 1 
^■).Thus we have 



12/3 2 J-! 



dx 



Po ~ P3X 



-PO -P3% 



-qo + \q\x + iO q o + \q\x + i0 



(D.4) 



In the second term we transform x — » — x and see that in the sum with the first term the 
imaginary part vanishes identically and the real part gives the same contribution two times. 
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Now the x-integration can be performed, and we have after transforming back to a general 
Lorentz system^ 



h = 



®P 2 y/(qu) 2 - m 2 

( -2pq + 2(pu){qu) 
\ y/ (qu) 2 — m 2 



{pu)m 2 — (pq)(qu) 
(qu) 2 — m 2 



In 



qu 



\J (qu) 2 — m? 



(D.5) 



The next integral is easier to calculate (use again f„ dr- 



~ I d*k5(k 2 



1 



gl a|fc«| _ i 



dr— 



1 



Now to I4: 



d 4 kS(k 2 



2 eP r - 1 



1 



■ 47T 



qu + -J (qu) 2 — m 2 

IT 2 7T 3 , , 

-2tt-— 2=-— (D.6) 



6/3 



3/3 2 



ki 



( e /3|H-l) (-2qk 
We proceed as in the calculation of I2 above and finally get 

^3 



■iO) 



h = 



1 



6/3 2 ^( qu ) 2 - w? 
There remains one integral to calculate: 



qu — y/ (qu) 2 — m 2 



d 4 kd(k 2 



qu + y/ (qu) 2 — m 2 
(ku) 2 



(e^lM-l) ( p k - iO)(qk - iO)' 



(D.7) 



(D. 



(D.9) 



Again we choose polar coordinates. The fco-integration is done with 8(k 2 ), and the re- 
integration decouples from the other variables and can be performed (using formulae as for 
the r-integration in I 2 ). Then we transform x — > —x, ip —* ip — ■n in the term with fco = — r, 
which leads to an integral . We transform its part J_ with ip — > <p + 2it to jj^. Now we 
can see that the term with fc = — r equals the complex conjugate of the term with k = r 
and thus the sum is two times its real part. We get 



2 6/5 2 



dx 



dp 



(Po - \p\x) (q - qi sin(¥>)\/T 



q 2 cos((^)Vl - x 2 - q 3 x) 
(D.10) 



The (^-integration now gives- 

f27T 

dip- 







gi sin(<^)vi — x 2 — (72 cos(y)vl — x 2 



q3x 



\/\q\ 2 x 2 



Thus we have 



h = 



&I3 2 



dx- 



(Po ~ IpIx)^/^ 2 ^ 2 Qoq3X + (m 2 + qf) 



2tt 

Zqoq?,x + (m 2 + ql) 
(D.ll) 



(D.12) 



Here we give some relations necessary to perform this transformation from the special system with 



u = (1,0), q 
\/ (qu) 2 — rri- 



(go, 0, 0, \q[) back to a general one: po — > p«, |p| — » J (pu) 2 



'la 



P3 



(pu)(qu)-pq 



qu, \q\ 



These relations can be verified straightforward. 



58 It has to be verified that the possible poles are integrable. This is done best by cutting them from 
the integration interval and then extend the result to this border values of the interval. In addition, using 
q 2 = m 2 > 0, \x\ < 1, etc. it can be shown that the radicand in the resulting formulae is always positive. 
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Having performed this integration (which can be done quite straightforward) , we transform 
the result back into a general Lorentz framep^] and thus have 



h =- 



In - 



m 4 ( qu 



(pu)(qu)-pq 
yj (pu) 2 — m 2 



12(3 2 v / (pq) 2 - to 4 [ \y/(pq) 2 -m*(qv 



(pu)(qu)-pq 
W (pu) 2 —m 2 



yj (pu) 2 —m- 



• (m 2 (pu) 2 — m 4 + (pq) 2 — (pq)(pu)(qu)) 



In 



/ y/ (pq) 2 — m i (qu + 



y/ (pu) 2 —m 2 

m 2 (pu) + (pq)(qu) \ 
m 2 (pu) — (pq)(qu) J 

(pu)(qu)-pq 



(m 2 (pu) 2 — m 4 + (pq) 2 — (pq)(pu)(qu)) 



V (p«) 2 



(D.13) 



m 4 ( gii + 
1 



(P")(9")-P'? 
-y/ (pu) 2 —m 2 

= (m 2 (pu) 



m 4 + (pq) 2 - (pq)(pu)(qu)) 



^/(puj 2 

m 2 (pu) - (pq)(qu) \ 
m 2 (pu) + (pq)(qu) J 



(m 2 (pu) 2 — m 4 + (pq) 2 — (pq)(pu)(qu)) 



Now we want to solve the linear system ( D.l ). We did that using Maple and get the general 
solution for A, G. Then we set q :=p — r\ with q 2 = (p — 77) 2 = m 2 and let it be simplified 
and factorised. Thus we get (since we will not need the solutions for B,C,D and G here, 
we just give the one for A): 



A 



8(pu) 2 I 3 m 2 



4m 4 h 



+ 4m 4 / 3 - ^(pufh + A(pufh + 8(pu) 2 m 2 I 2 



— 8(pufh(pi]) - 2(pu) 2 I 2 (pr 1 ) - Q(pufh(r,u) + &h(pu) 2 (prj) + 2m 2 I 2 (pi]) 
+ 6I 3 m 2 (pu)(rju) — Am 2 (pu) 2 I^rju) — 10m 2 (pu)l2(rju) + 8(pu)I 4 m 2 (pvi) + 
+ A(pu) 2 m 2 I^(prf) + 4/4771 (rju) — 4/5777 (prj) — &I 3 m 2 (prf) + 10(pu) 3 I 2 (r]u) 
+ 2(pu) 2 I 3 (r 1 u) 2 - 2(pu) 2 h(pr l ) 2 + 2m 2 I 2 (r 1 u) 2 + M b m 2 (pri) 2 - 

— 8(pu) 2 I 2 (i]u) 2 — 4(pu)m 2 15 (pi]) (i]u) + Am 2 (pu)Ii(riu) 2 — 6l4m 2 (pi])(i]u) - 
+ 10(pu) 2 I 4 (r)u)(pr/) - 4(pu)I 3 (pr])(r)u) + &(pu)I 2 (pr))(r)u) + 2I 3 (pi]) 2 - 

— 4(pu)h(pr/) 2 - I 2 (r]u) 2 (pr)) + 2h(pr]) 2 (rju) + 

+ 2(pu)I 2 (riuY + 2(pu)h(pri) 2 (r 1 u) - A(pu)h(riu) 2 (pij) - h(pv) 3 ) / 
I ( + 4?77 4 ( 777/) — 8m 2 (pri)(pu) 2 + A(pu) 4 (prj) + 2m 4 (riu) 2 — 

— 2(pu) 2 m 2 (r]u) 2 — 4: r m 2 (p r q) 2 + 4 (pi]) 2 (pu) 2 + 2(pu)m 2 (i]u) 3 — 

— m 2 (pi])(i]u) 2 — 2(pu) 2 (i]u) 2 (prj) + (prj) 3 



(D.14) 

Having a close look at the solution we see what is the problem of setting p = q <-> 77 = 
too early. The determinant of the linear system and the denominator in the solutions, resp., 
is zero for this values. Therefore we must not set 77 = but expand the various quantities 
involved in 'powers of 77'. Then we will see that all possibly divergent terms cancel and we 

59 See footnote ^. 
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are left with finite values for rj — > 0. This expansion is rather delicate and we have to pay 
attention not to throw away any term too early or even to forget some at all. A crucial thing 
with the linear system considered is the fact that it remains defined even if we set r/ = 
at the very beginning; we then can reduce the system to four equations and solve it. And 
we are in the lucky case that we only need to know A for which we get the same result this 
way. But this is not true in general and here it is also wrong for the other variables B, etc. 

To do these things we introduce some notation: set ff =: e\vecq e with egl, e<l, ff e G 
R, ff 2 = 1, so that the limit in M, rj — > 0, can be discussed as an ordinary limit in R, e — > 0. 
Pay attention to the fact that 770 depends on ff through q 2 — (p — rj) 2 = m 2 . We have 

? /o = Po — Poa/I — — Thus 7/0 is not simply proportional to e! This has to be 

V Po Po 

investigated for q 7^ and q = separately. But considering the first case it will be clear 
that the second follows by simply setting q = 0. 

Since the more direct but in general wrong method of setting p — q at the beginning 
leads to the same result we will not present here a detailed derivation of the calculations. We 
just give the general direction: All quantities involved in the numerator and denominator 
of A have to be expanded in e up to a suitable order. Especially the integrals J2> ••• ) Lj have 
to be expanded. This is a tedious but strightforward task. We emphasize especially the 
following: The structure of A is 



A _ L + L 1 e + L 2 e 2 + L 3 e 3 + ■ ■ ■ _ 
he 2 + hei + he± 

1 1 „ ( D - 15 ) 

= ^- = = (Lo + Lie + L 2 e 2 + L 3 e 3 + •••)■ 

he 2 (i + h e + l fe 2 ) + --- 

Now we expand the numerator, identify which term of A has to be known up to which 
order and calculate it straightforward. This leads to 



^- = + ^ - ))• (D - i6) 

E The Cross Section with General Mixed States 



The cross section involving the initial state ( 3.109|) and the corresponding mixture for the 



final state reads as follows ('0','1' and '2' in the matrix elements stand for contributions 
with zero, one or two bremsstrahlungs photons involved, the sums run from 1 to 4 and the 
hi, ...-integrations over ki,... < u>q, 'MS' indicates the mass shell and 'MSi' a small part 
of the mass shell centered around p. b + in the final state is always understood with the 
momentum p, in the initial one with momentum q.): 
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/ d*p [ d A q 

JMSi J MS 



+ |(6+|0+l + 2|6+)| 2 Co (g) + 



+ ]T J d 3 k 1 \(b + \0+l + 2\b+a+(k 1 ))\ 2 c 1 (q,k 1 ) t + 

i 

+ J2 I d 3 k 1 \(b + \0+l + 2\b+h+(k 1 ))\ 2 £ 1 (q,k 1 ) l + 
+ Y, / rf3fc i j d 3 k 2 \(b + \0+l + 2\b + a+{k 1 )a+{k 2 )}\ 2 c 2 (q, fei,fe)y + 
+ S/ rf3fcl / rf 3 fe|(6 + |0+l + 2|6+a+(fc 1 )a+(fc 2 ))| 2 2 2 ( g , fei,fe)y + 
+ Y,J rf3fc i / rf 3 fe|(6 + |0+l + 2|6+a+(fc 1 )a+(fc 2 ))| 2 2 2 ( g , k u k 2 )ij + 

+ E / ^i|(^ai)|0 + l + 2|6+)| 2 Co ( g ) + 

+ I Sh I d 3 k 1 \(b+a+(l 1 )\0 + l + 2\b+a+(k 1 ))\ 2 c 1 (q,k 1 ) l + 

+ EE / d * h I d 3 k 1 \(b+a+(l 1 )\0 + l + 2\b+&+(k 1 ))\ 2 c 1 (q,k 1 ) l + 

+ d3/ i / rf3fc i / d 3 k 2 \(b + a+(h)\0 + l + 2\b + a+(k 1 )a+(k 2 ))\ 2 c 2 (q,k 1 ,k 2 ) t3 + 

+ I d * h I d3kl I d 3 k 2 \(b+a+(l 1 )\0 + l + 2\b+a+(k 1 )&+(k 2 ))\ 2 £ 2 (q,k 1 ,k 2 ) l3 + 

m i,j J J J 

+ f d * h / d ' 3fcl / rf 3 fc 2 |(fo + a+ai)|0 + l + 2|6+a+(fc 1 )5+(fc 2 ))| 2 £ 2 (g,fc 1 ,A :2 ) 4 , + 
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+ E / ^i|(6 + a+(/i)|0 + l + 2|6+)| 2 Co ( (Z ) + 

m J 

+ I d * h / d 3 k 1 \(b+&+(l 1 )\0 + l + 2\b+a+(k 1 ))\ 2 c 1 (q,k 1 ) l + 

m i 

+ E]T I d\ [ d 3 k 1 \(b+&+(l 1 )\0 + l + 2\b+&+(k 1 ))\ 2 d 1 (q,k 1 ) l + 
mi'' J 

+ I d * h I d3kl [ d 3 k 2 \(b+a+(l 1 )\0 + l + 2\b+a+(k 1 )a+(k 2 ))\ 2 c 2 (q,k 1 ,k 2 ) l3 + 

m i,j J J J 

+ d * h [ d3kl /rf 3 fc2|(& + a+ai)|0 + l + 2|6+a+(fc 1 )5+(fc 2 ))| 2 £ 2 (g,fc 1 ,A :2 ) y - + 

m i,j J J J 

+ d3/l / d3kl [ rf 3 fc 2 |(& + a+ai)|o + i + 2|6+a+(fc 1 ) ( i+(fc 2 ))| 2 £ 2 (g,fc 1 ,A :2 ) 4 , + 

m i,j J J J 

+ / d%\{b+at n (h)a+{l 2 )\() + l + 2\b+)\ 2 c {q) + 

+ / / <? 1 * / d 3 k 1 \(b + a+(l 1 )a+(l 2 )\0 + l + 2\b+a+(k 1 ))\ 2 c 1 (q,k 1 ) l + 

m,n i J J J 

+ I Sh I ^ 2 I d 3 k 1 \(b+a+(l 1 a+(l 2 ))\0 + l + 2\b+a+(k 1 ))\ 2 £ 1 (q,k 1 )i + 

m,n i J J J 

d3h j d3/2 / rf3fci / d3fc2 

m,n ijj 

\(b + a+(l 1 )a+(l 2 )\0 + l + 2\b+a+(k 1 )a+(k 2 ))\ 2 c 2 (q,k 1 ,k 2 ) l ,+ 

! d3h I d3/2 / rf3fci / d3fc 2 

m,n i,j J J J J 

|(6+a+(Z 1 )a+(Z 2 )|0 + 1 + 2\b+a+ (k 1 )a+ (k 2 ))\ 2 £ 2 (q, k u k 2 ) l3 + 
+ d * h I rf3/ 2 / d 3 ki f d 3 k 2 

m,n i,j J J J J 

\(b+a+(l 1 )a+(l 2 )\0 + l + 2\b+a+(k 1 )af(k 2 ))\%(q,k 1 ,k 2 ) ij + 
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+ S/ dHl J d 3 h\(b + a+(l 1 )&+(l 2 )\0 + l + 2\b + )\ 2 c (q) + 

+ I dHl I I ^k 1 \{b+al l {h)h+{l 2 )\Q+l + 2\b+a+{k 1 ))\\ 1 {qMh 

m,n i •* 

+ l dHl I d3h I rf 3 ^i|(fe + «^(^i^(^))|0 + 1 + 2|b+a+(fc 1 ))| 2 5 1 ( 9 ,fc 1 )i 

m,n i J 

m,n ij J J J J 

|<6+a+(ii)a+(i 2 )|0 + 1 + 2|6+a+(fc 1 )a+(fc 2 ))| 2 C2 ( g , k u + 

m,n i.j J J J J 

|<6 + o+(ii)5+(i 2 )|0 + 1 + 2|6+a+(fc 1 )a+(fc 2 )> | 2 2 2 (g, k u k 2 ) ij + 

/ d3?i / rf3 ^/ d3fci / d3fe 

(6 + o+(ii)5+(i 2 )|0 + 1 + 2|6+a+(fc 1 )a+(fc 2 ))| 2 2 2 ( g , fci,fe)y + 



+ 



+ 



+ J2 I dH ^ /rf 3 ^|(b + a+(Zi)a+(? 2 )|0 + 1 + 2|fo+)| 2 c (g) + 

m,n 

I d3/ i / dH ^ f rf 3 fci|(6 + «™aiKta2)io + i + 2|fe+ a +(fci)>l 2 ci(?,fci) 4 + 

m,n i J J J 

+ I dHl f dH2 [ d 3 k 1 \(b + a+(ha+(l 2 ))\0 + l + 2\b + af(k 1 ))\ 2 c 1 (q,k 1 ) i + 

m,n i J J J 

[ d3h / d3?2 / d3ki / rf3fc2 

|(6+a+(« 1 )a+(/ 2 )|0 + l + 2|6+a+( fc i)4( fc 2))| 2 C2(9,fci,fc 2 ) ij + 

[ d3h / d3?2 / d3fci / rf3fc 2 

m,n j,j J J J J 

|(6+a+(Z 1 )a+(Z 2 )|0 + 1 + 2|6+a+(fc 1 )a+(fc 2 ))| 2 £ 2 (g,fc 1 ,fc 2 ), J + 

m,n i,j 

|(6+a+(« 1 )a+(/ 2 )|0 + l + 2|6+a l +(fc 1 )at(fc 2 ))| 2 g 2 ( g ,fci,fc 2 ) ij 



(E.l) 
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Now we collect the contributions that are not zero and of order not higher than four. Thus 
we are left with: 

/ d A p f d\\ 

JMSi JMS I 

+ \(b + \0\b + )\ 2 [c (q) + Y< J d 3 k lCl {qM)i + Y, J d 3 k 1 £ 1 (q,k 1 ) t + 

i i 

+ d3kl j d3fc 2C2(<Z,fcl,fc 2 ) tJ ■ +J2 j d3fcl / d3fc 25 2 ('Z,fcl,fc2) lJ ■ + 

i,3 i,j 

+ ^ J d 3 ki J d 3 k 2 C2(q,h,k2)ij^ + 

+ E / d3fcl [ + l(& + |l|& + a+(fci))| 2 ci(9, h)i + \(b+\l\b + d+(k 1 ))\ 2 c 1 (q, kt)i + 

+ \(b+a+(k 1 )\l\b+)\ 2 c Q (q)+Y / I d 3 k 2 \(b+\l\b+a+(k 1 ))\ 2 c 2 (q,k 1 ,k 2 ) l3 + 

i J 

+ W d 3 fc 2 |(6+|l|fe+a+(fc 1 ))| 2 5 2 (q,fc 1 ,fc 2 ) y + |(fe+a+(fc 1 )|l|6+>| 2 C0 (g) + 
i J 

+ f d 3 k 2 \(b+\l\b+a+(k 1 ))\ 2 £ 2 (q, fci, k 2 )ij + 

3 J 

+ [ d 3 k 2 \(b+\l\b+a+(k 1 ))\%(q,k 1 ,k 2 ) l3 + 

3 J 

+ W d 3 k 2 \(b+a+(k 1 )\l\b+)\ 2 c 1 (q,k 2 ) J + 

3 J 

+ J2f d 3 k 2 \{b + ~a+{k 1 )\l\b + )\ 2 c 1 {q,k 2 ) 3 + 

3 J 

+ £ I d 3 k 2 \(b+af(k 1 )\l\b+)\ 2 c 1 (q,k 2 ) j + 

3 J 

+ d 3 k 2 \(b + &+(k 1 )\l\b+)\ 2 £ 1 (q,k 2 ) J \ 



+ 
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EE / d3/ i / d 3 fci2Re[(fe + |0|6+)5 m ^(/ 1 -fc 1 )(6 + a+(Z 1 )|2|6+a+(fci)>*ci((?,fci) ! 

mi 

EE / d3/l / rf 3 ^i2Re[(fo + |0|6 + )5 m ^(/ 1 -fc 1 )(6 + a+(/ 1 )|2|6 + 5+(fc 1 ))*c 1 (( 7 ,fc 1 ) ! 

mi'' J 

EE / d3/i / dH2 1 d3ki I d3fc2 ' 

m,n i,j J J J J 

■ 2Re[(b + \0\b+)S mi S nj S(h - h)S(l 2 - k 2 ) ■ 

■ (fe + a+ai)|2|fo+a+(fc 1 ))*^ 2 -fc 2 )5„,c 2 (q,fc 1 ,fc 2 ). i ,] + 
EE / d3/l / ^ 2 J d3kl y"d 3 fc 2 2Re[(6 + |0|&+) ( 5 m4( 5„^(/i-fc 1 )(5(/ 2 -fc 2 )- 

■ {<6+o+(ii)|2|6+o+(fci)>*i«a 2 - fc 2 )5 nj + 
+ (6+a+G 2 )|2|6+at(A ;2 ))*^(/i-fc 1 )5 roi }c 2 (g,fc 1 ,fc 2 ) i 

EE/ d3/i / d3/2 / d3fci / d3 



m,n z, j 



j3 fc 2 • 



• 2Re 



(b + \0\b + )6 mi 5 nj 5(l 1 -k 1 )8(l 2 -k 2 )- 

■ (b + a+{h)\2\b+a+(ki))*8(l2 - h 2 )5 nj ^(q, fei, fc 2 )« 



(E.2) 



The products of identical ^-distributions showing up in some terms of this formula indicate 
that we should have done this calculation using wave packets instead of sharp momentum 
states in the description of the asymptotic states. Since the products of zero- and two-BS 



contributions play no role in our discussion since they are finite (cf. Sec. 3.4.3), we will not 
do this here explicitly, we just state that the analogous calculation involving wave packets 
leads to the same result with only one of the identical S's in each term. Thus we throw away 
one of them in the expression above. 
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